We present a highly stable bow-tie power enhancement cavity for critical second harmonic generation (SHG) into the UV using a Brewster-cut β-BaB2O4 (BBO) nonlinear crystal. The cavity geometry is suitable for all UV wavelengths reachable with BBO and can be modified to accommodate anti-reflection coated crystals, extending its applicability to the entire wavelength range accessible with non-linear frequency conversion. The cavity is length-stabilized using a fast general purpose digital PI controller based on the open source STEMlab 125-14 (formerly Red Pitaya) system acting on a mirror mounted on a fast piezo actuator. We observe 130 h uninterrupted operation without decay in output power at 313 nm. The robustness of the system has been confirmed by exposing it to accelerations of up to 1 g with less than 10% in-lock output power variations. Furthermore, the cavity can withstand 30 min of acceleration exposure at a level of 3 grms without substantial change in the SHG output power, demonstrating that the design is suitable for transportable setups.
I. INTRODUCTION
Today’s quantum optics experiments require a broad range of laser frequencies. Usually, several continuous wave laser systems are required per experiment1–3 which imposes tight bounds on the reliability and cost per system. Recent progress in the development of quantum sensors has led to systems outperforming classical devices, such as quantum gravimeters,4,5 and clocks exceeding 10−17 accuracy.2,6,7 New applications such as relativistic geodesy8 using transportable optical clocks9–11 or taking advantage of the long interaction times in atom interferometry experiments in a microgravity or even space environment12–14 have emerged. For these applications, quantum optics experiments need to be operated outside highly specialized laboratories, increasing the demands in terms of mechanical robustness of the optical setups. Small and at the same time reliable laser sources are available only for a restricted wavelength range.15,16 To reach other wavelengths, non-linear conversion processes, such as sum or difference frequency or harmonic generation,17–25 are typically employed. A common approach to generate the desired wavelength is second harmonic generation26 (SHG) of External Cavity Diode Lasers (ECDLs) or fiber lasers. The frequency doubling of infrared (IR) lasers has been demonstrated in single-crystal monolithic ring cavities.27–29 With the advent of commercially available periodically poled waveguide doublers, e.g., based on lithium tantalate or lithium niobate nonlinear crystals, wavelength conversion into the blue spectral range has been achieved in compact setups.30 Monolithic ring cavities have been demonstrated to produce blue light down to wavelengths of 429 nm.31–33 However, to our knowledge, neither single-crystal monolithic ring cavities nor modules are available for UV generation below 350 nm. For these UV wavelengths, SHG in nonlinear crystals such as Beta-Barium-Borate (β-BaB2O4, BBO), placed in an optical enhancement resonator, are typically employed.34 Commercially available systems usually come with restricted flexibility, e.g., lacking access to the SHG light internally reflected at the output facet of Brewster-cut nonlinear crystals. By contrast, self-built systems made from off-the-shelve components usually lack mechanical stability and reliability.
Traditionally, the length of the cavity is kept interferometrically stable by displacing one of the cavity mirrors in a proportional-integral (PI) feedback loop usually implemented using analog electronics. With the availability of fast general-purpose digital hardware such as Field Programmable Gate Arrays (FPGAs) and/or microcontroller-based systems, it has become popular to employ digital feedback controllers instead.35–39 These devices are more universally applicable, feature higher usability, and can easily be linked to the existing experimental control infrastructure.
Here we report on a mechanical monolithic bow-tie cavity design40 for critically phase-matched SHG in a BBO crystal and its general purpose digital PI locking electronics based on a modified STEMlab 125-14 application.41,42 The cavity is implemented in a robust monolithic mechanical support frame and equipped with a minimal set of high-quality adjustment screws that are accessible from outside. The cavity geometry is suitable for the generation of all wavelengths accessible with BBO in SHG ooe-processes, where the pump beam is polarized normal (ordinary, o) to the principal plane of the crystal and the SHG beam in the plane (extraordinary, e), starting at 204.8 nm.43 The internal reflections of the pump light (PL) at the fundamental frequency and SHG on the crystal’s facets are accessible through two additional windows. In order to prevent moisture-induced fogging of the crystal, a dry purging gas, such as nitrogen or oxygen, can be applied inside the sealed cavity.
We measure a locking bandwidth of 17 kHz and demonstrate 130 h continuous operation in lock for the conversion from 626 nm to 313 nm without substantial decay in the output power. Moreover, the cavity remains in lock while being exposed to accelerations in the vertical direction of around 1 g. A comparison of the output power before and after a 30 min high-acceleration shaking in one dimension demonstrates that the cavity can withstand a truck transport according to ISO 13355:2016 without significant changes in the alignment.
The paper is structured as follows: In Sec. II, important aspects of the theory of SHG and bow-tie power enhancement cavities are presented. In Sec. III, we describe the monolithic cavity design, and in Sec. IV, the setup of the complete system including the STEMlab-based locking electronics is described. The results on acceleration testing and long-term stability are presented in Sec. V.
II. THEORETICAL BACKGROUND
When light passes through a nonlinear medium, the SHG light of twice the frequency can be generated. The nonlinear conversion process scales with the square of the PL intensity. In the case of low conversion efficiency and therefore undepleted PL, the SHG power PSHG scales with the square of the pump light power PPL,
where κ is the conversion coefficient. Efficient SHG is only achieved if the wavevectors of the pump light and the SHG fulfill the phase-matching condition
In BBO, this can be achieved by taking advantage of the crystal’s birefringence with different refractive indices no and ne for different incident polarisations. By choosing the appropriate phase-matching angle θpm within the resulting index ellipsoid, nSHG = nPL can be achieved.44 For a SHG of 626 nm in BBO, the resulting phase-matching angle is
A. SHG power optimization
The single-pass conversion efficiency of cw light in a BBO crystal of a few mm length in general is rather low (on the order of κ ∼ 10−4 1/W). Therefore, an enhancement cavity for the PL is built around the nonlinear crystal.26,45,46 The optimization of the cavity SHG output power can be divided into two steps: first, the optimization of the power generated per single pass of the PL through the crystal47 and second, the optimization of the power enhancement.48
1. Single pass optimization
As mentioned above, the SHG power generated per infinitesimal crystal volume is proportional to the square of the PL intensity. Tighter focusing of the PL enhances the intensity at the focus, at the expense of lower intensity away from the focus. This trade-off for a given crystal length l has been investigated for the propagation of a circular Gaussian beam with minimum waist by Boyd and Kleinman,47 who derived an optimum focusing ratio of l/b = 2.84 in the absence of birefringence, where l is the crystal length and is the confocal parameter with k = 2π/λ. However, in the presence of birefringence, the phase vector and Poynting vector of the extraordinary SHG wave are in general not parallel, which results in a walk-off angle ϱ between the SHG and PL beam of
for a SHG of 626 nm in BBO. This effect is quantified by the walk-off parameter B,47
A strong walk-off (such as in BBO) requires weaker focusing of the pump light to optimize the spatial overlap between PL and SHG. Furthermore, a focused Gaussian beam exhibits a phase deviation Δk from an ideal plane wave that changes along the propagation direction, known as the Gouy effect.49 It is typically quantified by the parameter . For the SHG considered here, this means that the degree to which the phase-matching condition is fulfilled changes over the crystal length. Taking all these effects into account through integration over infinitesimal contributions of the PL to the SHG light when propagating through the crystal, the Boyd-Kleinman theory provides optimum parameters l/b and σ for single-pass SHG output power. For the SHG from 626 nm to 313 nm in a 10 mm long BBO crystal, we numerically obtain l/b = 1.42, σ = 0.75, and ω0 = 19 μm with a single-pass conversion efficiency of κ = 1.1·10−4. A similar analysis can be performed in the case of an elliptical focus inside the crystal,50 as is the case for astigmatic cavities and/or Brewster-cut crystals.
2. Cavity geometry
A nonlinear medium with low conversion efficiency placed in an optical cavity in which the PL power is enhanced by a factor on the order of 50 and can provide up to 2500 times more SHG output power than is possible in single pass configuration.
A common design for an optical enhancement cavity, referred to as “bow-tie,” consists of four mirrors and is schematically shown in Fig. 1. Compared to a linear cavity made of two mirrors, this geometry avoids the formation of a standing PL wave and thus reduces the likelyhood of photo-refractive effects. Furthermore, it allows for astigmatism compensation of Brewster-cut crystals by choosing a suitable distance between the two mirror pairs involved (see below). The incoupling mirror (1) fulfills the impedance matching condition (see below). Mirror (2) is mounted on a piezoelectric actuator with a maximum stroke sufficient to change the cavity length by more than one free spectral range (FSR) to enable length stabilization of the cavity to the PL wavelength. It is lightweight to achieve high mechanical resonance frequencies and thus a high feedback bandwidth. Mirrors (3) and (4) are concave, and their distance is chosen to provide a focus of the appropriate size at the center of the crystal. The fourth mirror also acts as an outcoupler for the SHG light. A fully monolithic bow-tie cavity design similar to the Nd:YAG NPRO design27 would suffer from a detrimental PL loss due to the high linear absorption of αPL ∼ 0.01 cm−1 for λ = 626 nm in BBO.51
Schematic layout of a bow-tie cavity. The resonator is pumped through mirror (1) with light focused on the incoupling waist ωic. A lightweight mirror (2) on a piezoelectric actuator is used to lock the length of the cavity. (3) and (4) focus the light to the waist ω0 inside the nonlinear crystal and mode match the beam for the next round trip.
Schematic layout of a bow-tie cavity. The resonator is pumped through mirror (1) with light focused on the incoupling waist ωic. A lightweight mirror (2) on a piezoelectric actuator is used to lock the length of the cavity. (3) and (4) focus the light to the waist ω0 inside the nonlinear crystal and mode match the beam for the next round trip.
3. Crystal shape
For maximum PL power enhancement and therefore a high SHG output power [cf. Eq. (1)], low PL loss inside the cavity from absorption or scattering is required. One major source of loss in circulating PL is reflection at the interfaces of the nonlinear medium and the surrounding gas. These reflections can be avoided either by the application of an antireflective coating on the surfaces or by cutting and mounting the crystal under Brewster’s angle to the incident beam. While anti-reflection (AR) coatings protect hygroscopic crystals such as BBO against moisture in the environment, they are prone to damage induced by UV or pump light. Furthermore, residual reflections on surfaces orthogonal to the PL beam can cause a backward traveling PL wave. Weighting these properties and given the option to operate the crystal in a dry environment, it seems advantageous to choose a crystal cut under Brewster’s angle θB for the PL. For a SHG of 626 nm in BBO, the resulting angle is
where nPL is the index of refraction for the PL. Since PL and SHG have orthogonal polarization in ooe SHG, the Brewster condition can only be fulfilled for the pump light. The corresponding fractional internal reflection of the SHG at the crystal-gas interface is given by Fresnel’s equations for s-polarized light that simplify to
4. Impedance matching
Assuming ideal mirrors and input beam mode matching, the entire pump light can be coupled into the cavity if the incoupling mirror transmission is chosen in such a way that the coupled pump light compensates all losses inside the cavity during a round trip. In the case of a SHG cavity, the latter consist of two parts: PL loss due to conversion and parasitic loss due to imperfect reflection inside the resonator and linear absorption in the nonlinear crystal and surrounding gas with absorption coefficient αPL.
The total round trip transmission of the cavity as a function of the single cavity mirror reflectivity RM and the crystal facet transmission TC without considering conversion into SHG is
Taking into account the loss through conversion with a factor , the optimal incoupling mirror reflectivity Rin for impedance matching is52
where the small correction factor t from the incoupling mirror has been applied.
5. Cavity parameters
Both the non-normal incidence of the cavity mode on the curved mirrors and crystal facets cut under Brewster’s angle introduce astigmatism in the cavity which needs to be accounted for in the design. The system forms two effective cavities in the xy (sagittal) and xz (tangential) planes which both need to be stable. Optimizing for simultaneous stability and conversion efficiency23,50 for a 10 mm BBO crystal and a 50 mm radius of curvature of the mirrors yields a mechanical cavity footprint which is inconveniently large to realize experimentally in a compact way. Instead, we choose the cavity parameters outlined in Table I. and are the saggital and tangential waists, respectively, which are larger than the optimum waist ω0 from Boyd-Kleinman. This has the added advantage of reducing crystal damage from high intensities at the crystal surfaces. In the end, the actual cavity geometry is defined by l, θB, the mirror to crystal distance dmc, θcm, and the geometric round-trip length lgeo.
Geometric parameters of the doubling cavity using a 10 mm BBO-crystal pumped with PPL = 0.5 W at λPL = 626 nm and r = 50 mm mirrors.
Phase-matching angle | θpm = 38.4° |
Waist saggital | ω0,xy = 25.3 μm |
Waist tangential | ω0,xz = 39.6 μm |
Brewster angle | θB = 59.0° |
Angle on curved mirrors | θcm = 15.7° |
Mirror to crystal distance | dmc = 25.1 mm |
Geometric round trip length | lgeo = 304.8 mm |
Input coupling waist | ωic = 166 μm |
Phase-matching angle | θpm = 38.4° |
Waist saggital | ω0,xy = 25.3 μm |
Waist tangential | ω0,xz = 39.6 μm |
Brewster angle | θB = 59.0° |
Angle on curved mirrors | θcm = 15.7° |
Mirror to crystal distance | dmc = 25.1 mm |
Geometric round trip length | lgeo = 304.8 mm |
Input coupling waist | ωic = 166 μm |
For this choice of parameters, both the sagittal and tangential cavities are stable simultaneously and the input coupling waist (between the two flat mirrors) is nearly circular with a waist of 166 μm. Compared to the geometry for optimal efficiency,23,50 the expected loss in conversion efficiency is only about 12%. This moderate loss is outweighed by the gain in mechanical stability and ruggedness due to the smaller footprint. The waists and Brewster’s angle have been calculated for a SHG of λ = 626 nm, but their dependence on the PL wavelength is rather weak. For λ = 534 nm PL, for instance, the waists change by −1.9 μm and −2.9 μm. Brewster’s angle increases by 0.1°. Therefore, the cavity geometry is practically suitable for frequency doubling of a spectral range of several hundred nm, as long as a 10 mm long BBO-crystal cut under the corresponding phase-matching angle is used. In conclusion, the shape of the monolithic housing (see below) is universal for SHG to the UV.
B. Length stabilization and its limit
During operation, the length of the cavity is stabilized to the PL wavelength via a voltage applied to the piezoelectric actuator using a PI controller. The required error signal is generated using the Hänsch-Coulliaud (HC) locking scheme53 that generates a dispersion-shaped error signal from the phase shift upon reflection of light from the cavity. P-polarized light is coupled into the cavity and experiences a varying phase shift near the cavity resonance, while s-polarized light does not couple into the cavity and thus obtains a constant π phase shift upon reflection. Compared to the Pound-Drever-Hall54 locking scheme, no sidebands are imprinted on the PL and thus are also absent in the SHG, which is an advantage for certain applications.
The locking bandwidth is fundamentally limited by the mechanical resonance frequency f* of the piezoelectric actuator and the attached mirror. It is given as55
where f0 is the mechanical resonance frequency of the piezoelectric actuator, m is its mass, and M is the mass of the attached mirror.
For the monolithic cavity design presented here, we expect only relatively small changes of the round-trip length. Therefore a short piezoelectric actuator with low weight and high resonance frequency seems to be sufficient to keep the cavity length stabilized under moderate environmental perturbations.
III. CAVITY DESIGN
The main design goals of the cavity are “turn-key” operation after being exposed to mechanical environmental conditions typical for long-distance truck and plane transportation. Furthermore, robust locking and operation in a non-lab environment is important to be a reliable part of a transportable quantum optics experiment.
Therefore, the main body is milled from a single Al block shown in Fig. 2. Two mirror holder front plates mounted via 3 μm screws (150 μm displacement per turn) and two springs each, all mounted directly on the main body, feature the four degrees of freedom (DOFs) required to close the beam path. In combination with the two other mirrors that are mounted on fine-threaded aluminum cylinders, the cavity round trip length can be adjusted on the few mm-scale while keeping the angles between the beams and therefore the astigmatism compensation constant. The mirrors are mounted using stress-free retaining rings to avoid birefringence caused by mechanical stress. This does not apply for the mirror behind the incoupler which is glued (Thorlabs 353NDPK Epoxy) on a piezoelectric actuator (Thorlabs AE0505D08F) to enable cavity length stabilization. In crystals used for SHG into the UV, crystal inhomogeneities and degradation effects have been observed. Compensation of these effects requires two translational DOFs parallel to the front surface of the crystal. With a third DOF, the waist position in the crystal along the direction of the PL propagation is adjusted. Two additional rotational DOFs allow the fine adjustment of the phase-matching and Brewster’s angle. All these DOFs are provided by a commercially available 5-axis-aligner (NEWPORT 9081-M) that is mounted directly in the main body. The crystal is mounted in two shells that are mounted to a lever connected to the aligner. These shells and the lever are equipped with channels that lead oxygen from a supply connected to the main body directly to the two facets of the crystal to prevent light-induced damage.56 Additionally, the tip of the lever can be actively temperature stabilized using a thick-film resistor (Vishay Sfernice RTO 20) and a sensor (TDK B57861S) connected via a standard D-sub feedthrough to a suitable PI controller to prevent condensate formation on the crystal facets that would likely cause degradation of hygroscopic crystals such as BBO. Finally, the cavity is equipped with FKM sealings to prevent unwanted substances from entering the enhancement cavity and cause degradation of the crystal and/or mirrors. Therefore, a gas outlet in the main body ensures a moderate gas exchange rate and overpressure inside the housing. The oxygen concentration can be measured by a sensor (Greisinger GOX 100) located in the lid of the main body.
Schematic cross section of the SHG cavity. All mechanical components are mounted directly to the monolithic housing. The incoupling and outcoupling mirrors are tiltable using micrometer screws while the other mirrors can be translated only. The temperature-stabilized BBO crystal is mounted on a 5-axis aligner. Ar-coated windows allow for optical access to the sealed cavity.
Schematic cross section of the SHG cavity. All mechanical components are mounted directly to the monolithic housing. The incoupling and outcoupling mirrors are tiltable using micrometer screws while the other mirrors can be translated only. The temperature-stabilized BBO crystal is mounted on a 5-axis aligner. Ar-coated windows allow for optical access to the sealed cavity.
IV. SETUP
Figure 3 shows a schematic overview of the complete setup for SHG at 313 nm. It consists of a sum frequency generation (SFG) setup for the 626 nm pump light, the SHG cavity setup, the locking electronics, and a computer for power recording. The general approach is similar to that of Wilson et al.23
Schematic overview of the complete setup for 313 nm generation. FCi: fiber coupler, Li: lens, Mi: mirror, PBS: polarizing beamsplitter cube, /: half/quarter wave plate, ND: neutral density filter. Left: pump light generation at 626 nm by SFG of two IR fiber lasers. Center: 313 nm generation in the resonant enhancement bow-tie cavity locked by the Hänsch-Coulliaud locking scheme (dashed beam path and digital multi purpose PI controller STEMlab-Box). For output power stability monitoring, a fraction of the pump and SHG light is picked off and measured via photodiodes.
Schematic overview of the complete setup for 313 nm generation. FCi: fiber coupler, Li: lens, Mi: mirror, PBS: polarizing beamsplitter cube, /: half/quarter wave plate, ND: neutral density filter. Left: pump light generation at 626 nm by SFG of two IR fiber lasers. Center: 313 nm generation in the resonant enhancement bow-tie cavity locked by the Hänsch-Coulliaud locking scheme (dashed beam path and digital multi purpose PI controller STEMlab-Box). For output power stability monitoring, a fraction of the pump and SHG light is picked off and measured via photodiodes.
A. Pump light setup
Up to 1.1 W, pump light (PL) at 626 nm is generated by sum frequency generation (SFG) of two fiber lasers at 1550 nm and 1051 nm in a periodically-poled lithium niobate (PPLN) crystal described in Ref. 57 and coupled in a polarization maintaining fiber at FC1. This fiber leads to a 12.7 mm thick (450 mm)2 aluminum breadboard with all the optical components for the 313 nm SHG. While this breadboard size has been chosen for convenience, it seems possible to reduce the overall footprint of the setup by a factor of three by optimizing the beam path layout. The breadboard as well as the aforementioned SFG is covered by two boxes made of aluminum. Both are mounted on a standard air damped optical table in a laboratory temperature stabilized at the 1 K level. The breadboard holding the cavity is placed on 5 mm thick viscoelastic damping rubber (Sorbothane®).
B. SHG setup
Behind the PL fiber coupler FC2, a telescope consisting of two lenses L1 and L2 with focal lengths f = +50 mm and f = +30 mm, respectively, mounted on single axis translation stages is used for the mode matching with the PL cavity mode. A 1% pickup beam is used to monitor the PL power. The two highly stable mirrors M1 and M2 provide the four DOFs required for coupling the PL into the cavity. Both mirrors are mounted using retaining rings to avoid birefringence due to mechanical stress. A polarizing beamsplitter (PBS) cube cleans the polarization of the PL and a subsequent λ/2 wave plate rotates it into the plane required for the SHG. Behind the outcoupling mirror of the SHG cavity, a cylindrical lens CL with focal length f = +100 mm is employed for correcting the astigmatism arising from the walk-off of the SHG light generated inside the BBO crystal. A second mirror (M3) of the same type as the outcoupler is used to filter the transmitted PL from the SHG beam. For the test presented here, a thermal powermeter measurement head is used to collect the SHG light. The PL reflected by the cavity is attenuated by a 3.0 absorptive neutral density (ND) filter, focussed by lens L3 with focal length f = +250 mm, and sent through a polarization analyzer consisting of a half-wave plate, a quarter-wave plate, and a Wollaston prism to generate two orthorgonally polarized beams of similar intensity carrying the information about the phase relation of the PL entering the cavity.53 These two beams hit a differential photodiode. The resulting signal is the error signal for the digital PI controller. The small fraction of the PL reflected on the front surface of the crystal due to imperfect polarization and scatter is focused on a photodiode to generate a signal as a measure for the circulating power.
C. Digital PI controller design
The error signal generated by the differential photodiode is fed into a digital PI controller based on a modified STEMlab 125-14 (formerly: Red Pitaya). A real-time digital PID controller algorithm is implemented in the FPGA as part of the hardware configuration. It operates independently from the embedded software on the STEMlab. The parameters of the digital controller are set either through a remote network connection or using a standalone user interface including a touchscreen and rotary encoder inputs. For the purpose of cavity locking, the D-part of the controller is set to zero. Besides the PID controller, a function generator was implemented. It can output sine, square, or triangular waveforms within frequency range of 0–50 MHz.
1. Hardware
STEMlab58 hardware is not open-source, but its processing core and peripherals give many options for low-level customization. The hardware platform is built around the Xilinx Zynq7010 All Programmable System-on-Chip (SoC),59,60 which combines a dual-core ARM Cortex A9 processor and a field programmable gate array (FPGA) in one chip. Its peripherals include 16 general-purpose input/output (GPIO) lines, various digital interfaces such as inter-integrated circuit (I2C), serial peripheral interface bus (SPI), and universal asynchronous receiver-transmitter (UART), and four 12-bit analog-to-digital converters (ADCs) and DACs operating at 100 kS/s. All of these signal lines are accessible on the extension connectors of the board. In addition, there are two high-speed communication interfaces based on Ethernet and universal serial bus (USB). The STEMlab provides two fast analog input and output channels, operating at a sampling rate of 125 MS/s and having 14-bit nominal resolution. They are implemented using fast ADCs (LTC2145)61 and DACs (AD9767),62 plus various analog signal conditioning circuits. Figure 4 shows an overview of the components in the STEMlab system that are most relevant for the implemented digital controller. The blue boxes indicate the two principal domains of the Zynq, where custom hardware and software is implemented.
The main characteristics of the fast analog frontend channels are given in Table II. They are the most critical interfaces because their conversion errors determine the achievable performance of the digital controller. Apart from the effects of quantization that would be present even in ideal converters, various other error sources in the ADCs and DACs have to be considered. Among them are static gain and offset errors, as well as temperature-dependent drifts.
To obtain a stand-alone controller unit based on the STEMlab, we added a user interface consisting of an LCD display with touchscreen and a rotary encoder. The voltage range of the DAC outputs was increased using non-inverting amplifiers with optional external low-pass filters for control applications with reduced bandwidth. Using amplifiers with gain G = 8, the output range of 2 Vpp was scaled to 16 Vpp, sufficient to cover two free spectral ranges of the monolithic doubling cavity, when applied to the cavity piezoelectric actuator. To reduce the output noise of the STEMLab, we implemented the modification described in Ref. 64, which eliminates the excess noise coupling from the digital supply rail to the analog output. This improvement led to a noise reduction from 3 to 4 LSB peak-to-peak to 1-2 LSB peak-to-peak. A certain variation of excess noise was observed between different units. In addition to the described circuits, the standalone box contained separate power supplies for the STEMlab board, its cooling fan, and the output amplifiers. The piezoelectric actuator inside the cavity was connected with a series resistor of 10 Ω, which, together with the capacitance of the actuator, forms a low-pass filter with a cutoff frequency of 21.2 kHz, preventing noise-induced oscillations at high frequency.
2. Software and FPGA configuration
The STEMlab embedded system uses the Linux operating system (presently the Ubuntu 16.04 distribution). The operating system (OS) image is loaded from a Secure Digital (SD) card which also contains the customized software and the FPGA configuration.
The open-source web applications provided with the STEMlab, as well as the PID controller application, can be accessed from any web browser via a network connection. The digitized input signals can be monitored in real time. All parameters of the digital PID controller can be set remotely through the web application, via serial connection, or using the local user interface when operating in standalone mode.
To achieve low-latency digital control, the calculations for the PID algorithm are implemented in the FPGA hardware. Hard real-time processing is possible due to the deterministic timing behavior of the FPGA. The original STEMlab PID application contains a Multiple Input Multiple Output (MIMO) PID controller which consists of four standard PID controllers with P, I, and D parameter settings and integrator reset control.65 The PID parameters can be set in the user interface and are then sent to the FPGA registers described in the register map.66 The implementation of the four discrete PID controllers can be approximately described by the equations for a continuous-time control system,
where m(t) is the actuating variable, e(t) is the error signal, and G(s) is the transfer function of the controller. KP, KI, and KD are the scaling factors for the P, I, and D parts of the controller. All three summands are calculated in parallel.
To match the specific requirements for quantum optics experiments, the standard PID controller was modified. The source code of the new implementation is free and can be downloaded from Ref. 67. Its structure is shown in Fig. 5. A master gain term, a second integral part, and an offset setting were added to the basic PID. The ranges of the control parameters were extended, and the original STEMlab signal generator application was integrated in the controller.
Block diagram of the modified STEMlab–PID in a feedback control loop.
A number of additional functions were implemented: the settings can be saved, the output voltage can be limited digitally, and an auxiliary input analog signal can be used to enable the controller. The latter function was implemented to allow automatic re-locking of the controller. For this the auxiliary intensity signal from the photodetector, “power diode” in Fig. 3 is used to discriminate between actual resonances of the PL TEM00-mode and parasitic signals. Moreover, an automatic integrator reset and a sample-and-hold mode were implemented. All additional functions are optional and run in parallel with the PID algorithm, without affecting its real-time performance.
To allow finer adjustment of the integrator settings, an optional decimation scheme was implemented, in which only every Nth sample (N = 2–1000) is summed.
D. Characteristics of the digital controller
1. Amplitude frequency response
Figure 6 shows the amplitude response with all ranges and step sizes of the control parameters, as well as their characteristic frequencies and slopes. The individually adjustable control parameters are the master gain G and the P-, I-, 2nd I-, and D-contribution of the control loop (Table III).
Amplitude response of the PID versus frequency. The ranges for the different control parameters are shown together with their effect on the response curve.
Amplitude response of the PID versus frequency. The ranges for the different control parameters are shown together with their effect on the response curve.
The physical meaning of the control parameters for the I- and D-part as 0 dB gain crossover frequencies becomes evident when we convert from angular to linear frequencies,
2. Limits
The closed-loop bandwidth of the digital controller is fundamentally limited by its group delay τ = 155 ns. The phase shift from the delay increases linearly with frequency and cannot be compensated in any real-time causal system. Its effect is visible in the Bode plots of Fig. 10. The additional phase shift is approximately −55.4° at 1 MHz, measured in-loop. This suggests a maximum closed-loop control bandwidth of a few MHz, determined by the total phase shift in the loop.
Another limitation exists due to the dynamic range bottlenecks at the ADC and the DAC. The effective resolution of the converters is lower than their nominal resolution of 14 bits. This was partially explained earlier for the DAC, where an excess digital noise contribution was observed. Figure 12 in the Appendix shows the voltage noise spectral density measured at the output of two STEMlab boards and two standalone boxes with and without an output amplifier. The integrated noise voltage within 1 Hz–1 MHz is measured to be 76.89 μV for the unmodified STEMlab board without amplifier stage and deactivated PID controller (A, red curve). With the offset resistors removed, 32.44 μV is measured (B, black curve). For the modified STEMlab installed inside the standalone box with activated PID controller, 120.12 μV is obtained (C, green curve). Finally, an integrated noise voltage of 392.66 μV is measured for the STEMlab board with offset resistors removed inside the standalone box and with amplifier stage added.
V. EXPERIMENT
During all experiments, the cavity was flushed with 99.95% pure oxygen at a rate of 0.022 l/min. A finesse Fmeas = 317 was measured with the Brewster-cut BBO crystal in place and mirrors with reflectivities of R1 = 98.50% and R2 = R3 = R4 = 99.95%.
A. Long term stability
For many applications in quantum optics, a powerstable UV source is required that reliably stays in lock over many hours. In order to demonstrate that our setup fulfills those requirements, it has been locked over 130 h during which the PL power and SHG power have been recorded every 1 s using two power meters (Thorlabs PM100D with S130VC Sensor and Ophir Vega with PD300-UV Sensor, cf. Fig. 3). Figure 7 shows the measured SHG output power scaled to the average SHG power (red) and to the instantaneous PL power (blue). The SHG power stayed on a level of about 186 mW, while the cavity was pumped with about 530 mW PL. The SHG power fluctuates by about ±7% peak-to-peak (3% rms) around its mean value on a few hour time scale, while the SHG power rescaled to the PL varies only by ±3% peak-to-peak (1% rms). This suggests that the SHG power fluctuations were mainly caused by pump power fluctuations, indicating that the alignment of the cavity did not change on relevant length or angular scales over the duration of the measurement.
Long-term SHG output power measurement. Red: measured SHG power, a power of “1” equals 186 mW. Blue: SHG normalized to the PL power squared as a measure of the conversion efficiency.
Long-term SHG output power measurement. Red: measured SHG power, a power of “1” equals 186 mW. Blue: SHG normalized to the PL power squared as a measure of the conversion efficiency.
To our knowledge, this is the longest published in-lock operation of a resonant SHG ring cavity. The output power of the SHG cavity fully agrees with the corresponding measurement and calculation presented in Ref. 23.
B. Lock performance
The performance of the cavity lock has been investigated by recording the power spectral density of the error signal generated by the differential photodiode (Fig. 3) using a vector signal analyzer (Agilent 89441A). Figure 11 shows the error signal for a loose lock achieved by choosing gain setting just sufficient enough to lock (blue) and for optimized settings (black). With the optimized lock, suppression by up to almost one order of magnitude for acoustic frequencies is achieved. A broad resonance at about 17 kHz sets the upper limit for the locking bandwidth in this setup (Fig. 11). This frequency is likely to be the first mechanical resonance of the small mirror, the piezoelectric actuator, and its holder. In a future setup, it could possibly be improved by applying the mount design presented in Refs. 68–70 since Eq. (10) results in an upper limit of f* ≈ 120 kHz for the piezoelectric actuator and mirror used here.
C. Acceleration sensitivity tests
The monolithic cavity described here is a prototype to be used for generating the clock transition probe light of a transportable Al+ clock to be used in relativistic geodesy campaigns. In order to prove the suitability of the doubling cavity for such an application, the entire breadboard was accelerated in the vertical direction while the cavity was locked. The acceleration was measured using an analog MEMS sensor (ADXL345) glued on the top lid of the cavity. Figure 8 shows the output power and vertical acceleration versus time. The measurement shows temporary SHG power drops of about 10% synchronous with the acceleration acting on the system. While an absolute acceleration of up to 1 g was obtained, the cavity stayed in lock and fully recovered the initial SHG power after the acceleration stopped.
Cavity SHG output power and acceleration versus time during acceleration excitation in the vertical direction. Red: SHG power, blue: vertical acceleration. For accelerations up to 1 g (gravity subtracted), the SHG output power fluctuates on the 10% level while the cavity stays in lock.
Cavity SHG output power and acceleration versus time during acceleration excitation in the vertical direction. Red: SHG power, blue: vertical acceleration. For accelerations up to 1 g (gravity subtracted), the SHG output power fluctuates on the 10% level while the cavity stays in lock.
In order to demonstrate the reliable operation of the system after transport, it was exposed to an acceleration profile typical for on road transportation. Therefore, the system was optimized first and then the SHG breadboard was placed on a multi-component acceleration exciter71 and exposed to a vibration power spectrum in the vertical direction according to ISO 13355:2016 with a total acceleration of 0.604 grms for 30 min. Afterwards, the SHG breadboard was reconnected to the rest of the system and locked without further realignment or optimization. The results are shown in Fig. 9. The obtained SHG power normalized to the pump power dropped to about 60%, while the SHG power normalized to the squared pump power circulating in the cavity remained almost constant. Since this indicates a drop in incoupling efficiency, the alignment of the mirrors M1 and M2 in Fig. 4 was optimized, resulting in approximately the same SHG power as before the acceleration exposure test. This strongly indicates that the alignment of the cavity itself was not substantially affected by the shaking test. Therefore the same test was repeated with five times the original acceleration (3.020 grms), also resulting in the same SHG output power after shaking and alignment correction of the mirrors M1 and M2. From these measurements, we conclude that the optical alignment of the cavity itself can withstand typical on road transport situations without any deterioration, while the mechanical stability of the input coupling optics and/or fiber connector for the PL needs to be improved.
SHG power before (left) and right after (middle) being exposed to a 30 min ISO 13355:2016 shaker test. Right: following optimization of the incoupling beam after the shaker test. Red: measured SHG power, blue: SHG normalized to the PL power squared, and pink: SHG normalized to the circulating PL power as a measure of the conversion efficiency.
SHG power before (left) and right after (middle) being exposed to a 30 min ISO 13355:2016 shaker test. Right: following optimization of the incoupling beam after the shaker test. Red: measured SHG power, blue: SHG normalized to the PL power squared, and pink: SHG normalized to the circulating PL power as a measure of the conversion efficiency.
VI. SUMMARY AND OUTLOOK
A mechanically stable monolithic enhancement cavity for SHG in the UV was demonstrated, including in-lock SHG power measurements during acceleration excitation. Less than 10% SHG power reduction during exposure of up to 1 g was observed with full recovery of the initial SHG power after acceleration stopped. It was shown that the cavity optical alignment can withstand 30 min of acceleration excitation with 3.020 grms. This is five times the acceleration amplitude as specified in the corresponding ISO 13355:2016 norm, demonstrating the suitability for transportable experiments. 130 h uninterrupted operation without decay in output power at 313 nm was demonstrated. During this time, the SHG power scaled to the pump power fluctuated by 1% rms. The basic design can easily be adapted to other resonator geometries in order to install crystals of different materials or shapes. The locking bandwidth of the setup presented here is 17 kHz. To obtain a higher locking bandwidth in the future, the design of the piezoelectrically actuated mirror mount demonstrated in Refs. 68–70 could be adopted. To reduce the risk of unwanted substances inside the cavity causing degradation of the crystal and/or mirrors in the next cavity generation, the housing can be left unanodized and metal sealings can be employed. For improved mechanical robustness and leak tightness, the crystal aligner can be made part of the cavity housing.
ACKNOWLEDGMENTS
We thank L. Klaus for performing the acceleration excitation measurements on the multi-component acceleration exciter and S. A. King for stimulating discussion. We acknowledge support from DFG through Grant No. CRC 1128 geo-Q, Project No. A03, CRC 1227 DQ-mat, and Project Nos. B03 and B06 and Leibniz Gemeinschaft through Grant No. SAW-2013-FBH-3.
APPENDIX: DETAILED SPECIFICATIONS AND PERFORMANCE CHARACTERISTICS OF THE DIGITAL CONTROLLER
STEMlab 125-14 analog front- and back-end specification.72
Inputs . | Outputs . | ||
---|---|---|---|
Number of channels: | 2 | Quantity: | 2 |
Bandwidth: | 50 MHz | Bandwidth: | 50 MHz |
Sample rate: | 125 MS/s | Sample rate: | 125 MS/s |
DAC resolution: | 14 bit | ADC resolution: | 14 bit |
Full scale voltage: | Selectable: ±1 V, ±20 V | Full scale power: | 9 dBm |
Minimal voltage sensitivity: | ±0.122 mV, ± 2.44 mV | 1 LSB: | 0.122 mV |
Input impedance: | 1 MΩ || 10 pF | Load impedance | 50 Ω |
(for full scale power): | |||
DC offset error: | <5% | DC offset error: | <5% |
Gain error: | LV: <3%, HV: <10% | Gain error: | <5% |
Input noise level: | −119 dBm/Hz | Slew rate limit: |
Inputs . | Outputs . | ||
---|---|---|---|
Number of channels: | 2 | Quantity: | 2 |
Bandwidth: | 50 MHz | Bandwidth: | 50 MHz |
Sample rate: | 125 MS/s | Sample rate: | 125 MS/s |
DAC resolution: | 14 bit | ADC resolution: | 14 bit |
Full scale voltage: | Selectable: ±1 V, ±20 V | Full scale power: | 9 dBm |
Minimal voltage sensitivity: | ±0.122 mV, ± 2.44 mV | 1 LSB: | 0.122 mV |
Input impedance: | 1 MΩ || 10 pF | Load impedance | 50 Ω |
(for full scale power): | |||
DC offset error: | <5% | DC offset error: | <5% |
Gain error: | LV: <3%, HV: <10% | Gain error: | <5% |
Input noise level: | −119 dBm/Hz | Slew rate limit: |
Limit values of the control parameters.
Parameter . | Maximum . | Step size . | Minimum frequency . | Maximum frequency . | Inaccuracy (%) . |
---|---|---|---|---|---|
Gain: | 10 | 0.01 | … | … | <10 |
Kp: | 10 | 0.01 | … | … | <10 |
Ki: | 80 Hz | 573 kHz | <17 | ||
Kd: | 58 ns | 0.01 ns | 2.74 MHz | 16 GHz | <6 |
Parameter . | Maximum . | Step size . | Minimum frequency . | Maximum frequency . | Inaccuracy (%) . |
---|---|---|---|---|---|
Gain: | 10 | 0.01 | … | … | <10 |
Kp: | 10 | 0.01 | … | … | <10 |
Ki: | 80 Hz | 573 kHz | <17 | ||
Kd: | 58 ns | 0.01 ns | 2.74 MHz | 16 GHz | <6 |
Bode diagram open loop. From top to bottom: proportional/integral/differential part.
Bode diagram open loop. From top to bottom: proportional/integral/differential part.
Error signal (=STEMlab input) power spectral density for two different PI settings. A (black): lock with optimized gain settings, B (blue): low-gain lock.
Error signal (=STEMlab input) power spectral density for two different PI settings. A (black): lock with optimized gain settings, B (blue): low-gain lock.
Output noise (=STEMlab output) power spectral density with terminated input. A (red): unmodified STEMlab without amplifier stage, PID deactivated. B (black): the same as A, but with offset resistors removed. C (green): STEMlab in the box, PID enabled, without amplifier, Kp = 1. D (blue): the same as C, but with amplifier (8×) and 100 KHz low-pass filter. Integrated output voltage noise (bandwidth 1 Hz–1 MHz): A = 76.89 μV, B = 32.44 μV, C = 120.12 μV, and D = 392.66 μV.
Output noise (=STEMlab output) power spectral density with terminated input. A (red): unmodified STEMlab without amplifier stage, PID deactivated. B (black): the same as A, but with offset resistors removed. C (green): STEMlab in the box, PID enabled, without amplifier, Kp = 1. D (blue): the same as C, but with amplifier (8×) and 100 KHz low-pass filter. Integrated output voltage noise (bandwidth 1 Hz–1 MHz): A = 76.89 μV, B = 32.44 μV, C = 120.12 μV, and D = 392.66 μV.