Laser spectroscopy of the [Formula: see text], [Formula: see text] transitions in stored and cooled relativistic C[Formula: see text] ions.

The [Formula: see text] and [Formula: see text] transitions in Li-like carbon ions stored and cooled at a velocity of [Formula: see text] in the experimental storage ring (ESR) at the GSI Helmholtz Centre in Darmstadt have been investigated in a laser spectroscopy experiment. Resonance wavelengths were obtained using a new continuous-wave UV laser system and a novel extreme UV (XUV) detection system to detect forward emitted fluorescence photons. The results obtained for the two transitions are compared to existing experimental and theoretical data. A discrepancy found in an earlier laser spectroscopy measurement at the ESR with results from plasma spectroscopy and interferometry has been resolved and agreement between experiment and theory is confirmed.

Introduction

The possibilities for laser spectroscopy experiments at the upcoming FAIR facility require the development of experimental tools and techniques adapted to the relativistic ion energies offered by the new accelerators. In this regard a commissioning beam-time with Li-like carbon ions at the experimental storage ring (ESR) at the GSI Helmholtz Centre in Darmstadt provided the possibility to test new laser systems and a novel detection system for extreme UV (XUV) photons. The laser systems will be required for laser cooling of bunched relativistic ion beams, as it is being planned for the SIS100 synchrotron at FAIR[1,2]. For that purpose, continuous-wave (cw) and pulsed laser systems are developed by groups at the TU Darmstadt[3] and the HZDR/TU Dresden[4], respectively. For the detection of forward-emitted XUV fluorescence photons that are created in the de-excitation of highly charged ions at relativistic velocities, a detector system adapted to the large Lorentz boost and Doppler shift of these photons was developed in Münster[5]. The design is optimized for wavelengths around 10 nm as emitted in a proposed ESR experiment studying the decay of the state in Be-like krypton ions[6] but is sensitive over a wide wavelength region spanning from the UV to soft X-rays.

In the measurements described in this paper the and transitions in C-ions are investigated which are, besides their relevance for validation of theoretical atomic structure calculations (see e.g.[7] or[8]), of high importance for the astrophysical community where they are used for diagnostics of hot plasmas in the interstellar medium and stellar atmospheres[9]. The wavelengths of the two transitions were already determined in a previous laser-cooling experiment at the ESR[10], but were in conflict with earlier results from plasma spectroscopy[11] and from interferometric measurements[9]. In this new measurement special emphasis is therefore placed on a careful calibration of all experimental components and a thorough analysis of the systematic uncertainties.

Experimental setup

The measurements were performed with a beam of -ions that was overlapped with an anti-collinear laser beam in one of the straight sections of the ESR to excite either the or the transition (see Fig. 1). For this purpose two laser systems were available: a pulsed laser system from HZDR/TU-Dresden[4] and a continuous-wave laser system from TU-Darmstadt[3], both with a wavelength of . As the cw laser system provided a higher output power and a more narrow linewidth than the pulsed laser, it was used during most of the beamtime, including the measurements detailed in this paper. According to the equation for the relativistic Doppler shiftwith the ion velocity , the Lorentz-factor and an angle between laser and ion beam of , the laser photons are blue-shifted to wavelengths around in the rest frame of the ions. In order to obtain a well-determined absolute ion energy ( MeV/u) and a small relative momentum spread (), the ion beam was cooled by the ESR electron cooler[13]. Following the principle of electron cooling, the precise velocity of the ions is determined by the acceleration voltage of the electron beam inside the cooler. Scans across the two transitions of interest were realized by varying the electron cooler voltage, thus changing the velocity of the ions stored inside the ESR and thereby the Doppler-shifted laser wavelength in the rest frame of the ions. The dependence of the ion velocity on the effective electron cooler voltage is given bywhere e and are the elementary charge and electron mass, respectively. The acceleration voltage of the electron cooler was provided by a Heinzinger HNC 320.000-10 power supply and set on a terminal in the main control room of the storage ring. Monitoring of the actual voltage present at the cooler was achieved using a Julie Research JRL HVA-100 high-precision voltage divider connected to a 6.5 digit Keysight 34465A digital voltmeter (DVM).

Figure 1

Schematic view of the ESR with an enlarged view of the optical detection region. The ion beam is circling clockwise in the ring while the laser beam is counterpropagating with the ions inside the straight section of the ESR (figure adapted from[12]).

Fluorescence photons, which are emitted when the resonance condition for one of the two transitions is met, were detected by the new XUV detection system located in the region of overlap between the laser and the ion beam. Due to the Lorentz boost, the fluorescence photons are mainly emitted in the forward direction with the photon wavelengths being Doppler shifted down to (see inlay in Fig. 1).

CW Laser system

The measurements discussed in the following sections were performed using the continuous-wave (CW) laser system originally developed in the scope of a PhD thesis[3,14] and further optimized as reported in Refs.[15,16]. The setup is divided into five main components (see Fig. 2). A diode laser with an external resonator (external cavity diode laser, ECDL) serves as master oscillator operating at 1028 nm. The ECDL setup of Littrow design features a blazed grating that can be tilted via a piezoelectric actuator so that the output frequency of the laser system is tunable. The laser light from the ECDL is then amplified in an Ytterbium-doped fiber-amplifier system. Subsequently, the amplified light is frequency-doubled twice. The first frequency doubling is achieved by guiding the light in a single pass arrangement through a magnesium-doped periodically-poled lithium-niobate (MgO:PPLN) crystal. From there, the emerging green light of approximately 514 nm is guided into the UV cavity, where it is again frequency-doubled with a -barium borate (BBO) crystal. The resulting beam with  nm wavelength is then transported to the overlap region of laser and ion beams. In order to measure the output power during the beam-time, a fraction of the laser light was diverted by a beam splitter at the output of the UV cavity. As seen in Fig. 2, it was continously monitored by a Gentec-EO power detector model XLP12-1S-H2-D0 and the data was stored on a local laptop. The laser system provides 4 W of infrared power before the first frequency doubling that results in about 1 W of power at 514 nm. After the BBO crystal an average output power of 20 mW in the UV was obtained. To stabilize the laser frequency, a small fraction of the beam from the ECDL was coupled into a single-mode fiber and guided to a wavelength meter (WLM) from HighFinesse, type WS7-60. In a feedback loop deviations from the desired laser frequency measured by the device were then used to actuate the blazed grating, counteracting the fluctuation. In this way we obtained a frequency stability of the infrared laser light of about 5 MHz at 291 THz. The beam diameter directly after the laser system was 3 mm and about 10 mm in the overlap region between laser beam and ion beam.

Figure 2

Schematic view of the continuous-wave laser system.

XUV detection system

For an effective detection of forward emitted XUV photons in the ESR a novel detection system was developed based on the idea of detecting secondary electrons produced by the incoming photons from a suitable cathode material[17,18]. The detector consists of a movable cathode plate that can be driven into the beam line via a pressured air motor. In order not to interfere with the ion beam, the plate has a central 30 mm wide slit. XUV photons hitting the cathode will produce mostly low energetic (typically a few eV) secondary electrons that will be guided electromagnetically onto a microchannel plate (RoentDek DET40 MCP) in Chevron configuration, which is mounted inside the vacuum. A similar detection system for optical photons, making use of a movable parabolic copper mirror and a photomultiplier outside the vacuum[19], was already successfully applied in a measurement of the hyperfine structure splitting in H-like and Li-like bismuth[12,20].

Figure 3, left panel, shows a schematic view of the setup as it was implemented for simulations using the software package SIMION 8.1[21]. The cathode is coated with a 300 nm thick layer of caesium iodide (see Fig. 3, right panel), vacuum evaporated onto the cathode surface in the target laboratory at GSI, to enhance the secondary electron yield at the wavelengths of interest. The yield of photocathodes with a 300 nm coating of CsI was measured by Henke et al.[22] and was found to vary between 0.3 and 2.5 electrons per incident photon over a photon energy range from 100 eV to 10 keV. For the UV region data exist for CsI coated MCPs that provide quantum efficiencies between 20 and 30% for wavelengths between 20 and 150 nm[23]. Although the quantum efficiency drops sharply for longer wavelengths, we were able to successfully test the detector using UV light at a wavelength of 265 nm[17].

Figure 3

Left: schematic view of the XUV detection system with movable cathode (1), electrodes (2), magnet coils (3), magnetic field lines (4) and MCP (5). Right: cathode plate with 300 nm CsI coating.

A magnetic field, produced by solenoid coils mounted outside the vacuum, is used to guide the electrons to the MCP. The simulations resulted in an optimum configuration with one solenoid coil providing the main guiding field whereas a second solenoid coil closer to the cathode is operated with a significantly smaller current in opposite direction, widening the magnetic field lines to cover the complete cathode plate (see blue field lines in Fig. 3, left panel). Ring electrodes (labeled (2) in the figure) were meant to provide additional electrostatic guiding of the electrons but were found to have only little effect on the collection efficiency of the system and were grounded during the experiment. The collection efficiency, i.e. the percentage of electrons that are successfully guided from the cathode plate onto the active diameter of the MCP, was found to reach 80% with optimal settings of magnetic fields and cathode-plate voltage[18].

Figure 4 shows the interconnections of the detection system with the voltage supplies of the electrodes, the current supplies for the two solenoid coils and the front-end electronics for readout of the MCP. The high voltage for the channelplates is generated using an ISEG NHQ 224M HV supply. Charge pulses amplified in the MCP front and back plates and detected in the MCP anode are capacitively decoupled from the high voltage lines and amplified using a CAEN N979 fast amplifier. After converting the analog pulses to logical NIM signals using a CAEN N840 Constant Fraction Discriminator the latter are recorded using a GSI VUPROM (VME Universal Processing Module) module[24] serving both as a multi-hit TDC and a scaler module.

Figure 4

Interconnections of the XUV detection system with power supplies and front-end electronics of the data acquisition system.

Calibration and systematic uncertainties

A careful calibration of the different systems utilized in the setup must be carried out to minimize the systematic uncertainties in the deduced transition wavelengths. A short overview of the most important calibration procedures is given in the following.

HV divider scale factor

The uncertainty of the scale factor of the JRL HVA-100 high voltage divider used to monitor the electron cooler voltage is one of the main systematic uncertainties in the transition wavelength analysis presented in the “Resonance analysis” section. For the HVA-100 divider, we have three sets of calibration data obtained with different reference dividers in 2013[25], 2014[26] and 2017[27] at voltages similar to the  kV applied during the beam-time. Applying a linear fit to the scale factor as function of time, a scale factor for the measurement period in 2016 is extracted:The relatively large uncertainty of 16 ppm is mainly due to missing temperature monitoring of the divider during the 2016 beam-time (the divider exhibits a dependence of the scale factor on the environmental temperature of  ppm/K[27]).

DVM gain calibration

A calibration of the Keysight 34465A DVM was performed in advance of the beam-time to determine a possible offset and gain factor of the device. The offset was found to be 0.000000(20) V. The gain is determined utilizing a Fluke 732B reference source, which provides an output voltage of 10.000076(50) V, and was measured to bewhere is the voltage value read from the DVM.

Voltage dependence on cooler current

The cooler voltage, as measured by the DVM, exhibits a small dependence on the cooler current. This dependence was measured at the two central cooler voltages used in the beamtime over a range from 0 to 250 mA cooler current. Up to a cooler current of 122 mA we observe a linear drop in the measured DVM voltages by V/mA after which it saturates at a value of V. This correction is used when reconstructing missing DVM measurements at different cooler currents.

Set voltage calibration

Owing to problems with the readout software of the Keysight voltmeter, the measured electron cooler voltages were not written to storage for larger parts of the beam-time. This especially concerns data taken for the transition. For these parts of the data, we have to reconstruct the DVM measurements from the cooler voltage set values . To do so, we use periods of the beam-time where both DVM data and control room set voltages are available to develop a suitable conversion formula. As these data were collected at different electron cooler currents, we calculate the corresponding DVM values at zero cooler current by applying the correction explained in the previous paragraph.

Besides a small offset of the set voltage and a scale factor, the formula has to take into account the finite 18 bit resolution with which the set voltages were transfered to the Heinzinger power supply of the electron cooler. The maximum output voltage of the supply is 320 kV so the least significant bit of the digital set voltage corresponds to . With this the conversion formula was found to bewhere the function returns the integer part of its operand. The relative deviation of measured DVM voltages and reconstructed values stays below 10 ppm for all available datapoints.

Electron cooler workfunction

The work function W of a material has a direct impact on the electrostatic potential produced in the vacuum when a voltage U is applied to the material: . Electrons in the ESR electron cooler are generated in the dispenser cathode of an electron gun. It consists of a tungsten filament coated with barium, having a nominal work function of  eV[28]. The beampipe of the electron cooler, where the ion beam is superimposed with the electron beam, consists of stainless steel with a nominal work function of  eV[29]. Work functions are typically determined under ideal conditions with ultrapure surfaces. Work functions encountered under standard conditions can, therefore, deviate by up to  eV from the literature values[30]. The acceleration voltage felt by the electrons inside the cooler is determined by the difference of the electrode potentials of cathode and beampipe given bywhere the difference between the cathode voltage and the beampipe voltage is the voltage applied to the electron cooler using the Heinzinger power supply.

In conclusion, the potential of the electron cooler is shifted by  V, effectively lowering the electron energies. Considering an electron cooler voltage of approximately kV, the uncertainty of the work function difference on this voltage amounts to 21 ppm, representing one of the major contributions to the error budget.

Space charge corrections

Electron and ion space charges, present in the overlap region of both beams inside the electron cooler, modify the electric potential that accelerates the electron beam and therefore impact the velocity of the cooled ions. Whereas the effect of the electron space charges can be measured (see “Resonance analysis” section), the potential change due to the ion current has been estimated from the beam parameters and the drift tube geometry to be  V/mA[31]. To assess the effect of the ion space charge, the analysis has been performed with and without the correction resulting in a  ppm change in the extracted wavelength, which was then used as a conservative estimate of the related systematic uncertainty.

Laser wavelength

Two different wavelength meters were used to determine and stabilize the frequency of the CW-laser. (1) A HighFinesse WS6-200 wavelength meter, calibrated to a HeNe-laser, was used to measure the absolute frequency of the ECDL driving the system before and after the experiment. (2) A HighFinesse WS7-60 wavelength meter was utilized for a continuous stabilization of the laser frequency. Both wavelength meters operate in a range of 330 to 1180 nm and mainly differ in terms of accuracy. The WS6-200 has a 1 sigma uncertainty of 67 MHz and the WS7-60 has a 1 sigma uncertainty of 20 MHz. The measurement of the ECDL frequency with the calibrated WS6-200 amounted to  THz. The frequency is quadrupled in a two-step process and thus, the resulting laser wavelength is given by  nm.

Angle between laser and ion beam

An additional uncertainty arises from the unknown angle between laser and ion beam. The effective laser wavelength in the ions rest frame directly depends on this angle according to Eq. (1). The position of the ion beam in the vacuum beam-pipe of the ESR can be determined using two sets (horizontal and vertical) of position-calibrated “scrapers” inside the vacuum, which are a length  m apart. These scrapers are simple metal plates that can be moved with sub-millimetre precision, blocking part of the ion beam. By moving them from the outside towards the center of the beam-pipe, the radius and the central position of the beam can be determined. The same scrapers are also used to measure the position of the laser beam. The centre points of laser and ion beam—at the positions of the scrapers - were at maximum a distance  mm apart. The misalignment angle between laser and ion beam is therefore estimated to be causing a corresponding uncertainty in the transition wavelengths of 0.000027 nm for both transitions.

Resonance analysis

The following sections provide an overview of the experimental procedure during the beamtime and the extraction of the resonance wavelengths for the and transitions.

Experimental background

There are two major sources of experimental background for the fluorescence detection system. One is laser stray light creating photoelectrons from the cathode plate that are subsequently guided to the MCP. This leads to a constant background whose magnitude depends on the position of the cathode plate inside the beam line. The second source of background is caused by the ion beam ionizing or exciting residual gas molecules inside the beam pipe, which can lead to electrons or photo-electrons, respectively, finally hitting the MCP. The magnitude of this background component is proportional to the ion beam current. During initial tests of the detector in this beam-time, it was found that the detector could not be used with the cathode plate positioned around the ion beam in the center of the beam pipe due to excessive background rates that saturated the MCP. Therefore, the cathode was positioned at a distance of  cm from the ion beam. The resulting background contributions were then found to be  kcps/mA from the ion beam and  kcps/mW from the laser system[31].

Voltage scans

To determine the exact wavelengths of the and transitions, scans of the electron cooler voltage were performed thus changing the velocity of the stored ions inside the ESR and thereby the Doppler-shifted laser wavelength in the rest frame of the ions. Fluorescence photons that were created when the resonance condition for either of the two transitions was met, were detected and recorded by the data acquisition system. Figure 5 shows on the left panel the raw MCP rates (blue curve) recorded while scanning the cooler voltage (red curve) in 1 V steps three times over the resonance. The maximum rates, corresponding to the resonance condition being fulfilled, occur for a cooler voltage of approximately  V. The observed decrease in the peak amplitudes is caused by the loss of ions in the ring. The beam current was continously monitored using a DC current transformer in the beam diagnostic system. During the measurements we observed a mean lifetime of the beam current of about 62 s.

Figure 5

Left: raw detector rates (blue curve) for a scan of the electron cooler voltage (red curve) over the resonance region of the transition. Right: fluorescence rates of the same voltage scan after background subtraction, normalization to ion current and laser power and MCP gain correction.

To extract the exact resonance voltage, the experimental background is subtracted from the data which are subsequently normalized to the measured ion current and laser power . Finally a gain correction compensating for saturation effects of the MCP at high event rates is applied to obtain the fluorescence rate The gain correction can, following methods by Giudicotti[32] and Gershman[33], be written aswhere the MCP efficiency is stated by the manufacturer to be for electrons with energies between 100 and 1000 eV[34] (the electrons from the cathode are accelerated by a positive voltage of 300 V applied to the front-plate of the MCP). corresponds to the rate for which the MCP gain drops to 50% of the unsaturated value and has been determined to be kcps[31]. The effect of the gain correction on the reconstructed fluorescence rate reaches a factor 2 at an MCP rate of  kcps. The corrected fluorescence rates, obtained by applying Eq. (7) to the raw rates, are displayed in Fig. 5, right panel. The similar height obtained for all three peaks provides confidence in the normalizing procedure.

Accumulating all events from a given experimental run that belong to the same voltage step and plotting the corrected fluorescence rates against the effective electron cooler voltage (containing all corrections discussed in the “Calibration and systematic uncertainties” section), one obtains a resonance curve as displayed in Fig. 6, left panel. By fitting a Voigt profile to the data points, the precise voltage of the peak center is determined, together with its associated uncertainty.

Figure 6

Left: Observed resonance curve for the transition at a cooler current of 218 mA. Right: Measured resonance voltages as a function of the cooler current for the transition.

Electron space charge correction

To account for the electron space charge inside the ESR cooler, which modifies the acceleration potential felt by the electrons and therefore the ion velocity, the transition measurements were performed for a series of electron cooler currents, ranging from 10 to 250 mA. Figure 6 displays the observed resonance voltages for the transition (right panel) as a function of the electron cooler current. Performing linear fits to the two datasets allows to extrapolate the measurements to zero cooler current, thereby eliminating the influence of the electron space charge effect. The resulting resonance voltages for the two transitions are and .

Transition wavelengths

From the resonance voltages of the two transitions at zero cooler current given above, the transition wavelength in the ions rest frame can be calculated according to Eqs. (1) and (2), yielding and , where the statistical uncertainties result from the fit uncertainties of the resonance voltages and the systematic uncertainties result from the different contributions discussed in the “Calibration and systematic uncertainties” section and given in detail in Table 1 for the two transitions.

Table 1

Contributions of the individual systematic uncertainties to the overall error budget.

Parameter	Value	Uncertainty	Unit	Error contribution [nm]
				^2S_1/2-^2P_1/2	^2S_1/2-^2P_3/2
Work function	2.84	1.42	eV	0.00082	0.00081
Scale factor	9999.698	0.155		0.00060	0.00060
DVM gain	1.0000006	0.0000131		0.00050	0.00050
Set voltage cal.	var.	0.000067	V	0.00038	–
e-cooler current	var.	2	mA	0.00026	0.00026
DVM offset	0.0	0.00002	V	0.00012	0.00011
Ion space charge	var.	max.	V	0.00005	0.00004
Laser frequency	291.43340	0.00007	THz	0.00004	0.00004
Angle laser—ion beam	0.0	0.036	°	0.00003	0.00003
MCP gain	var.	max.		0.00001	0.00003
Root summed squares				0.00123	0.00116

A variable parameter is labeled ‘var.’. If the uncertainty of a parameter is labeled ‘max.’, the corresponding contribution to the error budget of the wavelength has been estimated by omitting the correction from the analysis.

Table 2 provides a comparison of the results of this work to previous measurements and theoretical predictions. The wavelengths determined for the two transitions are in good agreement with the experimental results obtained from interferometry[9] and plasma spectroscopy[11] and, within the uncertainties, also with the given theoretical predictions. It was possible to improve the precision with respect to the previous ESR measurement[10] by about a factor of three.

Table 2

Overview of measured and calculated wavelengths for the and transitions obtained in this work and available from the literature.

	Year	Method	λ(^2P_1/2) [nm]	λ(^2P_3/2) [nm]
Experiment
This work	2020	Laser spectroscopy	155.0779(12)_sys(1)_stat	154.8211(12)_sys(2)_stat
Schramm et al.[10]	2005	Laser spectroscopy	155.0705(39)_sys(3)_stat	154.8127(39)_sys(2)_stat
Griesmann et al.[9]	2000	Interferometry	155.0781(2)	154.8204(1)
Bockasten et al.[11]	1963	Plasma spectroscopy	155.0774(10)	154.8202(10)
Theory
Yerokhin et al.[8]	2017	RCIa	155.080(29)	154.827(36)
Borschevsky[35]	2014	FSCCb	155.075	154.820
Tupitsyn et al.[10]	2003	n/a	155.0739(26)	154.8173(53)
Johnson et al.[36]	1996	RMBPTc	155.078	154.819
Kim et al.[37]	1991	RMBPT + MCDHFd	155.060	154.804

aRelativistic configuration interaction.

bFock space coupled cluster.

cRelativistic many body perturbation theory.

dMultiConfiguration Dirac–Hartree–Fock.

Bold values indicate the 1 sigma uncertainties of the wavelengths.

A further improvement in accuracy could be achieved in future measurements by tackling the three most important systematic uncertainties: (1) Measuring the work function difference between cooler electron gun and drift electrodes, e.g. by performing laser spectrocopy on ions with a well-known transition. (2) Using a freshly calibrated high-voltage divider for the measurement of the cooler voltage and monitoring or stabilizing its temperature during the experiment. (3) Using a 8.5 digit DVM to measure the cooler voltage. While the achievable accuracy of a future work function measurement is difficult to estimate, the accuracy of the scale factor of the divider used in this measurement could be improved to the < 5 ppm level assuming a temperature stability of C, as was shown in the 2013 calibration of the device[25]. An 8.5 digit precision DVM together with a 10 V reference source (like e.g. the Fluke 732C) would allow for a 1 ppm precise readout of the divider output voltage.

On the experimental side, the new cw UV-laser system provided a stable and reliable beam to the experiment and the novel XUV detection system proved to have a high detection efficiency for Lorentz-boosted fluorescence photons. To lower the background sensitivity of the XUV detector, dedicated shielding electrodes and a stepper-motor driven motion of the cathode are being implemented.