1. Introduction
The measurement of the magnetic field in laboratory and astrophysical plasmas represents one of the most studied parameter in plasma physics, and at the same time is one of the most crucial parameters for steady-state plasma operation [
1]. A variety of different techniques is used to determine the magnetic field. For instance, for laboratory plasmas the magnetic field can be obtained e.g., by optical emission spectroscopy (OES) by measuring the Zeeman splitting of an emission line [
2,
3], by Thomson scattering of laser light [
4], by inductive sensors [
5], or by accurate description of the line shapes of quasi-isotropic fields [
6]. In interstellar plasmas, only OES can be used to determine magnetic fields due to the distance between plasma source and observer [
7,
8,
9]. Especially the diagnostics that do not disturb the plasma are attractive. TDLAS is a widely and commonly used technique in physics to measure the density and temperature of atomic or molecular species in the plasma [
10,
11,
12,
13]. Zeeman splitting of TDLAS spectra was already reported before [
14]; however, we are not aware of measurements of the magnetic field using a TDLAS system on metastable levels of noble gases.
The knowledge of the magnetic field is also important to measure the behavior of metastable levels of argon (Ar) in the linear plasma of PSI-2 [
15,
16]. Metastable Ar has a potential role in the excitation channel [
17] of the emission of reflected hydrogen atoms at the plasma-surface interface. The current application and further development of the Doppler-shifted reflectance measurements (DSRM) diagnostic [
18] in plasma cleaning discharges with electron temperature on the order of a few eV or less depends crucially on the answer to the following question: Are the metastable levels of Ar, in comparison with Kr or other noble gases, responsible for the strong emission of Balmer lines [
19]? In this work, we demonstrate the first measurement of the magnetic field using tunable diode laser absorption spectroscopy (TDLAS) in the linear magnetized plasma of the PSI-2 device.
2. Materials and Methods
For the measurement of the transition, a TDLAS system is installed at the PSI-2, as shown schematically in
Figure 1. The PSI-2 device has a hollow plasma profile with a cylindrical geometry that is a result of the source geometry. As an example of the plasma profile, the electron density for two magnetic field configurations is shown in
Figure 2a. The laser beam path is positioned in a way that the absorption is measured through the center of the plasma column. The measurements were performed at a normal gas pressure between 0.01–0.02 Pa. The PSI-2 plasma is confined by six magnetic coils that produce a magnetic field in
z-direction. A more detailed description of the PSI-2 device is given in [
16]. The standard currents for the six coils are given in
Table 1 (values are taken from [
20]). In the following, the standard current is shortened with
. The measurement is performed approx. 1.6 m away from the plasma source. At this position, the standard magnetic field can be calculated by the Biot–Savart law (Equation (
1)) using the values given in
Table 1:
[
20].
where
is the current density of coil
i and
is the magnetic constant. The magnetic field along the
z-axis is shown in
Figure 2b. Two different magnetic field configurations were analyzed, namely
and
, by varying the current in the coils. The Biot–Savart law describes the dependence of coil current and the magnetic field as linear. Thus the expected magnetic field should be
for
and
for
.
The laser used in this experiment is a DLpro cw laser from Toptica (TOPTICA Photonics AG, 82166 Gräfelfing) with a nominal output power of
at
and a maximum current of 185 mA. The laser wavelength can be adjusted freely between 750 and 795
, and a line width equals
. During the measurement, the wavelength is scanned in a range of about
with the frequency of
around a central wavelength by changing the grating position inside the laser. The grating is controlled by a piezo crystal. The laser beam is coupled into a single mode fiber and split up with a ratio of 99:1. One percent of the beam intensity is guided into a wavemeter (WS-6, produced by HighFinesse) for in situ wavelength monitoring. The precision of the wavemeter is
. The wavelength resolution of the TDLAS system is limited by the wavemeter. The other
of the power is guided into the PSI-2 device using an additional
single mode optical fiber. The laser output power is measured with a power meter after the beam splitter and is determined to 10–20
by operating the laser at a current of
. The major power losses are due to the coupling between laser and fiber. Before the laser beam is guided into the PSI-2 device, a filter wheel with different gray filters is set into the optical path in order to control the laser power that is guided into the plasma. After the interaction with the plasma, a bandpass filter and a chopper wheel are positioned in the optical path. A polarization plate is also introduced between the vacuum window and the photo diode. The circular polarized
component (
) is observed using a linear polarizer with the axis being perpendicular to the z-axis of the PSI-2. The linear polarized
component (
) is observed using the same linear polarizer with the axis being parallel to the z-axis of the PSI-2. The bandpass filter (Thorlabs FB760-10, CWL =
, full-width-half-maximum (FWHM) =
) is used to block the plasma radiation at all wavelengths except the laser wavelength of
. The chopper wheel chops the signal with a frequency of
. The chopped beam intensity is measured with a photo diode (Thorlabs, PDA36A-EC-Si Switchable Gain Detector, 350–1100
) right behind the chopper wheel. Then, it is converted into an electrical signal and analyzed with a dual lock-in amplifier (OrtoLock 9502). The linearity of the photo diode output signal was measured in the laboratory using the laser at a wavelength of
and a set of ND filters (Thorlabs FW212CNEB). The results of the measurements for different amplifications of the photo diode are shown in
Figure 3. The linearity of the photo diode signal was demonstrated by changing the laser power from
to
. Only at low output signals below
, the photo diode demonstrated a higher signal, as expected. The results were confirmed by using three different sets of amplifications (40, 50, and
). For the
gain, the signal was saturated for the filters with transmission above
(1–4). For the reference signal of the lock-in amplifier, the frequency signal of the chopper wheel was used. The output of the lock-in amplifier as well as the voltage signal of the piezo crystal were measured with an oscilloscope (
Figure 4). Together with the measured signal of the wavemeter, the x-axis of the oscilloscope can be recalculated in units of wavelength for further computation of the absorption signal. Note that due to the integration time of the lock-in amplifier, a time offset
occurred between the absorption signal and the piezo voltage. The time offset
(
) was corrected during the analysis of the spectrum. The ratio between the scanning period (5 s) and the integration time was 15, which should be sufficient to avoid large distortions in the obtained absorption line profile. For each spectrum, four different measurements were performed to eliminate the errors due to e.g., transmission losses in the optical path, as explained in detail in [
11,
12]:
—Plasma on, Laser on
—Plasma off, Laser on
—Plasma on, Laser off
—Plasma off, Laser off
The total transmittance of the signal can then be calculated (Equation (
2)):
The experimental spectra shown in this work were calculated using only the measurements of
and
. Furthermore, the spectra were corrected with regard to the linear background. The latter was induced by scanning the wavelength of the laser (see
Figure 4).
To calculate the magnetic field out of the absorption spectrum, the theoretical Zeeman spectrum of the line needs to be known. The laser wavelength is tailored for the transition from
(the transition is given in Paschen notation). The relative intensities of the Zeeman lines
(weak field) can be calculated using the Wigner 3-j symbols, assuming that all sublevels are equally populated [
21]:
where
is the momentum quantum number of the upper level and
of the lower level, and
and
are the corresponding magnetic quantum numbers. Possible transitions and the corresponding theoretical Zeeman pattern are shown in
Figure 5. To fit the experimental spectrum, eight (
component) or five (
component) Gaussian functions are used, each for one Zeeman component (note that the intensity of the Zeeman component of the transition
is 0):
where
is the wavelength,
is the amplitude,
is the center position, and
is the width of the Zeeman component
i. All Zeeman components
i are characterized by the same Doppler broadening, e.g., the temperature T, as a consequence of the constraint
, is defined. Only the left component of
Figure 5 (left) is fitted freely; for the other components, further constraints are defined. The amplitude of the Gauss functions is constrained by the ratio of the relative intensity given by Equation (
3):
, where
is the theoretical relative intensity of the Zeeman pattern. The constraint for the central position follows from the energy-shift of the Zeeman splitting
, where
and
are the energy or the wavelength of the field-free transition. The energy shift can be calculated as
, where
is the Bohr magneton and
g is the Landé factor. The Landé factors are taken from the NIST database (for the upper level,
and for the lower level,
) [
22]. As a consequence, to fit the experimental data, only the amplitude
, the magnetic field
B, and the temperature
T are used as free parameters. For the temperature, the formula for Doppler broadening of spectral lines
is used [
23] in SI units, where
is the full-width-half-maximum (FWHM) of the Gaussian function,
is the central wavelength of the line,
k is the Boltzmann constant,
T is the temperature of the light-absorbing species,
M is the mass of the species, and
c is the speed of light in vacuum. The FWHM of the Gaussian function can be calculated from
by
. Note that for the analysis of the absorption spectra, the hyperfine structure of Ar is neglected [
24].
3. Results
The following results were obtained by fitting the ratio
, as described in
Section 2, either for the
or the
component. We note that the obtained absorption spectra and results are line-integrated, and the magnetic field varies as a function of the position on the
z-axis in the PSI-2. In
Figure 6a, the obtained ratio
at a magnetic field configuration of
is plotted. The measurement was performed without a polarization filter, thus the Zeeman splitting was not visible. However, when a polarization filter in front of the photo diode (
Figure 1) was used, the observed signal split up into a
(
Figure 6b) and a
component (
Figure 6c). The
spectrum was obtained by adjusting the polarizer in a perpendicular direction to the
z-axis of the PSI-2 and showed two absorption peaks, whereas the
spectrum where the polarization direction was adjusted parallel to the
z-axis of the PSI-2 showed only one peak of absorption. This figure demonstrates the strong impact of the Zeeman effect on the line shape of the Ar I line. Neglecting the Zeeman effect in the PSI-2 for noble gases such as Ne, Ar, or Kr results in an overestimation of the derived gas temperature if the Zeeman effect is comparable or even larger than the Doppler broadening. In comparison to that for H or D lines, the Zeeman effect produced at most
of the line width. Thus, for H or D lines in the PSI-2 device, the influence of the Zeeman effect on the determination of the temperature is negligible [
19].
The
spectrum was fitted as explained in
Section 2, where the magnetic field
B was used as one fit parameter. The splitting for the
component consisted of eight instead of four components (
component), and the splitting was larger, which reduced the errors of the fit considerably. The ratio
is depicted in black, the green curve shows the fitted spectrum, and the red or blue dashed lines are single
or
components. To fit the
component, the obtained magnetic field of the
component was used and held constant. The temperature of both fits was compared to validate the fit of the
component. Both components can describe the observed absorption spectrum well. The obtained magnetic field for the
setup was
, which is a deviation of
from the expected value of
. For the case of
, the obtained magnetic field was
, which deviates by
from the expected value of
. The neutral gas temperature for
was determined to be
(
component) and
(
component), and for
,
(
component) and
(
component). The errors for magnetic field
B and temperature
T were taken from the fitting routine.
4. Discussion
It is the first time that the TDLAS diagnostic is used to measure the magnetic field in a low density plasma. As described in
Section 2, the transition
of metastable Ar is used for the experiments. The fits are in good agreement with the experimental data. It is observed that an increase of the magnetic field by 50% leads to an increase in the splitting. The small asymmetry in the absorption spectrum (e.g.,
Figure 6) cannot be described by the fit because a symmetrical line profile is demanded. The asymmetry could be partially caused by the non-linear background of the obtained data but also the additional contribution of the charge-exchange between Ar ions and neutrals. The experimentally determined magnetic fields are systematically lower, and the discrepancy between the expected values in vacuum and those obtained experimentally is around
. Considering that a number of effects for the calculation of the theoretical value of the magnetic field are not considered, the agreement is rather good. On the one hand, the most probable reason for lower measured magnetic field values in comparison to the vacuum calculation is the collective motion of ions in the angular direction, which produces a magnetic field in the opposite direction compared to the external applied magnetic field. On the other hand, the obtained values of the magnetic field are derived from the line-of-sight integrated signal. The experimental values have to be analyzed for different lines-of-sight to derive the dependence on the observation angle. Eventually, using Abel inversion to obtain spatially resolved profiles of the magnetic field could lead to a better agreement. This will be investigated in the near future. In addition to the magnetic field, the temperature of the metastable Ar atoms was also determined. For the magnetic field of 35–53
, the gas temperature was found to be around 400–800
. The comparison between the obtained temperature of the
and
components shows that the temperatures are in the same order of magnitude, which validates both the magnetic field and the temperature determination. The two-fold increase in the gas temperature with magnetic field strength is probably related to an increase of the plasma density by a factor of five, as obtained by Langmuir probe measurements (
Figure 2a). The TDLAS measurements provide data not only for magnetic field measurements, but also for gas rotation. Therefore, the integration time of the lock-in amplifier must be taken into account by analyzing the plasma rotation. For instance, the integration time of
for the lock-in amplifier at the scanning time of
(
Figure 4) results in a measured velocity of the atoms on the order of
. The wavelength calibration using an Ar lamp could be the most reliable solution to define the optimal integration time of the lock-in amplifier. Furthermore, the optical thickness of the absorption measurements was not analyzed. However, due to the weakening of the laser beam and the relatively low values of absorption, this effect is not expected to be significant under our conditions. The question about the optical thickness could be partially answered by also monitoring the intensity of other Ar I lines, such as
, originating from the upper level of the absorption transition. Moreover, the experimental data will be extended by performing a radial scan of the plasma and compare the emission and absorption spectra of Ar at different experimental conditions. Finally, the TDLAS diagnostic could help to identify the density of metastable levels of Ar and their role in the plasma–solid interaction.