Microfabricated Atomic Vapor Cells with Multi-Optical Channels Based on an Innovative Inner-Sidewall Molding Process
Mingzhi Yu
,
Yao Chen
,
Yongliang Wang
,
Xiangguang Han
,
Guoxi Luo
,
Libo Zhao
,
Yanbin Wang
,
Yintao Ma
,
Shun Lu
,
Ping Yang
,
Qijing Lin
,
Kaifei Wang
,
Zhuangde Jiang
Existing microfabricated atomic vapor cells have only one optical channel, which is insufficient for supporting the multiple orthogonal beams required by atomic devices. In this study, we present a novel wafer-level manufacturing process for fabricating multi-optical-channel atomic vapor cells and an innovative method for batch processing the inner sidewalls of millimeter glass holes to meet optical channel requirements. Surface characterization and transmittance tests demonstrate that the processed inner sidewalls satisfy the criteria for an optical channel. In addition, the construction of an integrated processing platform enables multilayer non-isothermal anode bonding, the filling of inert gases, and the recovery and recycling of noble gases. Measurements of the absorption spectra and free-induction decay signals of xenon-129 (129Xe) and xenon-131 (131Xe) under different pump-probe schemes demonstrate the suitability of our vapor cell for use in atomic devices including atomic gyroscopes, dual-beam atomic magnetometers, and other optical/atomic devices. The proposed micromolding technology has broad application prospects in the field of optical-device processing.
With the advent of quantum technology, precision measurement devices based on quantum phenomena have been extensively studied [1], [2], [3]. Examples of such instruments include superconducting quantum interference devices (SQUIDs) [4], [5] and optically pumped magnetometers (OPMs) [6], [7]. Alkali metal atomic vapor cells, which are sealed with alkali metals and a specific volume of noble gas, are the core components of quantum devices such as atomic clocks [8], [9], [10], atomic magnetometers [11], [12], [13], and atomic gyroscopes [14], [15].
Traditional alkali metal atomic vapor cells are manufactured using the glass-blowing method [16], [17], [18], which is expensive, poses difficulties in batch production, and does not guarantee consistency. In the glass-blowing method, a stem of a few millimeters is attached to the vapor cell, increasing the quadrupole frequency shift when the atomic nuclear spin I ≥ 1 [19], [20]. This occurs with xenon-131 (131Xe), hydrargyrum-201 (201Hg), krypton-83 (83Kr), and neon-21 (21Ne), causing the primary systematic errors that limit the bias stability in nuclear magnetic resonance (NMR) gyroscopes [21]. Moreover, traditionally fabricated alkali metal atomic vapor cells are incompatible with microelectromechanical systems (MEMSs) and cannot be applied to chip-scale atomic devices. Wafer-level fabrication makes it possible to achieve MEMS alkali metal vapor cells with sufficient miniaturization, low power consumption, and low manufacturing costs [22]. In addition, in microfabricated MEMS vapor cells, having no stem can reduce the quadrupole frequency shift [20], [23], [24].
The MEMS alkali metal vapor cell was first proposed by Kitching et al. [25]. During the early development of the MEMS vapor cell [26], [27], [28], [29], [30], [31], it had only two transparency windows. In spin-exchange relaxation-free atomic magnetometers or atomic spin gyroscopes, multi-beam configurations are required to achieve better performance [11], [32], [33]. Therefore, developing multi-beam MEMS vapor cells for chip-scale atomic devices is necessary.
Based on a wafer-level glass-blowing procedure, Eklund et al. [34], [35] developed an MEMS spherical alkali metal atomic vapor cell. Although a multi-beam optical channel was enabled, the laser beam diverged because of the nonuniform wall thickness of the spherical vapor cell. Divergence of the laser beam may reduce spin polarization. Mirrors built into the cavities offer an alternative solution. Perez et al. [36], [37] fabricated a Bragg reflector on the inclined sidewall of a cavity using bulk silicon (Si) micromachining and multilayer plasma-enhanced chemical vapor deposition. The laser beam flowing through the cavity was redirected to the light source plane. However, the uneven thickness of the deposited Bragg reflector resulted in inconsistent reflectivity and higher angle accuracy requirements for laser incidence. Han et al. [38] built a vapor cell with an aluminum reflection sidewall; however, the aforementioned issues persisted, despite the relative simplicity of the production process. Moreover, in an NMR gyroscope, a vapor cell with reflective sidewalls disrupts the symmetry of the vapor cell, resulting in a frequency shift of 131Xe caused by nuclear quadrupole interaction. In summary, although researchers have investigated a number of solutions, the batch fabrication of atomic vapor cells with multiple optical channels based on the MEMS process still has the drawbacks of complex fabrication, irregular chamber shapes, the introduction of metal impurities, and high requirements for external laser incidence accuracy.
In this study, we develop a side-molding process to fabricate glass cavities with high-quality inner sidewalls in order to facilitate the mass production of MEMS multi-optical-channel alkali metal atomic vapor cells. In this process, a mold with a smooth surface is plugged into the glass cavity, and the mold and inner sidewalls squeeze each other at a specific viscoelastic temperature, such that the rough sidewalls of the glass cavity are deformed into an optical plane. Using this technique, an atomic vapor cell with multiple optical channels is fabricated, and its effectiveness is demonstrated through a series of optical tests. Platforms for testing the vapor cell are also developed to acquire the free-induction decay (FID) signal and the relaxation times of xenon-129 (129Xe) and 131Xe.
2. Methods
2.1. Structure design
The structure of the MEMS vapor cell is shown in Fig. 1. It is composed of top glass, Si, middle glass, Si, and bottom glass, which are connected using anodic bonding technology. The thicknesses of the top and bottom glasses are 500 μm, the thickness of the middle glass is 2 mm, and the Si is p-type (100)-oriented high-resistance monocrystalline Si with a thickness of 500 µm. SCHOTT BOROFLOAT®33 (BF33) glass is used.
The total size of the vapor cell is 5 mm × 5 mm × 4 mm, and the internal cavity size is 3 mm × 3 mm × 3 mm. The cavities contain rubidium (Rb), nitrogen gas (N2), 129Xe, and 131Xe. After precision-side molding, the vapor cell has three optical channels through which three orthogonal laser beams pass.
2.2. Inner sidewall molding
Optical glass molding technology is based on continuous and reversible thermal processes, with glass properties varying from the molten to the solid state [39]. Near the glass transition temperature, the glass and aluminum oxide (Al2O3) molds were heated and pressurized under oxygen-free conditions, and the glass was molded into optical components. Compared with the traditional molding process, which requires pressure to deform the glass, we used the thermal expansion of the mold and glass to cause them to squeeze each other. Subsequently, the high-precision surface of the Al2O3 mold was copied onto a glass surface. Thus, the glass surface satisfies the optical interface requirements.
Compared with the traditional molding process, this process has several advantages: First, it can remove the inner surface morphology, whereas traditional molding processes cannot. Second, the glass deformation is minimal (only a few tens of microns) when the thermal expansion of the material is used during extrusion molding; thus, the internal stresses generated in the glass are miniscule, and the optical properties of the glass, such as the refractive index, are maintained. Third, no molding pressure is applied, indicating that mold damage is minor, which is conducive to mold recycling and results in reduced cost.
Fig. 2 shows the fabrication process of the MEMS multi-optical-channel alkali metal atomic vapor cells. The treatment steps for the inner sidewall are shown in Figs. 2(a)-(f). Fig. 2(a) depicts the laser drilling of a glass wafer with a drilling size of 3 mm × 3 mm and a depth of 2 mm. Because the efficiency of the laser process for fabricating deep holes is much greater than that of wet or dry etching using MEMS technology, the laser process was used to fabricate deeper cavities. The parameters used for the laser drilling are listed in Table 1.
Fig. 3 shows the pictures of MEMS atomic vapor cell wafers after different process treatments. A glass wafer fabricated by means of laser drilling is shown in Fig. 3(a). After testing, the precision of the drilling size was greater than 5 μm, the surface roughness of the inner sidewall was less than 5 μm, and the inner sidewall steepness was 90.0° ± 0.2°. Thus, this treatment method has high processing efficiency, high dimensional accuracy, and high angular accuracy. Even though the edges of the holes have defects that are several tens of microns wide, they do not have any bumps; therefore, anodic bonding is unaffected.
Fig. 2(b) shows the etching of Si holes using inductively coupled plasma (ICP) etching. The size of the Si hole is slightly larger (by several microns) than that of the glass hole to prevent extrusion of the Si and α-Al2O3 molds during molding. Fig. 2(c) depicts the first anodic bonding of the etched Si wafers and the laser-drilled glass wafer. The bonding voltage was 800 V, and the temperature was 400 °C.
Fig. 2(d) shows that a layer of metal was sputtered on the peripheral side of the previously bonded wafer as a side lead to ensure proper voltage during anode bonding between the middle glass layer and the lower Si layer. A Si wafer was bonded to the other glass side, with the same bonding parameters as those previously described. The bonded Si-glass-Si wafer is shown in Fig. 3(b), and its cross-sectional view is shown in Fig. 3(c).
Fig. 2(e) shows the mold preparation process. First, the bonded Si-glass-Si wafer was etched in an aqueous solution of 5% HCl, 30% HF, and 10% 3-indolepropionic acid for 15 min to remove the burrs and smooth the inner surfaces of the laser-drilled holes. The bonded wafer was cleaned with deionized water and dried with N2. The α-Al2O3 molds were placed in the holes, as shown in Fig. 3(d). Wet etching improved the consistency of the sidewall roughness such that the mold would not be scratched by the glass sidewall when it was positioned, ensuring the reuse of the mold.
The dimensional accuracy of the α-Al2O3 molds was ± 1 μm, the surface shape accuracy peak to valley (PV) was greater than λ/4 (λ = 632.8 nm, the surface shape accuracy PV was measured using a laser interferometer whose light source was a He-Ne laser with a wavelength of 632.8 nm), and the surface roughness was less than 10 nm. The reasons for choosing α-Al2O3 as the mold are as follows: This material has good high-temperature resistance, wear resistance, and chemical stability; high hardness and elastic modulus; and very low glass reactivity and glass bonding propensity [40]. Compared with the tungsten carbide mold used in traditional molding, Al2O3 has a large thermal expansion coefficient, which is beneficial for molding during thermal expansion.
The viscosity-temperature curve of BF33 glass can be obtained from Schott’s official website [41], and the main viscosity points are listed in Table 2.
The distance between the mold and the inner sidewalls of the glass holes should be appropriate. A mold that is too large is not conducive to demolding, whereas a mold that is too small will lack sufficient mold pressure. During heating, the distance between the mold and the inner sidewalls of the glass holes (D1) changes due to thermal expansion, as given by the following equation:
where Tm is the molding temperature, which is generally between transition point Tg and softening point Ts. From Table 2, the range of Tm is 873-1123 K, while α1 = (7.85 ± 0.02) × 10−6 K−1 is the average linear expansion coefficient of α-Al2O3. L1 = (3000 ± 1) μm is the size of the α-Al2O3 molds, and α2 is the thermal expansion coefficient of BF33 glass. The specific value of α2 can be obtained from the official Schott website [41], where L2 = 2 mm is the total thickness of the middle glass, due to the nonlinear deformation of the glass above Tg = 798 K. In this configuration, D1 = 29.0-37.4 μm. Therefore, the mold size L1 is several microns smaller than the cavity size of the vapor cell, so molding can still be achieved. After the measurement, the actual size of the hole was (3000 ± 5) μm. Therefore, we fabricated five types of molds with side lengths of 2990, 2985, 2980, 2975, and 2970 μm to determine the appropriate mold size.
Fig. 2(f) shows the molding steps. The wafer bonded to the molds was placed in a furnace and molded under vacuum. As shown in Figs. 2(c) and (d), Si wafers were bonded in advance on both sides of the middle glass, effectively preventing nonlinear deformation of the middle glass at high temperatures, to ensure the flatness of the upper and lower sides of the bonded wafer. Thus, the mold pressure did not affect the subsequent anodic bonding.
The molding experiment was performed using the temperature curve shown in Fig. 4. First, the temperature was raised to the molding temperature of Tm at a rate of 20 °C·min−1. Then, the inner sidewalls of the middle glass holes were fully molded by remaining at Tm for 15 min. Afterward, the sample was cooled to the annealing temperature Ta at 5 °C·min−1 and held at that temperature for 30 min to eliminate the internal stresses of the glass. Finally, the sample was cooled to room temperature at 10 °C·min−1. Two temperature cycles were performed in each molding experiment. The mold and glass naturally separated, making demolding easy, due to the shrinkage of the mold and the glass after cooling. Fig. 3(e) shows the cross section of a Si-glass-Si wafer after the inner sidewall molding.
3. Results
3.1. Optimization of the Al2O3 mold size and molding temperature
As shown in Fig. 5, we studied the average surface roughness Sa of the inner sidewall of the middle glass hole at various mold sizes and temperatures to determine the appropriate optimal mold size and temperature Tm. Because the longitudinal coordinates were large, it was not possible to mark the roughness error for multiple samples with the same mold size and temperature. The average roughness values are shown in Fig. 5. The roughness of the inner sidewalls of the molded glass holes gradually decreased as the mold size and temperature increased. However, a mold that is excessively large can easily cause adhesion between the mold and the glass, resulting in molding failure. The 2985 μm mold adhered to the glass at temperatures greater than 775 °C, resulting in the inability to obtain roughness data. The 2990 μm mold could not be used, because it could not be smoothly placed in the glass holes.
In addition, adhesion was inevitable when the molding temperature exceeded the softening point Ts. When the mold size was too small, the molding process was insufficient, indicating that Sa did not meet the application requirements. In summary, we chose a 2980 μm mold for molding at 790 °C; this achieved the best roughness, because it eliminated adhesion between the mold and the glass hole, which will be discussed in detail in the next section.
3.2. Surface characterization of the molded inner sidewalls of the middle glass holes
Fig. 6 shows the inner sidewall surface morphologies of the glass holes before and after molding using laser scanning confocal microscopy. The size of the sampling area is 640 μm × 640 μm. The inner sidewall was ground glass before molding, as shown in Figs. 6(a) and (b). As shown in Table 3, the average surface roughness Sa of the sidewall surface was 5.653 μm, the average surface shape accuracy PV was approximately 67.304 μm, and the average kurtosis Sku was 6.236. The parameter Sku was used to determine the sharpness of the rough shape [42], as shown in Fig. 7. Many sharp peaks (i.e., burrs) were observed on the surface, as shown in Fig. 6(b).
As shown in Table 3, after wet etching, Sa was 3.630 μm, the PV was 38.379 μm, and Sku was 3.154. The sharp part of the surface burr was removed, and the surface roughness was relatively mild. Thus, wet etching can help improve the quality of molding and prevent burr scratching.
Figs. 6(c) and (d) show the surface morphology and three-dimensional (3D) topography after final molding. As shown in Table 3, Sa was 22 nm, the PV was 283 nm (less than λ/2), and Sku was 2.506.
An infrared-visible-ultraviolet spectrophotometer was used to measure the transmittance of the molded inner sidewalls of the glass holes in the near-infrared band, as shown in Fig. 8. The red curve at the top represents the transmittance of the polished BF33 glass, which has an average transmittance of approximately 94%. The five samples shown in Fig. 8 were collected from various parts of the wafer, ranging from the edge to the center. Their transmittance was not significantly correlated with the distribution location and fluctuated between 87% and 91% after measuring different batches of wafers, indicating that the molded inner sidewalls satisfied the requirements of the optical channel of the vapor cell.
3.3. Multilayer anodic bonding
As shown in Fig. 2(g), the peripheral edges of the two Si layers of the molded wafer were connected via magnetron sputtering with a metal. Then, the bottom glass was connected to the molded wafer via anodic bonding. Fig. 2(h) shows how the molten alkali metal was injected into the vapor cell using a micropipette in a N2 environment.
As shown in Fig. 2(i), the top glass and the molded wafer with the alkali metal were placed in the bonding cavity. The bonding cavity was vacuumed to 1 × 10−5 mbar (1 mbar = 100 Pa), and 5 Torr (1 Torr = 133.3224 Pa) 129Xe, 15 Torr 131Xe, and 600 Torr N2 were added to the bonding cavity. Subsequently, the upper plate of the anodic bonding equipment was pressed down to enable nonisothermal bonding. The temperatures of the upper and lower plates were 280 and 200 °C, respectively. The bonding voltage was 1200 V, and the duration was 120 min. The purpose of this step was to maintain the alkali metal at the bottom glass based on the temperature gradient to prevent it from affecting the bonding surfaces, thereby achieving a preliminary sealing of the vapor cell. Isothermal anodic bonding was also conducted. The temperatures of the upper and lower plates were both 280 °C, the bonding voltage was 1200 V, and the duration was 60 min to improve bonding strength and hermeticity.
A manufacturing platform with vapor cell bonding, gas filling, and noble gas recovery capabilities is shown in Fig. 9. The bonding cavity was placed on a platform in an inert gas glovebox. During the fabrication of the vapor cell, 129Xe and 131Xe, which are expensive, were selected as noble gases, and the volume of the bonding cavity was relatively large. We calculated that only about one thousandth of the Xe gas that filled the bonding chamber was sealed in the vapor cell each time. If noble gases are used only once, the processing costs will be high; therefore, it was necessary to design a gas-recovery system to reduce the manufacturing cost. In the recovery system, liquid nitrogen was used to cool the recovered gas mixture, the cheap nitrogen in the gas was liquefied, and the valuable Xe gas was stored in a cylinder for subsequent use.
A resin-bonded blade (i.e., a soft blade) with a 2000 mesh diamond surface was used for the wafer dicing. A high rotational speed and a low feed rate were used to cut the wafers with a good polishing effect on the cut surface. The surface roughness (Ra = (28 ± 5) nm) and face shape accuracy (PV = (180 ± 60) nm) were measured on the side of the cut wafer using a laser confocal microscope. The sample satisfied optical port quality requirements without additional polishing. Therefore, high efficiency of batch processing was ensured in this study.
The finished wafer containing Rb, 129Xe, 131Xe, and N2 is shown in Fig. 10(a). There were 236 vapor cells in the 4 in (1 in = 2.54 cm) wafer. A local magnification of the wafer is shown in Fig. 10(b), with the Rb clearly visible. As shown in Fig. 10(c), the segmented alkali metal atomic vapor cell has three optical channels, numbered 1, 2, and 3. The chip size is 5 mm × 5 mm × 4 mm, and the internal cavity size is 3 mm × 3 mm × 3 mm, meeting the requirements for the miniaturization and batch manufacturing of quantum devices such as atomic magnetometers and atomic gyroscopes with multiple laser channels.
3.4. Leakage rate
First, the hermeticity of the vapor cell was evaluated based on the leakage rate. The following method based on the MILSTD-883E standard method 1014.10 was used to measure the leakage rate: The vapor cell was placed in a sealed chamber, which was pressurized with helium (He) to pressure of exposure (PE) = 5.17 × 105 Pa and maintained at this pressure for T1 = 7200 s while keeping the temperature constant at 25 °C. Subsequently, the vapor cell was removed from the chamber and purged with N2 to remove the surface He. Next, the vapor cell was quickly placed in a helium mass spectrometer (ASM310; PFEIFFER VACUUM, Germany) to measure the amount of diffused He. The volume of the vapor cell was approximately V = 27 mm3, and the leakage rate of the vapor cell was calculated. The measured leak rate detected by the helium mass spectrometer leak detector was R1 = 4.4 × 10−10 Pa·m3·s−1, and the equivalent standard leak rate was L = 2.28 × 10−10 Pa·m3·s−1, as determined using Eq. (2).
where R1 is the measured leak rate of the tracer gas (He) through the leak, L is the equivalent standard leak rate, P0 = 1.01 × 105 Pa is the atmospheric pressure, MA = 28.7 g is the molecular weight of air, M = 4 g is the molecular weight of He, T1 = 7200 s is the time of exposure to PE, T2 = 1800 s is the dwell time between the release of pressure and leak detection, and V = 27 mm3 is the internal volume of the vapor cell.
The measured leakage rate was significantly higher than the actual leakage rates in air and nitrogen during real operation due to He permeation of the BF33 glass substrate. During the actual operation, the lifetime of the atomic gas chamber must also be assessed. There are several ways to mitigate He permeation, including switching to aluminosilicate glass instead of borosilicate glass [43]. However, this method is expensive. Other options include sputtering an alumina coating on the chamber walls [44] or using gallium phosphide as a sealing material for anode bonding [45]. These methods can significantly prolong the lifetime of atomic vapor cells.
3.5. Absorption spectra and FID signal
A schematic of the test system is shown in Fig. 11(a). This system was used to measure the absorption spectrum of the vapor cell and the FID signals of the different pumping detection schemes. This study aimed to evaluate the effects of the inner sidewalls of molded glass holes on the performance of a vapor cell.
Two non-magnetic heating coils were attached to both sides of the vapor cell in order to heat it to a specified temperature. During the measurement, the vapor cell was placed in a magnetic shielding barrel composed of four layers of permalloy and one layer of manganese-zinc (Mn-Zn) ferrite. As shown in Fig. 11(b), the magnetic fields inside the shields were controlled by three-axis magnetic coils with a set of solenoid coils for the longitudinal field and two sets of cosine-theta coils for the transverse field. In addition, a 795 nm circularly polarized light beam was passed through the vapor cell in the Z-direction to pump the Rb atoms. The Xe nuclear spins were hyperpolarized via spin-exchange optical pumping with Rb atomic spins. Another 780 nm linearly polarized light beam was passed through the vapor cell in the Y-direction to detect the atomic spin signal. Performing experiments using a free-space laser beam is difficult because the vapor cell is small. Therefore, a polarization-maintaining fiber is used to guide the pump and probe lasers into the vapor cell.
The vapor cell was heated to 393 K to obtain the absorption spectrum. After the temperature of the vapor cell was stabilized, the frequency of the laser was adjusted to set the wavelengths of the pump and probe lasers to approximately 795 and 780 nm, respectively. The output voltages of the photodiode at different detuning frequencies were recorded to obtain the absorption curves, as shown in Fig. 12.
The width of the absorption curve is mainly due to the natural broadening , Doppler broadening , and pressure broadening [46]. The natural broadening of the absorption spectral lines can be calculated using Eq. (3).
where τ is the natural lifetime of the excited state of an alkali metal atom, which is approximately 25-35 ns. The natural linewidth for the transition of the D1 and D2 lines of alkali metal atoms is 4-6 MHz. The effect of natural broadening on the absorption spectrum is negligible compared with a pressure broadening of more than 10 GHz.
The Gaussian function of the atomic absorption spectrum under the Doppler effect is given as follows:
where v is the pump laser detuning frequency, v0 = 377 THz is the resonance frequency of Rb.
where kB is the Boltzmann constant, T = 393 K is the temperature, m is the mass of Rb, and c is the speed of light. After the calculation, = 414 MHz.
The absorption spectrum under pressure broadening conforms to the Lorentz line distribution, which is expressed as follows:
where a and b are fitting coefficients.
As shown in Eq. (7), the absorption data should be fitted using the Voigt function, which is a convolution of a Gaussian function and a Lorentzian function. The Gaussian function accounts for the Doppler broadening, whereas the Lorentzian function estimates the pressure broadening caused by the collision of buffer gases [47].
where v′ is the integral variable of the convolution, denoting all possible positions in the range (0, ∞).
In this study, the buffer gas pressure was greater than 0.6 atm (1 atm = 1.01325 × 105 Pa); thus, the influence of Doppler broadening on the absorption spectrum can be ignored [48], [49]. For each Rb isotope, the ground-state hyperfine splitting is greater than 3 GHz. The sum of the four Lorentzian quantities in the form of Eq. (6) was used to fit the absorption spectral data and estimate the pressure broadening . At 393 K, was performed by means of Lorentz fitting. From the 18.3 GHz∙amg−1 [48] (the amg is the numerical density unit of a gas, with 1 amg = 2.6867811 × 1025 m−3, defined as the number of molecules of an ideal gas per unit volume at 273.15 K and 1 atm pressure) pressure broadening coefficient of the Rb D2 line in N2, the corresponding pressure at room temperature was found to be approximately 595 Torr, which was close to the buffer gas pressure that filled with 620 Torr during the fabrication process in Section 3.3.
The FID method was used to measure the spin precession signals and relaxation times of 129Xe and 131Xe. The test temperature was set at 393 K. First, a static magnetic field B0 = 510 nT was applied along the pump axis—that is, the Z-axis. This led to a Larmor precession of the nuclear spin along the Z-axis and the precession frequency, where is the nuclear gyromagnetic ratio of 129Xe or 131Xe. Then, a π/2 pulsed magnetic field was applied along the Y-direction to reflect the magnetization vector of the nuclear spin to the XY plane. A signal proportional to the realignment of the nuclear spin along the Z-axis was detected by the probe beam using an Rb magnetometer. Fig. 13 shows the FID and fast Fourier transform (FFT) signals for different pump-detect schemes. The measured FID signals are shown in Figs. 13(a) and (c). Eq. (8) can be used to accurately fit the obtained precession signal by constructing a sinusoidal signal-fitting function for the four frequencies and combining the natural decay exponential process of the spin-precession signal of each frequency component [20], [50].
where Mi is the amplitude of the transverse magnetization of Xe for each vibration component, the attenuation constant 1/Ti of the FID signal amplitude corresponds to the transverse relaxation time of Xe, as shown in Fig. 13, Ti is the transverse relaxation time of 129Xe or 131Xe, t0 is the initial test time, fi is the precession frequency of 129Xe or 131Xe with a nuclear quadrupole shift, φi is the initial phase of each component, and d is the fitting coefficient.
Figs. 13(b) and (d) show the frequency spectrum obtained by the FFT presented in Figs. 13(a) and (c). The FFT data were then fitted to a quadruple Lorentzian function, which is expressed as follows:
where Ai is the amplitude of each vibration component, vi is the resonance frequency of each vibration component, Γi is the resonance peak width of each vibration component, respectively, and b is used for the case where the decay center of the FID signal is not zero. The FID and FFT signal fitting codes are available in Appendix A.
Figs. 13(a) and (b) show the signals obtained using channel 1 as the pump channel and channel 2 or 3 as the probe channel, respectively. The precession frequency of some atoms shifted because of electric quadrupole splitting during the collision between 131Xe and the vapor cell wall. After fitting, the relaxation times T2 of 129Xe and 131Xe were 2.04 and 4.47 s, respectively. The Larmor precession frequencies of 129Xe and 131Xe are 6.04 and 1.73 Hz, respectively. The linewidths of the frequency peaks were 0.61 and 0.36 Hz, and the electric quadrupole frequency shift of 131Xe was 52.24 mHz. Figs. 13(c) and (d) show the signal obtained using channel 1 as the probe channel, channel 2 or 3 as the pump channel, and the frequency spectrum obtained by the FFT. The T2 values of 129Xe and 131Xe under this working scheme are 2.09 and 4.92 s, respectively. The Larmor precession frequencies were 6.05 and 1.74 Hz, and the linewidths of the frequency peaks were 0.58 and 0.36 Hz, respectively. The electric quadrupole frequency shift of 131Xe is 43.31 mHz.
4. Discussion
We did not observe the common electric quadrupole splitting peak of 131Xe in the FFT signal spectrum, because the electric quadrupole frequency shift was very small for observation; therefore, the peak of 131Xe is the combination of three resonance peaks. This phenomenon can be explained as follows:
(1) Compared with vapor cells made via the traditional glass-blowing process, which must have a residual stem after the alkali metal injection, the vapor cell used in this study was stemless. The stem breaks the symmetry of the vapor cells, particularly in small cells [19].
(2) The highly symmetric cubic MEMS atomic vapor cell cavity had a significant suppression effect on the electric quadrupole frequency shift of 131Xe [24], [51], resulting in a small frequency shift of the energy level. Thus, we observed only a single peak in the 131Xe signal.
(3) Side molding gives the vapor cell a smooth inner surface, further suppressing the electric quadrupole frequency shift compared with the rough inner surface of a conventional MEMS vapor cell cavity formed by means of dry etching [52].
The nuclear spin of 131Xe is 3/2, and the magnetic dipole interacts with the electric quadrupole. In addition, the nuclear spin of 129Xe is 1/2, and only magnetic dipole interactions are allowed. The NMR energy levels of 131Xe and 129Xe are 4 and 2, respectively. In terms of spin arrangement, fewer energy levels can ensure greater polarizability. Theoretically, compared with 129Xe, 131Xe exhibits a smaller magnetization and wider NMR width [20]. However, the NMR linewidth of the 131Xe processed in this study was much smaller than that of 129Xe. It is possible that the polarization gradient of the pump light causes the relaxation time of 129Xe to be shorter than that of 131Xe, leading to a smaller resonance linewidth for 131Xe.
The relaxation time of 131Xe observed in this experiment was small because of the wall relaxation of the 131Xe atoms. The quadrupolar interaction between the cell wall and the 131Xe nuclear spin was the main relaxation process. Filling a small amount of hydrogen (H2) or deuterium (D2) into the vapor cell to form RbH or RbD films can effectively prolong the nuclear spin-relaxation time of 131Xe, which will be the focus of our next study.
The FID signals of 129Xe and 131Xe obtained by the two pump-probe schemes show that the multi-optical-channel MEMS alkali metal vapor cell meets the requirements for multi-beam atomic devices. We used optical magnetometer signal amplitudes as a metric to compare different orientations and to estimate the optimal orientation of light coupling to the atomic sensors. The signal amplitudes of the two pump-probe schemes after the tests are shown in Figs. 13(a) and (c). When channel 1 was the pump channel and channel 2 or 3 was the probe channel, the signal amplitude of the channel 1 pumping was approximately 1.3 times that of the channel 2 or 3 pumping. This difference in the signal amplitude was attributed to a decrease in the lateral transmittance, resulting in a decrease in the pumping efficiency. In the future, we will further optimize the molding process to improve the lateral transmittance.
5. Conclusions
In this study, the batch manufacturing of MEMS alkali-metal atomic vapor cells with multiple optical channels was accomplished using multilayer anodic bonding and precise side-molding technology. A highly transparent vapor-cell sidewall was created using a process that included laser drilling, wet etching, and precision side molding. The required laser wavelength transmittance was approximately 89.5%, the average surface roughness was 22 nm, and the average surface shape accuracy PV was 283 nm. Alkali-metal injection, noble-gas recycling, and recovery systems were developed, and a test platform with a double-path laser was constructed. The absorption spectrum demonstrated the precise control of gas pressure. The developed multi-optical-channel alkali metal vapor cell satisfies the requirements of multi-optical-channel atomic devices, based on the results of the measurement and a comparison of the FID signals of 129Xe and 131Xe of the processed vapor cell under various pump-probe schemes.
Acknowledgments
This work was supported in part by the National Key Research and Development Plan (2022YFB3203400), the National Natural Science Foundation of China (62103324 and U1909221), and the Natural Science Foundation of Shaanxi (2022JQ-554).
Compliance with ethics guidelines
Mingzhi Yu, Yao Chen, Yongliang Wang, Xiangguang Han, Guoxi Luo, Libo Zhao, Yanbin Wang, Yintao Ma, Shun Lu, Ping Yang, Qijing Lin, Kaifei Wang, and Zhuangde Jiang declare that they have no conflict of interest or financial conflicts to disclose.
M.W.Doherty, N.B.Manson, P.Delaney, F.Jelezko, J.Wrachtrup, L.C.L.Hollenberg. The nitrogen-vacancy colour centre in diamond. Phys Rep, 528 (1) (2013), pp. 1-45.
[4]
R.C.Jaklevic, J.Lambe, A.H.Silver, J.E.Mercereau. Quantum interference effects in Josephson tunneling. Phys Rev Lett, 12 (7) (1964), pp. 159-160.
[5]
R.C.Jaklevic, J.Lambe, J.E.Mercereau, A.H.Silver. Macroscopic quantum interference in superconductors. Phys Rev, 140 (5A) (1965), p. A1628.
[6]
W.E.Bell, A.L.Bloom. Optical detection of magnetic resonance in alkali metal vapor. Phys Rev, 107 (6) (1957), pp. 1559-1565.
[7]
E.B.Alexandrov, V.A.Bonch-Bruevich. Optically pumped atomic magnetometers after three decades. Opt Eng, 31 (4) (1992), pp. 711-717.
[8]
K.W.Martin, G.Phelps, N.D.Lemke, M.S.Bigelow, B.Stuhl, M.Wojcik, et al. Compact optical atomic clock based on a two-photon transition in rubidium. Phys Rev Appl, 9 (1) (2018), 014019.
[9]
Z.L.Newman, V.Maurice, T.Drake, J.R.Stone, T.C.Briles, D.T.Spencer, et al. Architecture for the photonic integration of an optical atomic clock. Optica, 6 (5) (2019), pp. 680-685.
[10]
A.Godone, F.Levi, C.E.Calosso, S.Micalizio. High-performing vapor-cell frequency standards. Riv Nuovo Cimento, 38 (3) (2015), pp. 133-171.
[11]
I.K.Kominis, T.W.Kornack, J.C.Allred, M.V.Romalis. A subfemtotesla multichannel atomic magnetometer. Nature, 422 (6932) (2003), pp. 596-599.
V.Shah, S.Knappe, P.D.D.Schwindt, J.Kitching. Subpicotesla atomic magnetometry with a microfabricated vapour cell. Nat Photonics, 1 (2007), pp. 649-652.
[14]
T.W.Kornack, R.K.Ghosh, M.V.Romalis. Nuclear spin gyroscope based on an atomic comagnetometer. Phys Rev Lett, 95 (23) (2005), 230801.
[15]
S.S.Sorensen, D.A.Thrasher, T.G.Walker. A synchronous spin-exchange optically pumped NMR-gyroscope. Appl Sci, 10 (20) (2020), p. 7099.
[16]
W.E.Bell, A.L.Bloom, J.Lynch. Alkali metal vapor spectral lamps. Rev Sci Instrum, 32 (6) (1961), pp. 688-692.
[17]
M.V.Balabas, D.Budker, J.Kitching, P.D.D.Schwindt, J.E.Stalnaker. Magnetometry with millimeter-scale antirelaxation-coated alkali-metal vapor cells. J Opt Soc Am B, 23 (6) (2006), pp. 1001-1006.
[18]
N.Castagna, G.Bison, G. diDomenico, A.Hofer, P.Knowles, C.Macchione, et al. A large sample study of spin relaxation and magnetometric sensitivity of paraffin-coated Cs vapor cells. Appl Phys B, 96 (4) (2009), pp. 763-772.
[19]
Z.Wu. Wall interactions of spin-polarized atoms. Rev Mod Phys, 93 (3) (2021), 035006.
[20]
Y.K.Feng, S.B.Zhang, Z.T.Lu, D.Sheng. Electric quadrupole shifts of the precession frequencies of 131Xe atoms in rectangular cells. Phys Rev A, 102 (4) (2020), 043109.
[21]
T.G.Walker, M.S.Larsen. Spin-exchange-pumped NMR gyros. E. Arimondo, C.C. Lin, S.F. Yelin (Eds.), Advances in atomic, molecular, and optical physics, 65, Elsevier, Cambridge (2016), pp. 373-401.
Z.Wu, W.Happer, M.Kitano, J.Daniels. Experimental studies of wall interactions of adsorbed spin-polarized 131Xe nuclei. Phys Rev A, 42 (5) (1990), pp. 2774-2784.
[24]
E.A.Donley, J.L.Long, T.C.Liebisch, E.R.Hodby, T.A.Fisher, J.Kitching. Nuclear quadrupole resonances in compact vapor cells: the crossover between the NMR and the nuclear quadrupole resonance interaction regimes. Phys Rev A, 79 (1) (2009), 013420.
[25]
J.Kitching, S.Knappe, L.Hollberg. Miniature vapor-cell atomic-frequency references. Appl Phys Lett, 81 (3) (2002), pp. 553-555.
[26]
L.A.Liew, S.Knappe, J. Moreland, H.Robinson, L.Hollberg, J.Kitching. Microfabricated alkali atom vapor cells. Appl Phys Lett, 84 (14) (2004), pp. 2694-2696.
[27]
S.Knappe, V.Gerginov, P.D.D.Schwindt, V.Shah, H.G.Robinson, L.Hollberg, et al. Atomic vapor cells for chip-scale atomic clocks with improved long-term frequency stability. Opt Lett, 30 (18) (2005), pp. 2351-2353.
[28]
L.A.Liew, J. Moreland, V.Gerginov. Wafer-level filling of microfabricated atomic vapor cells based on thin-film deposition and photolysis of cesium azide. Appl Phys Lett, 90 (11) (2007), 114106.
[29]
M.Hasegawa, R.K.Chutani, C.Gorecki, R.Boudot, P.Dziuban, V.Giordano, et al. Microfabrication of cesium vapor cells with buffer gas for MEMS atomic clocks. Sens Actuators A, 167 (2) (2011), pp. 594-601.
[30]
V.Maurice, J.Rutkowski, E.Kroemer, S.Bargiel, N.Passilly, R.Boudot, et al. Microfabricated vapor cells filled with a cesium dispensing paste for miniature atomic clocks. Appl Phys Lett, 110 (16) (2017), 164103.
[31]
R.Vicarini, V.Maurice, M. AbdelHafiz, J.Rutkowski, C.Gorecki, N.Passilly, et al. Demonstration of the mass-producible feature of a Cs vapor microcell technology for miniature atomic clocks. Sens Actuators A, 280 (2018), pp. 99-106.
[32]
M.Bulatowicz, R.Griffith, M.Larsen, J.Mirijanian, C.B.Fu, E.Smith, et al. Laboratory search for a long-range T-odd, P-odd interaction from axionlike particles using dual-species nuclear magnetic resonance with polarized 129Xe and 131Xe gas. Phys Rev Lett, 111 (10) (2013), 102001.
[33]
M.E.Limes, D.Sheng, M.V.Romalis. 3He-129Xe comagnetometery using 87Rb detection and decoupling. Phys Rev Lett, 120 (3) (2018), 033401.
[34]
E.J.Eklund, A.M.Shkel. Glass blowing on a wafer level. J Microelectromech Syst, 16 (2) (2007), pp. 232-239.
[35]
E.J.Eklund, A.M.Shkel, S.Knappe, E.Donley, J.Kitching. Glass-blown spherical microcells for chip-scale atomic devices. Sens Actuators A, 143 (1) (2008), pp. 175-180.
[36]
M.A.Perez, U.Nguyen, S.Knappe, E.A.Donley, J.Kitching, A.M.Shkel. Rubidium vapor cell with integrated Bragg reflectors for compact atomic MEMS. Sens Actuators A, 154 (2) (2009), pp. 295-303.
[37]
M.A.Perez, J.Kitching, A.M.Shkel. Design and demonstration of PECVD multilayer dielectric mirrors optimized for micromachined cavity angled sidewalls. Sens Actuators A, 155 (1) (2009), pp. 23-32.
[38]
R.Han, Z.You, F.Zhang, H.Xue, Y.Ruan. Microfabricated vapor cells with reflective sidewalls for chip scale atomic sensors. Micromachines, 9 (4) (2018), p. E175.
[39]
W.Ming, Z.Chen, J.Du, Z.Zhang, G.Zhang, W.He, et al. A comprehensive review of theory and technology of glass molding process. Int J Adv Manuf Technol, 107 (5,6) (2020), pp. 2671-2706.
[40]
Y.Zhang, G.Yan, K.You, F.Fang. Study on α-Al2O3 anti-adhesion coating for molds in precision glass molding. Surf Coat Technol, 391 (2020), 125720.
[41]
schott.com [Internet]. Shanghai: SCHOTT AG; c2023[cited 2020 Aug 20]. Available from: https://www.schott.com/zh-cn/products/borofloat
[42]
keyence.com.cn [Internet]. Shanghai: KEYENCE CORPORATION; c2023 [cited 2020 Aug 20]. Available from: https://www.keyence.com.cn/ss/products/microscope/roughness/surface/tab01_e.jsp
[43]
A.T.Dellis, V.Shah, E.A.Donley, S.Knappe, J.Kitching. Low helium permeation cells for atomic microsystems technology. Opt Lett, 41 (12) (2016), pp. 2775-2778.
[44]
S.Karlen, J.Gobet, T.Overstolz, J.Haesler, S.Lecomte. Lifetime assessment of RbN3-filled MEMS atomic vapor cells with Al2O3 coating. Opt Express, 25 (3) (2017), pp. 2187-2194.
[45]
N.Dural, M.V.Romalis. Gallium phosphide as a new material for anodically bonded atomic sensors. APL Mater, 2 (8) (2014), 086101.
[46]
M.D.Rotondaro, G.P.Perram. Collisional broadening and shift of the rubidium D1 and D2 lines (52S12→52P12, 52P32) by rare gases, H2, D2, N2, CH4 and CF4. J Quant Spectrosc Radiat Transfer, 57 (4) (1997), pp. 497-507.
[47]
R.M.Noor, M.H.Asadian, A.M.Shkel. Design considerations for micro-glassblown atomic vapor cells. J Microelectromech Syst, 29 (1) (2020), pp. 25-35.
[48]
M.V.Romalis, E.Miron, G.D.Cates. Pressure broadening of Rb D1 and D2 lines by 3He, 4He, N2, and Xe: line cores and near wings. Phys Rev A, 56 (6) (1997), pp. 4569-4578.
[49]
H.Zheng, W.Quan, X.Liu, Y.Chen, J.Lu. Measurement of atomic number of alkali vapor and pressure of buffer gas based on atomic absorption. Spectrosc Spectral Anal, 35 (2) (2015), pp. 507-511.
[50]
T.M.Kwon, J.G.Mark, C.H.Volk. Quadrupole nuclear spin relaxation of 131Xe in the presence of rubidium vapor. Phys Rev A, 24 (4) (1981), pp. 1894-1903.
[51]
P.A.Heimann, I.A.Greenwood, J.H.Simpson. Quadrupole perturbation effects upon the 201Hg magnetic resonance. Ⅱ.relaxation due to an anisotropic perturbation. Phys Rev A, 23 (3) (1981), pp. 1209-1214.
[52]
R.Butscher, G.Wackerle, M.Mehring. Nuclear quadrupole interaction of highly polarized gas-phase 131Xe with a glass surface. J Chem Phys, 100 (9) (1994), pp. 6923-6933.