Yan Wu is currently pursuing his master’s degree in Department of Modern Physics at University of Science and Technology of China. His primary research is related with experimental high energy physics at BESIII
Qinsong Zhou is currently pursuing his Ph.D. degree in School of Physical Science and Technology at Lanzhou University. His primary research is related with theoretical high energy physics. He has published several papers at Physical Review C, Physical Review D, and The European Physical Journal C
In this work, a combined analysis is performed on the processes of e+e−→ωπ0π0, e+e−→ωπ+π−, and e+e−→ωη to study possible ω excited states at approximately 2.2 GeV. The resonance parameters are extracted by simultaneous fits of the Born cross section line shapes of these processes. In the fit with one resonance, the mass and width are fitted to be (2207±14) MeV/c2 and (104±16) MeV, respectively. The result is consistent with previous measurements. In the fit with two resonances, the mass and width for the first resonance are fitted to be (2160±36) MeV/c2 (solution I), (2154±12) MeV/c2 (solution II) and (141±74) MeV (solution I), (152±77) MeV (solution II), respectively. The mass and width for the second resonance are fitted to be (2298±19) MeV/c2 (solution I), (2309±6) MeV/c2 (solution II) and (106±77) MeV (solution I), (99±23) MeV (solution II), respectively. The result is consistent with the theoretical prediction of ω(4S) and ω(3D). The intermediate subprocesses in e+e−→ωπ+π− are analyzed using the resonance parameters of the previous fits in this work. In the fit with one resonance, the fitting result of Γe+e−RBR is partially consistent with the previous result. In the fit with two resonances, the fitting result of Γe+e−RBR is of the same order of magnitude as the theoretical prediction. This work may provide useful information for studying the light flavor vector meson family.
Graphical Abstract
Simultaneous fit of e+e−→ωπ0π0, ωπ+π−, and ωη with one resonance.
Abstract
In this work, a combined analysis is performed on the processes of e+e−→ωπ0π0, e+e−→ωπ+π−, and e+e−→ωη to study possible ω excited states at approximately 2.2 GeV. The resonance parameters are extracted by simultaneous fits of the Born cross section line shapes of these processes. In the fit with one resonance, the mass and width are fitted to be (2207±14) MeV/c2 and (104±16) MeV, respectively. The result is consistent with previous measurements. In the fit with two resonances, the mass and width for the first resonance are fitted to be (2160±36) MeV/c2 (solution I), (2154±12) MeV/c2 (solution II) and (141±74) MeV (solution I), (152±77) MeV (solution II), respectively. The mass and width for the second resonance are fitted to be (2298±19) MeV/c2 (solution I), (2309±6) MeV/c2 (solution II) and (106±77) MeV (solution I), (99±23) MeV (solution II), respectively. The result is consistent with the theoretical prediction of ω(4S) and ω(3D). The intermediate subprocesses in e+e−→ωπ+π− are analyzed using the resonance parameters of the previous fits in this work. In the fit with one resonance, the fitting result of Γe+e−RBR is partially consistent with the previous result. In the fit with two resonances, the fitting result of Γe+e−RBR is of the same order of magnitude as the theoretical prediction. This work may provide useful information for studying the light flavor vector meson family.
Public Summary
We perform a combined analysis on e+e−→ωπ0π0, e+e−→ωπ+π−, and e+e−→ωη to study possible ω excited states around 2.2 GeV.
In the fit with one resonance, the mass and width are consistent with previous measurements.
In the fit with two resonances, the masses and widths are consistent with theoretical predictions.
The fitting result of Γe+e−RBR is partially consistent with experimental result and theoretical prediction.
Synchrotron radiation light sources have become important devices for scientific research. As a main part of the accelerator radio frequency (RF) system, which provides energy to electrons, superconducting cavities have been widely used in the last 50 years because of their low ohmic power loss[1, 2]. The SOLEIL-type cavity was named because it was first used in the 3rd generation light source SOLEIL[3]. Two 352.2 MHz SOLEIL-type superconducting modulers were installed on its storage ring in 2006 and 2008. They provided an accelerating voltage of 3 MV (750 kV/cell) and satisfied the power requirement for a beam current of 500 mA at a beam energy of 2.75 GeV[4]. Currently, SOLEIL cavities have the potential to operate under 500 mA beam current using a single crymodule with some improvements in their solid-state power amplifiers (SSPAs)[5].
The Hefei Advanced Light Facility (HALF), a large-scale scientific apparatus in the Fourteenth Five-Year Plan of China, is a soft X-ray and VUV light source based on a diffraction-limited storage ring (DLSR), which has a beam energy of 2.2 GeV and an emittance goal of less than 100 pm·rad. It is now under construction at the National Synchrotron Radiation Laboratory (NSRL). Based on KEK-B[6] and BEPC-II[7], a KEK-type superconducting cavity has been developed and will be used in HALF to meet the 350 mA beam current requirement[8]. In the future, with the increasement of beam current and number of insert parts, exploring SOLEIL-type superconducting cavities, which have two fundamental power couplers, may be meaningful. The linear section length of a 499.8 MHz SOLEIL module with two fundamental power couplers may be shorter than that of the two KEK modules because a higher power level could be easily realized.
All modes except the fundamental mode that exist in the cavity are called higher-order modes (HOMs). HOMs are excited when relativistic electrons passing through the cavity and must be extracted efficiently to prevent coupled bunch instability (CBI). HOM damping concepts may vary but are principally based on coaxial and waveguide couplers as well as beam line absorbers or any combination[9]. SOLEIL-type superconducting cavities were equipped with coaxial HOM couplers on their center pipe to reduce the linear section length and to have a moderate opening at the end beam pipes. Four HOM couplers were mounted on a 352 MHz SOLEIL cavity to extract the HOMs. Six HOM couplers were used to meet the HOM damping requirement of a 1500 MHz SUPER3HC 3rd harmonic cavity[10].
In this paper, the finite element analysis software CST Studio Suite was used to complete the research. First, the HOM properties of a bare 499.8 MHz SOLEIL-type superconducting cavity were analyzed. Second, coaxial HOM couplers were assumed to act as hooks to optimize their quantity and location, which has a great influence on the multipacting analysis. Then, the RF transmission properties and multipacting properties of the HOM couplers were optimized. Fourth, detailed HOM couplers were added to the model for the confirmation of the design and for overall optimization. Finally, a primary thermal simulation was presented to explore the cooling requirements of HOM couplers.
2.
Materials and methods
2.1
Cavity analysis
A 499.8 MHz SOLEIL-type superconducting cavity was analyzed by CST Suite Studio based on the 352 MHz SOLEIL cavity. The main parameters are listed in Table 1.
Table
1.
Main parameters for the 499.8 MHz SOLEIL-type superconducting cavity.
As shown in Fig. 1, the SOLEIL-type cavity have cells linked with a large beam pipe in between, called the center beam pipe (CBP), and terminated on smaller side beam pipes (SBPs). The radius of the CBP was specially designed to ensure that all the HOMs could be extracted into the CBP while the fundamental mode was still localized in cavity cells. When HOM couplers were assembled on CBP, this structure provided strong coupling to HOMs and very weak coupling to the fundamental mode.
Figure
1.
A 499.8 MHz SOLEIL-type superconducting cavity.
Theoretically, HOMs between the fundamental mode frequency f0 and the cutoff frequency of the SBP fc,SBP should be extracted by coaxial HOM couplers on the CBP. HOMs with frequencies above fc,SBP were assumed to be absorbed at the beam pipe absorbers or lossy metal tapers assembled outside the crymodule. In our case, as the radius of the SBP was set to 91.5 mm, the cutoff frequency was calculated to be 965.95 MHz in TE11 and 1261.87 MHz in TM01. Therefore, HOMs with frequencies between 499.8 MHz and 1261.87 MHz are of particular concern.
To find HOMs that may have a strong effect on the beam, the wake impedance is an accurate and convenient approach[11, 12]. It can be calculated by 2D software such as ABCI or 3D software such as the wakefield solver in CST Particle Studio. Using the cavity model in Fig. 1, two SBP ports were set as waveguide ports with five modes. The beam was set at a 1 mm bias from the cavity center, parallel to the red line in Fig. 1, and had a 34 mm bunch length. The wakelength was set to 800 m. The wake impedance was calculated and is shown in Fig. 2. The peak at approximately 499.8 MHz represented the fundamental mode, and the other peaks represented HOMs that may be dangerous. For HOMs that have a higher frequency, the electromagnetic field modes in the cavity cell and CBP may be different. They were named EM (frequency in MHz). For example, mode EM(716) has a field similar to that of TM111 in cavity cells, while the distribution in CBP is TE112. All HOM impedances were notably high and will be further increased with increasing simulation wakelength. This phenomenon can be attributed to the absence of HOM couplers in the model depicted in Fig. 1, which allows the HOMs to sustain resonance within the cavity.
Figure
2.
Wake impedance calculation results for a bare 499.8 MHz SOLEIL-type superconducting cavity: (a) longitudinal and (b) transverse.
According to the coupling method, coaxial-type HOM couplers can be divided into two categories, probe couplers or loop couplers, using electrical coupling or magnetic coupling, respectively[13]. Since the loop coupler could be designed to be dismountable and sometimes easier to use to form a filter structure, it was chosen as the basis for our HOM couplers.
2.2.1
Impedance threshold
The wake impedance threshold to cause CBI should not be exceed both longitudinal and transverse. It can be analytically calculated with the widely used Eqs. (1) and (2)[3, 9].
fHOMZth∥=2QsEατ∥I0,
(1)
f0Zth⊥=2Eβ⊥τ⊥I0,
(2)
where fHOM and f0 are the HOM resonance frequency and revolution frequency, E and I0 are the beam energy and average current, Qs is the synchrotron tune, α is the momentum compaction factor, β⊥ is the β function at the cavity location in the x or y direction, and τ∥ and τ⊥ are the longitudinal and transverse radiation damping times, respectively. These formulas were obtained by equating the radiation damping time with the respective instability rise time.
2.2.2
Coupler quantity and location
To achieve strong coupling, the plane of the hook should be as perpendicular as possible to the magnetic field line and close to the strong field. As shown in Fig. 3, the magnetic field of transverse HOMs such as TE112 tends to localize around the cavity, and longitudinal HOMs such as TM015 have a strong field around the center of the CBP. In our prediction, two kinds of couplers were needed: T-type couplers with hooks perpendicular to the beam axis and close to the cells, as shown in Fig. 3a, and L-type couplers parallel to the beam axis and mounted around the center of the CBP, as shown in Fig. 3b, which handle transverse and longitudinal HOMs, respectively.
Figure
3.
HOMs magnetic field distribution: (a) TE112 distribution, (b) TM015 distribution.
To estimate the quantity of HOM couplers and their location, couplers were assumed to be hooks added to the cavity model. Otherwise, the large mesh number and unacceptable calculation time will make the optimization impossible. A model with 4 T-type and 2 L-type couplers is shown in Fig. 4; L-type couplers are located around the center, and T-type couplers are located near the cavity cell. Couplers with different orientations were considered to absorb transverse HOMs with different polarizations.
Figure
4.
Cavity model with the coupling hooks, (a) the hooks, (b) the port setting.
We used the CST wakefield solver to optimize the model. During the optimization, the number of couplers and the distance between the couplers and cavity were varied to obtain a smaller impedance. All coupler ports were set as waveguide ports with TEM mode, and two SBP ports were set as waveguide ports with five modes. The beam was set to a 1 mm offset from the axis, and a long simulation wakelength of 1000 m was used to ensure that all the HOMs were extracted. A beam length of 34 mm was chosen to calculate the impedance accurately. Finally, 4 T-type couplers were symmetrically localized 137.5 mm from the cells and rotated 90 degrees on one side. At least 2 L-type couplers should be located around the center, as shown in Fig. 4, 362.5 mm from the cells.
The wake impedance results for the optimized model are shown in Fig. 5. The red lines are the impedance thresholds calculated by Eqs. (1) and (2) with the HALF parameters. The longitudinal impedance was well suppressed, and the maximum transverse impedance remained approximately 20 kΩ/m below the threshold. In this phase, a negative effect on the transverse impedance caused by the asymmetry of the L-type couplers was discovered but was considered acceptable for reducing the number of couplers.
Figure
5.
Wake impedance of the optimized model: (a) longitudinal and (b) transverse.
The coaxial line does not have the same cutoff property as the wakeguide, so coaxial HOM couplers need a filter structure at the fundamental mode frequency to prevent coupling of the fundamental mode, which in our case is 499.8 MHz. Filter structures can be realized by folding coaxial lines or adding extra rods or plates[14, 15]. The optimized equivalent circuits and models of our L-type and T-type couplers are shown in Fig. 6.
Figure
6.
Equivalent circuits and the optimized model, (a) L-type coupler, (b) T-type coupler.
Both couplers have one filter structure: a capacitor plate around the support rod for L-type couplers and a lengthened hook facing the outer conductor for T-type couplers, as shown in the red boxes of Fig. 6a and b, respectively. Since T-type couplers are close to the cavity cell, a filter structure around the coupler hook may block the fundamental mode from entering the coupler, reducing extra heat loss. In the L-type coupler case, the fundamental mode fields were weak around the coupler, and the lengthened hook did not work as a filter after 90 degrees of rotation, so the hook was shortened, and a capacitor plate was added around the support rod.
The transport properties of these couplers were analyzed by beampipe-coupler models using CST microwave (MW) studio. As shown in the left part of Fig. 7, the TE11 or TM01 mode was launched at the side of the beam pipe, where it was set as a waveguide port named port1 or port3. Port2 was set to the TEM mode. The S21 results are shown in the right part of Fig. 7. Since the fundamental mode TM010 excites the TM01 mode at the beam pipe port, the TM01 mode S21 should receive special attention. After the optimization of the coupler geometry, both couplers exhibited good reflection properties—over −80 dB for L-type couplers and −75 dB for T-type couplers—at 499.8 MHz.
Figure
7.
Beampipe-coupler models and S21 results for the (a) L-type coupler and (b) T-type coupler.
Longitudinal HOMs were extracted from L-type couplers. As shown by the blue line in Fig. 2, the longitudinal HOMs were located between 0.80 GHz and 1.26 GHz, with the strongest occurring at approximately 1.00 GHz. The TM01 launched S21 for the L-type coupler in this range was mostly above –20 dB, which predicted good transmission. On the other hand, the transverse HOMs were located between 0.55 GHz and 1.10 GHz, and the T-type couplers also exhibited good transmission properties in two polarization directions.
2.4
Multipacting analysis of HOM couplers
Multipacting (MP) is a phenomenon caused by electrons that resonate with electromagnetic fields. When some primary electrons gain energy from the field and collide with conductors, secondary electrons may be generated. If secondary electrons act similar to primary electrons and form more secondary electrons, the number of electrons will increase exponentially, and then the MP will occur. MP results in extra heat loss, quench of the superconducting, or even breakdown of HOM couplers.
The particle in cell (PIC) solver in CST particle studio was used for MP analysis because it has several advantages, such as convenient particle settings and better computer memory processing compared with the tracking solver, while maintaining comparable accuracy[16]. The electromagnetic field was calculated by a model including the cavity and HOM couplers using the eigenmode solver of CST MW Studio and then input into the PIC solver. The model and electromagnetic field are presented in Fig. 8. Since the fundamental mode is below the cutoff frequency of the CBP, fields were localized inside cavity cells. The brown parts in Fig. 8 were alumina ceramics.
Figure
8.
Electromagnetic field model, (a) cavity model with couplers, (b) fundamental mode electric field.
When the cavity is in operation, the electromagnetic field in the HOM couplers could be as strong as the field in cavity cells. Thus, niobium was chosen as the material for the inner conductor of HOM couplers, enabling the formation of a superconductive state, as illustrated in the blue parts of Fig. 9a.
Figure
9.
MP analysis model and true SEY used: (a) model and (b) SEY curve.
In the PIC solver, the Furman emission model was used to simulate the MP effect. The backscattered electrons, rediffused electrons and true electrons were all considered in this model. A niobium shell must be added outside the cavity to make the cavity boundary reflectable, as shown in Fig. 9a. The secondary emission yield (SEY) of niobium and copper was obtained from the offset in CST particle studio, as shown in Fig. 9b.
2.4.1
T-type coupler multipacting analysis
Close to the cavity, the T-type couplers were analyzed first. All the inner conductor surfaces were set as particle sources, as shown in Fig. 10a. Over 2000 primary electron launch points were uniformly distributed on the surface. The electron launch angle was set to a random range from 0° to 89.9° with an initial energy of 2 eV. In the time domain, the launch current was set as a 2 ns rectangular pulse, covering the first period of the fundamental mode resonance. Each calculation was run for 15 ns, and the space-charge effect was neglected.
Figure
10.
MP analysis for T-type couplers: (a) particle source area, (b) particle monitoring where MPs occurred, and (c) <SEY> results.
The time-averaged secondary emission yield (<SEY>) was used to indicate whether MP occurred, and it could be calculated by Eq. (3):
<SEY>=<Ie><Ic>,
(3)
where <Ie> and <Ic> represent the time-averaged emission current and collision current, respectively, calculated by electrons emitted from or colliding with the entire conductor surface and averaged from the simulation time of 4 ns to 15 ns. When <SEY> is greater than 1, the emission current is greater than the collision current, which always represents a high possibility that MP may occur.
Initially, the materials of the upper coaxial line and the hook were set as copper and niobium, respectively, but MP occurred at an acceleration gradient of approximately 4 MV/m. Particle monitor revealed that the MP was located at the capacity part, as shown in the upper panel of Fig. 10b. After the upper part was replaced by niobium, the MP was suppressed well between the 2 MV/m and 16 MV/m acceleration gradients.
2.4.2
L-type coupler multipacting analysis
Contrary to what we predicted, L-type couplers have a weak field strength inside but are more likely to experience MP. Similar to T-type couplers, the entire inner conductor surface was set as a particle source, as shown in Fig. 11a. All primary electron settings and the <SEY> calculations were consistent with the T-type coupler analysis.
Figure
11.
MP analysis for L-type couplers: (a) particle source area, (b) particle monitoring where MPs occurred, and (c) <SEY> results.
After optimization, MP still occurred over 10 MV/m at the capacity plate facing the outer conductor, as shown in the upper panel of Fig. 11b. An inner conductor capacity plate might also lead to an MP under a higher acceleration gradient. As a storage ring cavity, the acceleration gradient was always set between 2 MV/m and 6 MV/m, rarely more than 8 MV/m. In this range, optimized L-type couplers had a low probability of having an MP effect. The HALF required a 1.5 MV RF voltage, and the acceleration gradient could be calculated to be 2.48 MV/m per cell using a 499.8 MHz SOLEIL-type cavity. It could be predicted that no multipacting will occur under an operating acceleration gradient in L-type or T-type HOM couplers.
3.
Results and discussion
3.1
Wake impedance recalculation
After RF and MP analysis, the structure of HOM couplers is finally determined. To ensure that these couplers could dampen HOMs efficiently, the wake impedance was recalculated with a model including detailed HOM couplers and the cavity, as shown in Fig. 12a.
Figure
12.
The cavity model and its impedances: (a) model with 4 T-type and 2 L-type couplers, (b) longitudinal impedance, and (c) transverse impedance.
In this model, the coupling hooks were replaced by 4 T-type couplers and 2 L-type couplers in detail. The beam length was set to 50 mm after a tradeoff between accuracy and calculation time, and the beam position was set to 1 mm biased from the beam axis. To ensure that all the HOMs with frequencies less than 1.3 GHz could be extracted, two SBP ports were set as waveguide ports with five modes, and all the HOM coupler ports were set as waveguide ports with three modes. The impedance results are presented in Fig. 12b and c, and the longitudinal impedances were well damped. Unfortunately, the transverse impedance greatly exceeds the threshold.
To ensure that the transverse impedance is below the threshold, the asymmetry effect of L-type couplers was considered. Two extra L-type couplers were added opposite to the former L-type couplers shown in Fig. 13a, all of which are located 362.5 mm from the cavity. During the calculations, all beam parameter settings were inherited from the previous model. The addition of two L-type couplers resulted in a significant reduction in both longitudinal and transverse impedance, as shown in Fig. 13b and c.
Figure
13.
The cavity model and its impedances: (a) model with 4 T-type and 4 L-type couplers, (b) longitudinal impedance, and (c) transverse impedance.
Finally, a preliminary thermal analysis of HOM couplers was performed. Using CST Microwave Studio and Mphsics Studio, a bidirectional electromagnetic-thermal calculation was studied. The model is presented in Fig. 14a. Since the cells of the SOLEIL-type cavity were completely immersed in 4.2 K liquid helium and the CBP part was in vacuum. The cavity cells were assumed well cooled and simplified with two 4.2 K thermal anchors on both sides of the CBP[3, 4]. T-type couplers were also assumed to be cooled by 4.2 K liquid helium because they were close to the cavity cell. Additionally, a 5 K thermal anchor was added at the top of each HOM coupler for thermal calculations. The 4.2 K and 5 K thermal anchors are shown in blue and red, respectively, in Fig. 14a. The thickness of the CBP niobium wall was 4 mm. The temperature-dependent heat conductivity and electrical conductivity of niobium were obtained from Fermilab documents[17]. The surface resistance of the superconductive niobium was set to 20 nΩ at 4.2 K and 1 GHz, with a residual resistance of 10 nΩ. This value exhibits variations with both temperature and frequency according to the following formula[2]:
Figure
14.
Thermal analysis: (a) boundary conditions, (b) temperature distribution under a 10 MV/m acceleration gradient, and (c) temperature distribution under a 3 MV/m acceleration gradient and 2 kW HOM power.
where A is a constant factor, T is the temperature, f is the frequency, Δ is the half-energy gap of niobium, k is the Boltzmann constant and R0 is the residual resistance. Other materials, such as ceramics, were obtained from the CST library as a default.
In this model, the TM01 mode waveguide ports on both sides of the CBP represented the power launched by the TM010 mode in the cavity cells, as shown in Fig. 14a. The initial temperature distribution was calculated without RF power and input into microwave studio. Thereafter, new temperature distribution was computed based on the RF losses derived from the initial distribution. This iterative process was repeated until stable results were achieved. Two simulation conditions are presented as follows.
Since MP may not appear until a 10 MV/m acceleration gradient is reached, the electromagnetic fields were scaled to 10 MV/m at 499.8 MHz for thermal analysis first. As shown in Fig. 14b, the hook parts of all HOM couplers remained superconducting. The maximum temperature was 8.3 K at the ceram part of the T-type couplers. With such a small temperature increase, there should be no thermal breakdown until MP occurs.
The second condition was established to simulate the temperature distribution under an acceleration gradient of 3 MV/m, taking HOM power into account, assuming the module was in operation. The HOM power could be calculated by the loss factor:
PHOM=(k||−kFM)Ib2Nfrev,
(5)
where k|| is the total loss factor and kFM is the loss factor for the fundamental mode. Ib is the beam current, N is the number of bunches, and frev is the revolution frequency. The loss factor for the fundamental mode kFM can be calculated by Eq. (6).
kFM=ωFM4(RQ)FMexp(−ω2FMσ2zc2),
(6)
where ωFM and (R/Q)FM is the angular frequency and the R over Q value of the fundamental mode, σz is the longitudinal beam length and c is the speed of light.
Using the model shown in Fig. 1, when the beam length was 5 mm, the total loss factor and the fundamental mode loss factor was computed to be 0.6 V/pC and 0.14 V/pC, respectively. With the parameters of HALF, the HOM power could be predicted to be less than 1 kW. However, the tapers and beam pipe HOM absorbers outside the model were proven to be important contributors to losses at short beam lengths[18]. Consequently, we assumed a total HOM power of 2 kW. On each side of the CBP, the 499.8 MHz and 1 GHz fields were scaled to fit a 3 MV/m gradient and 1 kW HOM power, respectively. The simulation results are shown in Fig. 14c. The maximum temperature was 5.2 K. The temperature distribution of the L-type couplers was nearly identical to that of under 10 MV/m gradient. While the T-type couplers had a lower temperature related to the decrease in the acceleration gradient. As a summary, the HOM power has minor influence on the temperature distribution. There are no risks of thermal breakdown under this condition.
4.
Conclusions
SOLEIL-type superconduting cavities have two fundamental power couplers, which may easily provide a higher power level. Two kinds of coaxial HOM couplers were designed for a 499.8 MHz SOLEIL-type superconducting cavity. The geometry, quantity, and location of these couplers were optimized. The RF transfer characteristics and MP properties were good until 10 MV/m was achieved. The HOM damping requirement of HALF can be greatly satisfied using 4 L-type and 4 T-type couplers. A preliminary thermal analysis showed that there should be no thermal breakdown until MPs occur. This work will provide preliminary research for a 499.8 MHz SOLEIL-type superconducting cavity.
Acknowledgements
This work was supported by the National Natural Science Foundation of China (12035013), Joint Large-Scale Scientific Facility Funds (U1732263), and National Key R&D Program of China ( 2020YFA0406403).
Conflict of interest
The authors declare that they have no conflict of interest.
We perform a combined analysis on e+e−→ωπ0π0, e+e−→ωπ+π−, and e+e−→ωη to study possible ω excited states around 2.2 GeV.
In the fit with one resonance, the mass and width are consistent with previous measurements.
In the fit with two resonances, the masses and widths are consistent with theoretical predictions.
The fitting result of Γe+e−RBR is partially consistent with experimental result and theoretical prediction.
Wang J Z, Chen D Y, Liu X, et al. Costructing J/ψ family with updated data of charmoniumlike Y states. Phys. Rev. D,2019, 99 (11): 114003. DOI: 10.1103/PhysRevD.99.114003
[2]
Wang J Z, Qian R Q, Liu X, et al. Are the Y states around 4.6 GeV from e+e– annihilation higher charmonia? Physical Review D,2020, 101: 034001. DOI: 10.1103/PhysRevD.101.034001
[3]
Wang J Z, Sun Z F, Liu X, et al. Higher bottomonium zoo. The European Physical Journal C,2018, 78 (11): 915. DOI: 10.1140/epjc/s10052-018-6372-1
[4]
Pang C Q, Wang J Z, Liu X, et al. A systematic study of mass spectra and strong decay of strange mesons. The European Physical Journal C,2017, 77 (12): 861. DOI: 10.1140/epjc/s10052-017-5434-0
[5]
Song Q T, Chen D Y, Liu X, et al. Charmed-strange mesons revisited: Mass spectra and strong decays. Physical Review D,2015, 91: 054031. DOI: 10.1103/PhysRevD.91.054031
[6]
Song Q T, Chen D Y, Liu X, et al. Higher radial and orbital excitations in the charmed meson family. Physical Review D,2015, 92: 074011. DOI: 10.1103/PhysRevD.92.074011
[7]
Wang J Z, Chen D Y, Song Q T, et al. Revealing the inner structure of the newly observed D2*(3000). Physical Review D,2016, 94: 094044. DOI: 10.1103/PhysRevD.94.094044
[8]
Particle Data Group, R. L. Workman R L, Burkert V D, et al. Review of particle physics. Progress of Theoretical and Experimental Physics,2022, 8: 083C01. DOI: 10.1093/ptep/ptac097
[9]
Pang C Q, Wang Y R, Hu J F, et al. Study of the ω meson family and newly observed ω-like state X(2240). Physical Review D,2020, 101: 074022. DOI: 10.1103/PhysRevD.101.074022
[10]
Barnes T, Close F E, Page P R, et al. Higher quarkonia. Physical Review D,1997, 55: 4157–4188. DOI: 10.1103/PhysRevD.55.4157
[11]
Ebert D, Faustov R N, Galkin V O. Masses of light mesons in the relativistic quark model. Modern Physics Letters A,2005, 20: 1887–1893. DOI: 10.1142/S021773230501813X
[12]
Ebert D, Faustov R N, Galkin V O. Mass spectra and Regge trajectories of light mesons in the relativistic quark model. Physical Review D,2009, 79: 114029. DOI: 10.1103/PhysRevD.79.114029
[13]
Wang X, Sun Z F, Chen D Y, et al. Nonstrange partner of strangeonium-like state Y(2175). Physical Review D,2012, 85: 074024. DOI: 10.1103/PhysRevD.85.074024
[14]
Anisovich A V, Baker C A, Batty C J, et al. I=0, C=−1 mesons from 1940 to 2410 MeV. Physics Letters B,2002, 542: 19–28. DOI: 10.1016/S0370-2693(02)02303-1
Omega Photon Collaboration, Atkinson M, Axon T J, et al. Photon diffractive dissociation to ρρπ and ρπππ states. Zeitschrift für Physik C Particles and Fields,1988, 38: 535–541. DOI: https://doi.org/10.1007/BF01624357
[17]
Aubert B, Bona M, Boutigny D, et al. The e +e– → 2(π +π –) π0, 2(π+π–) eta, K+K–π+π–π0 and K +K –π +π –η cross sections measured with initial-state radiation. Physical Review D,2007, 76: 092005. DOI: 10.1103/PhysRevD.76.092005
[18]
Lees J P, Poireau V, Tisserand V, et al. Study of the reactions e +e− → π +π −π0π0π0and π+π−π0π0η at center-of-mass energies from threshold to 4.35 GeV using initial-state radiation. Physical Review D,2018, 98: 112015. DOI: 10.1103/PhysRevD.98.112015
[19]
Lees J P, Poireau V, Tisserand V, et al. Resonances in e+e− annihilation near 2.2 GeV. Physical Review D,2020, 101: 012011. DOI: 10.1103/PhysRevD.101.012011
[20]
The BESIII collaboration, Ablikim M, Achasov M N, et al. Measurement of e +e −→ωπ +π − cross section at √s = 2.000 to 3.080 GeV. Journal of High Energy Physics,2023, 1: 111. DOI: https://doi.org/10.1007/JHEP01(2023)111
[21]
M. Ablikim, Achasov M N, Adlarson P, et al. Measurement of the e +e −→ ωπ0π0 cross section at center-of-mass energies from 2.0 to 3.08 GeV. Physical Review D,2022, 105: 032005. DOI: 10.1103/PhysRevD.105.032005
[22]
Ablikim M, Achasov M N, Adlarson P, et al. Observation of a resonant structure in e+e−→ωη and another in e+e− → ωπ0 at center-of-mass energies between 2.00 GeV and 3.08 GeV. Physics Letters B,2021, 813: 136059. DOI: 10.1016/j.physletb.2020.136059
[23]
Zhou Q S, Wang J Z, Liu X. Role of the ω(4S) and ω(3D) states in mediating the e +e −→ωη and ωπ0π0 processes. Physical Review D,2022, 106: 034010. DOI: 10.1103/PhysRevD.106.034010
[24]
Wang J Z, Wang L M, Liu X, et al. Deciphering the light vector meson contribution to the cross sections of e +e– annihilations into the open-strange channels through a combined analysis. Physical Review D,2021, 104: 054045. DOI: 10.1103/PhysRevD.104.054045
Figure
1.
A comparison of resonance parameters of these reported ω states with masses around 2.2 GeV[14–22].
Figure
2.
Solution I of the simultaneous fit of e+e−→ωπ0π0, ωπ+π−, and ωη with one resonance. Black dots with error bars are from experimental data, and the error bars include both statistical and systematic uncertainties. The blue solid curve represents the total fit. The red dashed curve represents the contribution of resonant part. The green dashed curve presents the contribution of non-resonant part. The pink curve represents the interference between the two parts.
Figure
3.
Solution II of the simultaneous fit of e+e−→ωπ0π0, ωπ+π−, and ωη with one resonance. The symbols’ descriptions are consistent with Fig. 2.
Figure
4.
Solution I of the simultaneous fit of e+e−→ωπ0π0, ωπ+π−, and ωη with two resonances. Black dots with error bars are from experimental data, and the error bars include both statistical and systematic uncertainties. The blue solid curve represents the total fit. The red dashed curve represents the contribution of ω(4S). The cyan dashed curve represents the contribution of ω(3D). The green dashed curve presents the contribution of non-resonant part. The pink curve represents the interference between the resonant and non-resonant parts.
Figure
5.
Solution II of the simultaneous fit of e+e−→ωπ0π0, ωπ+π− and ωη with two resonances. The symbols’ descriptions are consistent with Fig. 4.
Figure
6.
Solution I of the simultaneous fit of the intermediate modes in e+e−→ωπ+π− with one resonance. The symbol descriptions are consistent with Fig. 2.
Figure
7.
Solution II of the simultaneous fit of the intermediate modes in e+e−→ωπ+π− with one resonance. The symbol description is consistent with Fig. 2.
Figure
8.
Solution I of the simultaneous fit of the intermediate modes in e+e−→ωπ+π− with two resonances. The symbols’ descriptions are consistent with Fig. 4
Figure
9.
Solution II of the simultaneous fit of the intermediate modes in e+e−→ωπ+π− with two resonances. The symbols’ descriptions are consistent with Fig. 4.
References
[1]
Wang J Z, Chen D Y, Liu X, et al. Costructing J/ψ family with updated data of charmoniumlike Y states. Phys. Rev. D,2019, 99 (11): 114003. DOI: 10.1103/PhysRevD.99.114003
[2]
Wang J Z, Qian R Q, Liu X, et al. Are the Y states around 4.6 GeV from e+e– annihilation higher charmonia? Physical Review D,2020, 101: 034001. DOI: 10.1103/PhysRevD.101.034001
[3]
Wang J Z, Sun Z F, Liu X, et al. Higher bottomonium zoo. The European Physical Journal C,2018, 78 (11): 915. DOI: 10.1140/epjc/s10052-018-6372-1
[4]
Pang C Q, Wang J Z, Liu X, et al. A systematic study of mass spectra and strong decay of strange mesons. The European Physical Journal C,2017, 77 (12): 861. DOI: 10.1140/epjc/s10052-017-5434-0
[5]
Song Q T, Chen D Y, Liu X, et al. Charmed-strange mesons revisited: Mass spectra and strong decays. Physical Review D,2015, 91: 054031. DOI: 10.1103/PhysRevD.91.054031
[6]
Song Q T, Chen D Y, Liu X, et al. Higher radial and orbital excitations in the charmed meson family. Physical Review D,2015, 92: 074011. DOI: 10.1103/PhysRevD.92.074011
[7]
Wang J Z, Chen D Y, Song Q T, et al. Revealing the inner structure of the newly observed D2*(3000). Physical Review D,2016, 94: 094044. DOI: 10.1103/PhysRevD.94.094044
[8]
Particle Data Group, R. L. Workman R L, Burkert V D, et al. Review of particle physics. Progress of Theoretical and Experimental Physics,2022, 8: 083C01. DOI: 10.1093/ptep/ptac097
[9]
Pang C Q, Wang Y R, Hu J F, et al. Study of the ω meson family and newly observed ω-like state X(2240). Physical Review D,2020, 101: 074022. DOI: 10.1103/PhysRevD.101.074022
[10]
Barnes T, Close F E, Page P R, et al. Higher quarkonia. Physical Review D,1997, 55: 4157–4188. DOI: 10.1103/PhysRevD.55.4157
[11]
Ebert D, Faustov R N, Galkin V O. Masses of light mesons in the relativistic quark model. Modern Physics Letters A,2005, 20: 1887–1893. DOI: 10.1142/S021773230501813X
[12]
Ebert D, Faustov R N, Galkin V O. Mass spectra and Regge trajectories of light mesons in the relativistic quark model. Physical Review D,2009, 79: 114029. DOI: 10.1103/PhysRevD.79.114029
[13]
Wang X, Sun Z F, Chen D Y, et al. Nonstrange partner of strangeonium-like state Y(2175). Physical Review D,2012, 85: 074024. DOI: 10.1103/PhysRevD.85.074024
[14]
Anisovich A V, Baker C A, Batty C J, et al. I=0, C=−1 mesons from 1940 to 2410 MeV. Physics Letters B,2002, 542: 19–28. DOI: 10.1016/S0370-2693(02)02303-1
Omega Photon Collaboration, Atkinson M, Axon T J, et al. Photon diffractive dissociation to ρρπ and ρπππ states. Zeitschrift für Physik C Particles and Fields,1988, 38: 535–541. DOI: https://doi.org/10.1007/BF01624357
[17]
Aubert B, Bona M, Boutigny D, et al. The e +e– → 2(π +π –) π0, 2(π+π–) eta, K+K–π+π–π0 and K +K –π +π –η cross sections measured with initial-state radiation. Physical Review D,2007, 76: 092005. DOI: 10.1103/PhysRevD.76.092005
[18]
Lees J P, Poireau V, Tisserand V, et al. Study of the reactions e +e− → π +π −π0π0π0and π+π−π0π0η at center-of-mass energies from threshold to 4.35 GeV using initial-state radiation. Physical Review D,2018, 98: 112015. DOI: 10.1103/PhysRevD.98.112015
[19]
Lees J P, Poireau V, Tisserand V, et al. Resonances in e+e− annihilation near 2.2 GeV. Physical Review D,2020, 101: 012011. DOI: 10.1103/PhysRevD.101.012011
[20]
The BESIII collaboration, Ablikim M, Achasov M N, et al. Measurement of e +e −→ωπ +π − cross section at √s = 2.000 to 3.080 GeV. Journal of High Energy Physics,2023, 1: 111. DOI: https://doi.org/10.1007/JHEP01(2023)111
[21]
M. Ablikim, Achasov M N, Adlarson P, et al. Measurement of the e +e −→ ωπ0π0 cross section at center-of-mass energies from 2.0 to 3.08 GeV. Physical Review D,2022, 105: 032005. DOI: 10.1103/PhysRevD.105.032005
[22]
Ablikim M, Achasov M N, Adlarson P, et al. Observation of a resonant structure in e+e−→ωη and another in e+e− → ωπ0 at center-of-mass energies between 2.00 GeV and 3.08 GeV. Physics Letters B,2021, 813: 136059. DOI: 10.1016/j.physletb.2020.136059
[23]
Zhou Q S, Wang J Z, Liu X. Role of the ω(4S) and ω(3D) states in mediating the e +e −→ωη and ωπ0π0 processes. Physical Review D,2022, 106: 034010. DOI: 10.1103/PhysRevD.106.034010
[24]
Wang J Z, Wang L M, Liu X, et al. Deciphering the light vector meson contribution to the cross sections of e +e– annihilations into the open-strange channels through a combined analysis. Physical Review D,2021, 104: 054045. DOI: 10.1103/PhysRevD.104.054045