thanks: A footnote to the article title

Radio-frequency cavity field measurements through free falling bead

Xiaonan Du    Lars Groening GSI Helmholtzzentrum für Schwerionenforschung GmbH, Darmstadt, Germany
(June 27, 2024)
Abstract

This paper introduces a novel bead-falling measurement method for the precise and efficient mapping of electromagnetic fields within radio-frequency (RF) cavities, which are crucial components in the design of accelerators. The traditional bead-pull method, while effective, involves mechanical complexities and is prone to errors from wire perturbations. The innovative method reported here leverages the simplicity and accuracy of free-falling beads to mitigate these issues. This technique eliminates the need for a wire-pulley system, thereby simplifying the experimental setup and reducing potential mechanical errors. We detail the development and operational principles of this new method, including the design of a compact, portable measurement device that integrates a bead/droplet release system and a bead detection system linked to a Vector Network Analyzer (VNA). The device has been tested with a three-gap buncher cavity and a scaled Alvarez-type cavity, demonstrating its ability to perform rapid, reliable field measurements under various challenging conditions, including low signal-to-noise ratios and environmental vibrations. The results confirm the method’s superiority in precision and operational efficiency, potentially setting a new standard for RF cavity diagnostics and tuning.

preprint: APS/123-QED

I Introduction

In the evolving landscape of accelerator physics and the design of resonant cavities, precise measurement and tuning of electromagnetic fields of rf-cavities are paramount. The bead-pull method, a technique steeped in the principles of Slater’s small signal perturbation theory [1, 2], emerges as a cornerstone in this domain. By introducing a small bead into a resonant cavity, this method allows for the quantification of shifts in the cavity’s resonant frequency, attributing these shifts to the perturbative presence of the bead. Such frequency variations are indicative of the electric and magnetic field strengths at the bead’s location, offering insights into the field distribution within the cavity.

In a typical bead-pull setup, X-Y linear motion systems are placed on each side for aligning the axis of the wire (fishing line) with the cavity axis. A motorized linear actuator with a specified stroke is used to move the bead along the measurement path inside the cavity by pulling the wire [3]. In some instances, an optical encoder provides feedback for the bead’s location, which is calibrated to ensure precise positioning and movement.

This method not only underpins the characterization of longitudinal and transverse electromagnetic fields within cavities but also facilitates the iterative tuning required to achieve desired field distributions. Its application spans the assessment of cavities crafted from various materials, employing objects like spheres and needles to invoke perturbations. The advent of automated test stands, notably those developed at prestigious institutions like CERN, marks a significant advancement in the method’s implementation. These test stands, equipped with motors to maneuver the bead across three axes, enable comprehensive field mapping, ensuring the cavity’s optimization for its intended application.

The bead-pull technique’s relevance is underscored by its role in the manufacturing and tuning of high-frequency resonant cavities, a critical component of the burgeoning accelerator industry. With the demand for these cavities on the rise, the method’s efficacy in facilitating the fine-tuning of electromagnetic field distributions becomes increasingly important. The common measurement procedure for Radio-Frequency Quadrupoles (RFQs) consists of measuring each quadrant and comparing field flatness along the length of the structure. Obtaining the desired field distribution is an iterative process. bead pull measurements were a continuous part of the process, critical for achieving the RFQ’s precise electromagnetic field requirements [4, 5].

However, the bead-pull method is not without its limitations. The wire and pulley system typically requires customization for the cavity to be measured, and its final dimensions depend on the cavity under test. The need for mechanical adjustments, dealing with potential issues like bead vibration or wire over-latch, and ensuring stable pulling can make the process more time-consuming and labor-intensive. Additionally, the wire can cause perturbations, particularly concerning the potential interference of the wire system with certain (dipole) modes of the cavity, thereby rendering the proportion of higher-order field components in the measured field distribution inaccurate [6].

In pursuit of an innovative approach, we have been developing what we term the bead-falling-measurement method. This technique leverages the free fall of a bead to induce perturbations, thereby eliminating the need for a pulley-wire system and its associated complications. This method presents a promising alternative to the conventional bead-pull measurement, offering an avenue for further refinement and efficiency in cavity diagnostics and tuning. As the field of accelerator technology continues to advance, the evolution of measurement techniques like the bead-pull and bead-falling methods will play a critical role in shaping the future of resonant cavity measurement and tuning.

The next section illustrates the method’s working principle and outlines its advantages towards the conventional technique. It is followed by description of the applied from-the-shelf hardware. The fourth section is dedicated to the applied software. Measurements done with a buncher cavity and a model of an accelerating cavity are presented afterwards. Final conclusions evaluate the method’s capabilities and its further potential. Two appendices are on issues related to free-fall including air resistance and to the reduction of noise.

II conceptual design and advantages

Similar to bead-pull measurements, our measurement also requires the use of a bead to create perturbations and employs a Vector Network Analyzer (VNA) to capture Radio-Frequeny (RF)-measurement data. Since the bead is in free-fall through the cavity, the acceleration causes the measured perturbation-time curve to not linearly translate to a perturbation-position curve. Therefore, a setup is needed to release the bead steadily, and then through measurement and calculation, correlate the time-position relationship accurately. Figure 1 shows the conceptual design of this measurement method, and subsequent sections will elaborate on the implementation details.

Refer to caption
Figure 1: Scheme of bead-fall RF-measurements. The cavity is connected to two of the VNA ports. The device includes the bead/droplet release system and the bead detector. It is positioned such that the falling bead is aligned to the beam axis of the cavity. Beads are released from the top and fall freely through the bead-detector and subsequently through the cavity.

The new method overcomes several technical challenges inherent in traditional bead-pull setups and brings numerous benefits listed below:

  • well-defined perturbation position: The bead can be very small (e.g., 0.5 to 1 mm in diameter), ensuring the perturbation position to be precisely defined. This allows the bead to pass through small gaps, like those between RFQ vane tips [5] or X-band cavities [7] without significant field variation over its size.

  • no perturbation from wire: Perturbation from the nylon (ϵr=4subscriptitalic-ϵ𝑟4\epsilon_{r}=4italic_ϵ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = 4) wire is inevitable in the bead-pull method, though it is not an issue using the new method.

  • fast measurements: With the free-falling duration of less than a second, measurements can be conducted rapidly. This enables real-time field monitoring at about 2 Hz, allowing tuners to be adjusted based on continuously updated field distributions, thus making the cavity tuning process more efficient.

  • material versatility: use of different materials for the perturbation, such as liquid drops, facilitates fast iterative bead release, adding flexibility to the measurement process.

  • well-defined path: The bead’s path is determined by the launching device’s position and velocity, which can be easily controlled, eliminating the need for a wire pulling device and, thus, avoiding wire and bead vibration as well as perturbation from materials like nylon.

  • simplified setup: The absence of a wire and pulley system not only simplifies the experimental setup but also significantly saves time during preparation.

  • accessible measurements in complex locations: The possibility of the bead’s path being a parabola, may allow for measurements in cavity locations not accessible by traditional bead-pull measurements.

Integrating these advantages, the new method streamlines the experimental setup, enhances measurement speed and accuracy, offers flexibility in measurement, reduces mechanical interference, and enables real-time monitoring. These improvements make the cavity tuning process significantly more efficient and adaptable to varying experimental conditions.

One of the most promising applications of this device is the tuning of RF-cavities. For multi-gap cavities, such as RFQs and Drift Tube Linacs (DTLs), a complex field distribution tuning process is always necessary to address deviations in the actual electric field distribution caused by mechanical machining errors from the design value (typically equal electric field peaks in all cells). This tuning process usually requires multiple field distribution measurements using the bead pull measurement technique, thereby iteratively adjusting the field distribution closer to the design value.

First, the device introduced in this paper can significantly simplify the experimental preparation phase, as it integrates all necessary functionalities into a small box, requiring only a single alignment on top of the cavity. Second, it enables each field distribution measurement to be completed in a very short time (less than 0.5 s), and continuous measurements can be conducted at these intervals to achieve real-time monitoring of the field distribution. Any changes in the field distribution caused by moving the tuner can be immediately presented to the experimenters. This can greatly accelerate the entire field tuning process.

III Hardware

After continuous exploration, experimentation, and iterations, the design of the first measurement device has been finalized, which utilizes a series of components as shown in Fig. 2. These components are easily accessible, inexpensive, and mature products.

Refer to caption
Figure 2: Overview of required components.

Figure 3 showcases the device formed after integrating all the components. It is evident that the entire device’s volume is only half the size of a letter-sized sheet of paper, making it portable and lightweight compared to the traditional wire-pulley system, and it is universal for all kinds of cavities, except the cavity is not suitable for vertical standing. For the issue of alignment, the current solution is to use a 3D-printed cap (designed according to the cavity) to position the bead. The positioning accuracy depends on the precision of 3D printing, being ensured for typical measurements. The following subsection provides detailed decriptions of each hardware component individually.

Refer to caption
(a) Inner structure overview of the device.
Refer to caption
(b) The device is placed on top of the cavity to be measured, connected to a VNA and a monitor.
Figure 3: All components are fixed to the 3D printed case (a). Measurement device installed on top of the cavity to be measured (b).

III.1 bead release system

The perturbator needs to be released with the release point precisely positioned relatively to the chamber (cavity) being tested in both, height and horizontal position. When releasing the perturbator, it is crucial to ensure that it does not acquire any horizontal velocity that could affect its falling path, nor should it have any random initial vertical velocity that could cause deviations in the time-position relationship. For the two types of perturbators, solid beads and liquid droplets, different release mechanisms should be designed. As a starting point, a medical drip infuser has been used to generate continuous droplets of water as perturbators. The blunt tip syringe needles of the infuser can smoothly produce water droplets in a manner that meets the requirements, with the size of the droplets controllable by selecting the size of the needle, and the rate of droplet formation easily regulated using a valve. Additionally, we have designed and 3D printed a fixture part to conveniently secure the needle to the exterior of the device.

III.2 linear motion system

The current design utilizes a linear motion system composed of a roller screw to achieve one-dimensional movement. If necessary, this system can be expanded to a two-dimensional motion system. In some cases, where only the electromagnetic field along a single path of the beam axis needs to be measured, the linear motion system can be omitted.

III.3 bead detection system

A laser diode and a light sensor are positioned directly opposite to each other and are separated by about two centimeters. They form a photoelectric barrier, refered to as trigger location. The setup is arranged such that the laser diode emits a beam that travels across the gap and is received by the light sensor. Without bead, the light sensor outputs a high voltage level indicating the uninterrupted reception of the laser beam.

The bead is released about 2 cm from above the gap. Prior to the bead reaching the gap, a digital light intensity sensor module, which can switch between Binary Output (BO) and Analog Output (AO), continuously outputs ”1” due to the unobstructed laser light. When the bead reaches the gap and obstructs the laser beam, the amount of light reaching the light sensor decreases below a predefined threshold, causing it to output ”0”.

The Raspberry Pi is programmed to monitor the output from the light sensor in real-time. Upon detecting a transition to a low voltage level, indicative of the bead obstructing the laser beam and thus being at the gap position, the Raspberry Pi triggers a VNA also connected to it to commence measurements. This setup allows for precise timing and initiation of measurements based on the physical interruption of a laser beam by a falling object. The total cost of the hardware mentioned so far is below300US$300𝑈��currency-dollar~{}300~{}US\$300 italic_U italic_S $.

III.4 computer (Raspberry PI 4B)

The system employs the compact Raspberry Pi 4B (RPI), both as a controller for the VNA and as a computing unit capable of running software for data collection, processing, and visualization. Thanks to its small size (85 mm×\times×56 mm), the Raspberry Pi can easily be integrated into a portable device.

III.5 Vector Network Analyzer

A crucial point to note during measurements is that the VNA should have a sufficiently high (Intermediate frequencies) IF bandwidth. To avoid undersampling issues, it is necessary to collect enough data points during the short time the perturbator falls through the cavity, especially since the falling speed increases towards the rear end of the cavity, resulting in fewer data points being allocated. The experiments utilized the Rohde &\&& Schwarz ZNB-4 VNA, which has a maximum IF bandwidth of 10,000, ensuring adequate sampling. For situations where the falling speed is exceptionally fast or the VNA’s IF bandwidth is low, leading to undersampling, an alternative solution is available. This involves repeating the falling measurements multiple times, each time initiating the VNA sweep with a slightly different delay to offset the sampling points. Combining the data points from these multiple measurements can resolve the issue of undersampling.

IV software

The software component of the newly developed device is designed to facilitate a seamless interaction between the user and the measurement content and process through its user interface (UI), as detailed in Fig. 4. This software is pivotal in orchestrating the device’s operations, offering several key functionalities that enhance the measurement capabilities and data processing efficiency.

Refer to caption
Figure 4: The user interface of the software.

Firstly, the software allows users to customize measurements via input through the UI, tailoring the process to meet specific requirements. It engages in communication with a VNA to dispatch measurement instructions and retrieve data, thereby streamlining the measurement workflow. A significant feature includes the conversion of phase difference-time curves into electric field strength-position curves, which are then visualized through plots to present the measurement outcomes comprehensively.

Additionally, the software is capable of monitoring optical detectors to detect signals triggered by the bead passing through trigger points, subsequently prompting the VNA to initiate a measurement instantly. It also controls the linear motion system to adjust the bead release position, ensuring precise manipulation of the measurement environment.

A focal point of the software’s design is its ability to automate repeated measurements and process the data efficiently. This capability is critical for rapidly acquiring a large volume of data within a short timeframe, enabling the use of Singular Value Decomposition (SVD) for noise filtering. Such an approach not only significantly enhances the quality of the measured electric field distribution but also allows for noise analysis. The software’s automation feature facilitates continuous droplet release, inducing perturbations that enable real-time updates of the field distribution measurements. This is particularly beneficial for the tuning process of cavities, where the ability to quickly adjust and measure changes is crucial. The flow chart of the software is shown in Fig. 5.

InitailVNA conectordinary frequency and Q measuretimer<<<120s?wait until bead detected\pgfmathresultptdo frequency track,reset timersend start command to VNAsuccess?\pgfmathresultptplot on UI, append data\pgfmathresultptNmaxsubscript𝑁𝑚𝑎𝑥N_{max}italic_N start_POSTSUBSCRIPT italic_m italic_a italic_x end_POSTSUBSCRIPT reachedor stop buttonpushed?stand byyesnoyesnonoyes
Figure 5: Flow chart of the measuring procedure.

Moreover, the software incorporates additional considerations to ensure measurement accuracy and reliability. It is programmed to re-measure the resonance frequency at set intervals (e.g., every 120 seconds) to mitigate the impact of frequency drift on the results. It also includes mechanisms to identify and discard measurements that are compromised due to delays or timeouts in commands sent to the VNA, which could render the results unusable. This level of meticulousness in design and functionality underscores the software’s integral role in enhancing the device’s performance and reliability in complex measurement scenarios.

V measurements on a buncher cavity

A first set of measurements has been done using a three-gap buncher cavity intended to be installed into an existing accelerator section. It is depicted in Fig. 6 and its parameters are listed in Table 1.

Refer to caption
Figure 6: Three-gap buncher cavity of spiral-type, to be operated at 36 MHz.
Table 1: Parameters of the buncher BB3
Parameter Unit Value
RF-frequency MHz 36.136
Gaps # 3
Gap length mm 12.7–13.5
Drift tubes # 2
Drift tube length mm 38.1–40.8
Aperture mm 100
Tank diameter mm 420
Tank length mm 525.7
Q-factor 1000

The cavity’s orientation aligns the beam axis w.r.t. local gravity, with the measurement apparatus mounted atop the cavity along the beam axis. Due to low Q value (1000) of this buncher with long spirals, in the measured data, the magnitude of noise in the measurement data is close to the magnitude of the signal itself. Though fast iterative measurements allow to collect big mount of data within short time.

Before the device is put into application, two fixed parameters need to be determined. The first one is the speed v0subscript𝑣0v_{0}italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT of the bead when passing the trigger location. The second parameter to be determined is the system delay tdsubscript𝑡𝑑t_{d}italic_t start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT, designated as, which is the time delay between the bead passing the trigger location and the VNA generating the first measurement point. A cavity with three equidistant gaps is suitable for calibrating these two parameters. The distance between the three gaps of the cavity is known and taken as a constant L𝐿Litalic_L. Conducting a measurement will produce a curve with three distinct peaks corresponding to the three gaps. The speed of the falling bead at all measurement points can be inferred from the distances between these three peaks, thereby deducing the speed of the bead at the trigger location and hence tdsubscript𝑡𝑑t_{d}italic_t start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT. The bead is dropped from different heights into the cavity at different speeds such that more reliable determinations of v0subscript𝑣0v_{0}italic_v start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and tdsubscript𝑡𝑑t_{d}italic_t start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT can be obtained by averaging.

Figure 7 depicts an example for such measurements. After these two parameters have been determined, the device is ready for measurements on any other type of cavities without needing to repeat the aforementioned process. In the bead falling measurement, the distance between the device and the cavity entrance is a known constant. Therefore, the measured time-phase shift curves can be translated into position-field strength curves.

Refer to caption
Figure 7: Phase shift data measured from the bead launched from four different heights L0subscript𝐿0L_{0}italic_L start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, L1subscript𝐿1L_{1}italic_L start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, L2subscript𝐿2L_{2}italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and L3subscript𝐿3L_{3}italic_L start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT into the buncher cavity.

The noisy data set obtained from the perturbation method is compiled into 1D arrays, making it suitable for processing with SVD. It should be noted that SVD is performed only after multiple iterations and sufficient data collection; it does not play a role in the real-time field distribution monitoring at 2 Hz mentioned earlier. Therefore, to remove noise from the updates of real-time field distribution, other algorithms should be considered. However, before applying SVD, the measured frequency variation data must be converted into phase shift data by subtracting it with a reference phase. The reference phase for each iteration is determined by taking the average phase of the points near both boundaries. Each iteration provides a 1D array of 3000 points corresponding to time steps. After 500 iterations, a 2D data array (500×30005003000500\times 3000500 × 3000) is obtained. Figure 8 shows a summarizing plot.

Refer to caption
Figure 8: Measured data taken with the buncher cavity with 3000 discrete time points within 0.4 second each, corresponding to a length of about 1.2 m along the cavity axis.

After applying SVD (see Appendix B), the data is decomposed into several components, with the first three singular vectors shown in Fig. 9. The first singular vector is essentially considered to be the noise-free essence of the data. It is treated and compared with simulated field distributions.

Refer to caption
Figure 9: First three singular vectors obtained from SVD (buncher cavity), ordered by the magnitude of their corresponding singular values.

Comparison between measured field and Micro-Wave-Studio (MWS) simulations [8] is shown in Fig. 10. There are some deviations as a result of the finite accuracies of the machining, eigenmode position definition, and the measurement itself.

Refer to caption
Figure 10: Measured electric field(denoted by E^abssubscript^𝐸𝑎𝑏𝑠\hat{E}_{abs}over^ start_ARG italic_E end_ARG start_POSTSUBSCRIPT italic_a italic_b italic_s end_POSTSUBSCRIPT along the axis of the buncher cavity compared with simulation.

Afterwards, 2D field measurements have been performed on the same cavity. The positioning of the bead is achieved using a 3D-printed cap that fits the entrance of the cavity, with a slit constraining the position of the droplets. The length of the slit covers the off-axis area intended for measurement, and its width is just sufficient to allow a water droplet to pass through. In each iteration, the needle is moved a small constant distance, and a 1D field along the path of free falling is measured. As the needle moves from one end of the slit to the other, the measurement of the entire 2D electric field is completed. Figure 11 show the measured 2D array data in the form of a bitmap. This data is the phase shift as a function of time, measured at different off-axis distances. it is converted to 2D field distribution, according to time-height relation presented in appendix A.

Refer to caption
Figure 11: Original 2D array data from the buncher cavity with even time interval of dt=0.4s/3000𝑑𝑡0.4𝑠3000dt=0.4s/3000italic_d italic_t = 0.4 italic_s / 3000.

Then, with proper translation and scaling, the measured field(denoted by E^abssubscript^𝐸𝑎𝑏𝑠\hat{E}_{abs}over^ start_ARG italic_E end_ARG start_POSTSUBSCRIPT italic_a italic_b italic_s end_POSTSUBSCRIPT) can be compared with corresponding field obtained from simulation(denoted by Eabssubscript𝐸𝑎𝑏𝑠E_{abs}italic_E start_POSTSUBSCRIPT italic_a italic_b italic_s end_POSTSUBSCRIPT), as shown in Figure 12.

Refer to caption
Figure 12: The distribution E^abssubscript^𝐸𝑎𝑏𝑠\hat{E}_{abs}over^ start_ARG italic_E end_ARG start_POSTSUBSCRIPT italic_a italic_b italic_s end_POSTSUBSCRIPT represents the experimentally measured two-dimensional field, while Eabssubscript𝐸𝑎𝑏𝑠E_{abs}italic_E start_POSTSUBSCRIPT italic_a italic_b italic_s end_POSTSUBSCRIPT denotes the corresponding field obtained from simulations. Both fields are normalized to their respective maxima, specifically the peak value observed at the middle gap. The term E^absEabssubscript^𝐸𝑎𝑏𝑠subscript𝐸𝑎𝑏𝑠\hat{E}_{abs}-E_{abs}over^ start_ARG italic_E end_ARG start_POSTSUBSCRIPT italic_a italic_b italic_s end_POSTSUBSCRIPT - italic_E start_POSTSUBSCRIPT italic_a italic_b italic_s end_POSTSUBSCRIPT quantifies the deviation of the measured values from the simulated ones.

VI measurement on a 1:3 aluminum model cavity of Alvarez type

The measurement process was similarly applied to a model cavity, specifically a 1:3 scale aluminum replica of a 10-gap Alvarez-type cavity. This model primarily served to explore an efficient stabilization strategy, focusing on tilt sensitivity [9, 10]. It is depicted in Fig. 13 and its parameters are listed in Table 13.

Refer to caption
Figure 13: Scaled (1:3) aluminum model of an Alvarez-type cavity.
Table 2: Parameters of the scaled (1:3) Alvarez-type cavity model.
Parameter Unit Value
RF-frequency MHz 325.224
Gaps # 10
Gap length mm 12.7–13.5
Drift tubes # 9
Drift tube length mm 38.1–40.8
Drift tube diameter mm 60.0
Aperture mm 10.0
Tank diameter mm 634.5
Tank length mm 525.7
Q-factor 500

A phase shift measurement was chosen rather than resonant frequency shift as the phase is significantly easier to detect. The measurement exhibited a very low signal-to-noise ratio, approaching 2:1, due to three main reasons: 1) The cavity’s Q-value of 530 was low because it was made from aluminum. 2) The testing environment included background noise from surrounding equipment, such as vacuum pumps, contributing to the overall noise level. 3) The coupling loop and pickup were made of loose copper wires, which, under the influence of background mechanical vibrations, created relative vibrations with the cavity. These adverse factors are usually strenuously avoided in conventional bead pull measurements. However, our bead-fall method and the corresponding device made it possible to conduct measurements in such an extremely unfavorable environment. Figure 14 plots the measured phase shifts.

Refer to caption
Figure 14: Measured phase shift data from the DTL model cavity, with a signal-to-noise ratio of approximately 2. Different colors represent multiple phase shift measurements conducted using the DTL. Only five representative results from thousands are presented.

Despite the initial data being heavily laden with noise, attributed to the soft coupler, we successfully discerned the field distribution. The ability to perform rapid, iterative measurements facilitated the accumulation of a substantial dataset in a brief timeframe. This wealth of data enabled us to employ SVD techniques to filter out noise and accurately ascertain the field distribution, as depicted in Fig. 15.

Refer to caption
Figure 15: Measured electric field along the beam axis of the DTL model cavity.

Alvarez-type cavities are characterized by heavy drift tubes suspended by support rods from the upper side of a horizontally positioned cavity cylinder. However, during measurements, the cavity is placed vertically, leading to gravity causing the support rods to bend downwards, which in turn leads to unpredictable changes in the distance between the drift tube gaps. Additionally, the cavity was not precisely aligned. These factors resulted in a measured field distribution that could not be directly compared with theoretical values. Nonetheless, the measured field distribution clearly displayed the gap fields of ten cells, providing a reliable basis for tuning the field flatness.

VII conclusion

This study introduces an innovative bead-falling measurement technique, establishing a groundbreaking alternative to traditional RF-field mapping methods, particularly the bead-pull approach. The technique is underscored by its streamlined operational efficiency, reducing mechanical complexity and facilitating rapid, successive measurements. Employing SVD for meticulous noise filtration, this method delineates precise field distributions from noisy data, demonstrating exceptional accuracy.

A key highlight of our research is the development of a new measurement device characterized by its compactness, portability, and universal applicability across a diverse range of cavity types. This device, integral to the bead-falling measurement technique, has been rigorously tested under challenging conditions, including low signal-to-noise ratios and environmental vibrations. Its robust design and operational resilience were evident in the accurate discernment of field distributions in both a three-gap buncher cavity and a 10-gap Alvarez-type cavity model, confirming its efficacy in adverse settings.

The device’s minimalist design simplifies the preparatory phase of experiments, significantly reducing setup times and facilitating swift measurement cycles. This advancement not only enhances the practicality of field distribution measurements but also extends the scope of application to various cavity configurations without the need for extensive modifications or customization. The measurement process’s efficiency, coupled with the device’s versatility, underscores the potential for widespread adoption in accelerator physics and resonant cavity diagnostics.

Looking ahead, we aim to further refine this technique and its associated hardware, broadening its application spectrum and leveraging its capabilities to optimize cavity diagnostics and tuning. The introduction of the bead-falling measurement method and its accompanying device offers a potent tool for the precise evaluation and enhancement of electromagnetic field distributions in resonant cavities.

Appendix A consideration of air resistance during free fall

A spherical bead with the cross section area A𝐴Aitalic_A which falls at velocity of |v|𝑣|v|| italic_v | experiences the friction force from air having the density ρ𝜌\rhoitalic_ρ=1.2 kg/m33{}^{\text{3}}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT of

Fair=12ρACd|v2|,subscript𝐹𝑎𝑖𝑟12𝜌𝐴subscript𝐶𝑑superscript𝑣2F_{air}=\frac{1}{2}\rho AC_{d}|\vec{v}^{2}|\,,italic_F start_POSTSUBSCRIPT italic_a italic_i italic_r end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_ρ italic_A italic_C start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT | over→ start_ARG italic_v end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT | , (1)

where Cdsubscript𝐶𝑑C_{d}italic_C start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT is the dimensionless drag coefficient being equal to 0.47 for spheres. The falling speed as a function of time reads as

v(t)=2mgρACdtanh(tgρCdA2m).𝑣𝑡2𝑚𝑔𝜌𝐴subscript𝐶𝑑𝑡𝑔𝜌subscript𝐶𝑑𝐴2𝑚v(t)=\sqrt{\frac{2mg}{\rho AC_{d}}}\tanh{(t\sqrt{\frac{g\rho C_{d}A}{2m}})}\,.italic_v ( italic_t ) = square-root start_ARG divide start_ARG 2 italic_m italic_g end_ARG start_ARG italic_ρ italic_A italic_C start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT end_ARG end_ARG roman_tanh ( italic_t square-root start_ARG divide start_ARG italic_g italic_ρ italic_C start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT italic_A end_ARG start_ARG 2 italic_m end_ARG end_ARG ) . (2)

Defining k:=12ρACdassign𝑘12𝜌𝐴subscript𝐶𝑑k:=\frac{1}{2}\rho AC_{d}italic_k := divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_ρ italic_A italic_C start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT, the height as a function of time turns into

h(t)=h0+mln(coshgkmt)k𝑡subscript0𝑚𝑔𝑘𝑚𝑡𝑘h(t)\,=\,h_{0}\,+\,\frac{m\ln{(\cosh{\sqrt{\frac{gk}{m}}t})}}{k}italic_h ( italic_t ) = italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + divide start_ARG italic_m roman_ln ( roman_cosh square-root start_ARG divide start_ARG italic_g italic_k end_ARG start_ARG italic_m end_ARG end_ARG italic_t ) end_ARG start_ARG italic_k end_ARG (3)

where h0subscript0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the initial height, m𝑚mitalic_m is the spheres mass, and g𝑔gitalic_g is the acceleration from gravity. With advancing time, the velovity converges towards the stable velocity of

vs=2mgρCDA.subscript𝑣𝑠2𝑚𝑔𝜌subscript𝐶𝐷𝐴v_{s}=\sqrt{\frac{2mg}{\rho C_{D}A}}.italic_v start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = square-root start_ARG divide start_ARG 2 italic_m italic_g end_ARG start_ARG italic_ρ italic_C start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT italic_A end_ARG end_ARG . (4)

Within relevant heights of few meters, translations from the coriolis force can be neglected.

Appendix B Singular value decomposition in data process

SVD, or Singular Value Decomposition, is a mathematical technique used in linear algebra to decompose a matrix into three distinct matrices, specifically designed to identify the intrinsic geometric structure of the data. The decomposition is expressed as:

A=USVT,𝐴𝑈𝑆superscript𝑉𝑇A=USV^{T}\,,italic_A = italic_U italic_S italic_V start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT , (5)

where A𝐴Aitalic_A is the original matrix, U𝑈Uitalic_U is an m×m𝑚𝑚m\times mitalic_m × italic_m orthogonal matrix whose columns are the left singular vectors of A𝐴Aitalic_A, S𝑆Sitalic_S is an m×n𝑚𝑛m\times nitalic_m × italic_n rectangular diagonal matrix with non-negative real numbers on the diagonal, known as the singular values of A𝐴Aitalic_A, VTsuperscript𝑉𝑇V^{T}italic_V start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT (the transpose of V) is an n×n𝑛𝑛n\times nitalic_n × italic_n orthogonal matrix whose columns are the right singular vectors of A𝐴Aitalic_A. The singular values in are typically arranged in descending order. This ordering is crucial because the first few singular values (and their corresponding singular vectors) often capture the most significant modes of variation in the data, whereas the smaller singular values tend to correspond to noise or less important details.

One of the most significant advantages of SVD in the context of noisy data is its ability to perform noise reduction or data smoothing. By truncating the S𝑆Sitalic_S matrix to retain only the largest singular values (and ignoring the smaller ones that often correspond to noise), SVD can reconstruct a version of the original matrix that preserves the most critical structures while filtering out noise. This process is known as low-rank approximation.

SVD is particularly useful in signal processing and statistics for noise reduction. By identifying the singular values, which are essentially the strengths of the signals within the data, it allows for the separation of the ”signal” (important information) from the ”noise” (unimportant information). This is because the singular values corresponding to noise tend to be significantly smaller than those corresponding to the signal. By truncating or ignoring the smaller singular values (and their corresponding vectors), we can approximate the original data with less noise. Details including successfull application of SVD to simulated RF-cavity fields can be found in [11].

Appendix C Coriolis force

The Coriolis force is an apparent force that arises from the rotation of the Earth. It is an inertial force described mathematically by the Coriolis effect. To understand it intuitively, consider that the Earth rotates from west to east; this rotation influences the movement of air, water, and even objects in motion across its surface.

When an object falls freely (like the bead or droplet dropped from a height), it not only moves downward due to gravity but also moves horizontally due to its initial horizontal component of velocity (if any) and due to Earth’s rotation. Here’s how the Coriolis force affects it:

The magnitude of this deflection depends on several factors: The height from which the object falls (greater height means more time to be deflected). The latitude (effect is maximal at the poles and zero at the equator). The object’s speed of descent.

Consider a freely falling object at a latitude ϕitalic-ϕ\phiitalic_ϕ. The Coriolis acceleration acsubscript𝑎𝑐a_{c}italic_a start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT can be approximated by the formula:

ac=2Ωvsin(ϕ)subscript𝑎𝑐2Ω𝑣italic-ϕa_{c}=2\Omega v\sin(\phi)\,italic_a start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = 2 roman_Ω italic_v roman_sin ( italic_ϕ )

where: ΩΩabsent\Omega\approxroman_Ω ≈7.27 \cdot 10-5-5{}^{\text{-5}}start_FLOATSUPERSCRIPT -5 end_FLOATSUPERSCRIPT rad/s is the Earth’s angular velocity, v𝑣vitalic_vis the velocity of the falling object, and ϕitalic-ϕ\phiitalic_ϕ is the latitude at which the object is falling.

For an object starting from rest and falling a distance hhitalic_h, the time t𝑡titalic_t of fall can be calculated using:

t=2hg.𝑡2𝑔t=\sqrt{\frac{2h}{g}}\,.italic_t = square-root start_ARG divide start_ARG 2 italic_h end_ARG start_ARG italic_g end_ARG end_ARG .

The acceleration acsubscript𝑎𝑐a_{c}italic_a start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT leads to a deflection which can be integrated over time to find the eastward displacement due to the Coriolis force alone. This deflection is given by

d=12act2.𝑑12subscript𝑎𝑐superscript𝑡2d=\frac{1}{2}a_{c}t^{2}\,.italic_d = divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_a start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT .

Considering the latitude ϕitalic-ϕ\phiitalic_ϕ at the GSI laboratory of 50oo{}^{\text{o}}start_FLOATSUPERSCRIPT o end_FLOATSUPERSCRIPT and assuming a launching height of 2 m, the eastward displacement is less than 0.15 mm, being negligible.

References

  • Maier and Slater [1952] L. C. Maier and J. C. . Slater, Field strength measurements in resonant cavities, Journal of Applied Physics 23, 68 (1952).
  • Waldron [1960] R. Waldron, Perturbation theory of resonant cavities, Proceedings of the IEE - Part C: Monographs 107, 272 (1960).
  • Hahn et al. [2011] H. Hahn, W. Xu, P. Jain, and E. C. Johnson, Hom identification and bead pulling in the brookhaven erl, in Proceedings of the 2011 SRF Conference, Chicago, IL, USA, edited by S. Arora (Brookhaven National Laboratory, Upton, NY 11973, USA, 2011) p. THPO041, supported by Brookhaven Science Associates, LLC under Contract No. DE-AC02-98CH10886 with the U.S. DOE and DOE grant DE-SC0002496 to Stony Brook University.
  • Koubek et al. [2017] B. Koubek, A. Grudiev, and M. Timmins, rf measurements and tuning of the 750 mhz radio frequency quadrupole, Phys. Rev. Accel. Beams 20, 080102 (2017).
  • Berrutti et al. [2014] P. Berrutti, T. N. Khabiboulline, V. Poloubotko, G. Romanov, J. Steimel, V. Yakovlev, D. Li, and J. Staples, Pxie rfq bead pull measurements, in Proceedings of the 27th Linear Accelerator Conference (LINAC2014), edited by C. Carli (JACoW, Geneva, Switzerland, 2014) p. TUPP047.
  • Goudket et al. [2008] P. Goudket, C. Beard, P. McIntosh, G. Burt, A. Dexter, and R. M. Jones, Comparison of stretched-wire, bead pull and numerical impedance calculations on 3.9 ghz dipole cavities, in Proceedings of the European Particle Accelerator Conference (EPAC08), edited by P. Pierini (Genoa, Italy, 2008) p. MOPP128.
  • Argyropoulos et al. [2018] T. Argyropoulos, N. Catalan-Lasheras, A. Grudiev, G. Mcmonagle, E. Rodriguez-Castro, I. Syrachev, R. Wegner, B. Woolley, W. Wuensch, H. Zha, V. Dolgashev, G. Bowden, A. Haase, T. G. Lucas, M. Volpi, D. Esperante-Pereira, and R. Rajamäki, Design, fabrication, and high-gradient testing of an x𝑥xitalic_x-band, traveling-wave accelerating structure milled from copper halves, Phys. Rev. Accel. Beams 21, 061001 (2018).
  • Dassault Systèmes [2023] Dassault Systèmes, CST studio suite (2023).
  • Seibel [2017] A. Seibel, Untersuchungen zu einer neuen Alvarez-Struktur für den GSI Post-Stripper, doctoralthesis, Johann Wolfgang Goethe-Universität Frankfurt am Main (2017).
  • Du et al. [2017] X. Du, L. Groening, S. Mickat, A. Seibel, and O. K. Kester, Field stabilization of alvarez-type cavities, Phys. Rev. Accel. Beams 20, 032001 (2017).
  • Du and Groening [2018] X. Du and L. Groening, Compression and noise reduction of field maps, Phys. Rev. Accel. Beams 21, 084601 (2018).
  • Bethe and Schwinger [1943] H. A. Bethe and J. Schwinger, Perturbation theory for cavities, N.D.R.C. Report D1-117 (Cornell University, 1943).
  • Dolgashev et al. [2021] V. A. Dolgashev, L. Faillace, B. Spataro, S. Tantawi, and R. Bonifazi, High-gradient rf tests of welded x𝑥xitalic_x-band accelerating cavities, Phys. Rev. Accel. Beams 24, 081002 (2021).
  • Steele [1966] C. Steele, A nonresonant perturbation theory, IEEE Transactions on Microwave Theory and Techniques 14, 70 (1966).
  • Wegner et al. [2014] R. Wegner, W. Wuensch, G. Burt, and B. Woolley, Bead-pull measurement method and tuning of a prototype clic crab cavity, in Proceedings of LINAC2014, edited by C. Carli (JACoW, Geneva, Switzerland, 2014) p. MOPP035.
  • Heilmann et al. [2018] M. Heilmann, X. Du, L. Groening, M. Kaiser, S. Mickat, A. Seibel, and M. Vossberg, Scaled alvarez-cavity model investigations for the unilac upgrade, in 9th International Particle Accelerator Conference (IPAC2018), edited by S. Koscielniak (JACoW Publishing, Vancouver, BC, Canada, 2018) p. TUPAF079.
  • Podlech et al. [2007] H. Podlech, U. Ratzinger, H. Klein, C. Commenda, H. Liebermann, and A. Sauer, Superconducting ch structure, Phys. Rev. ST Accel. Beams 10, 080101 (2007).