MD Isabelle Mosnier
email: isabelle.mosnier@aphp.fr

MD Olivier Sterkers

PhD Yann Nguyen

PhD Ghizlene Lahlou

Benefits in noise from sound processor upgrade in thirty-three cochlear implant users for more than 20 years

Keywords: long-term benefit, hearing in noise, upgrading, speech processor

come

Benefits in noise from sound processor upgrade in thirty-three cochlear implant users for more than 20 years

Introduction

Deafness has been established to play a key role on cognitive functions all along the life, as early deafness restoration permits development of language [START_REF] Kral | Neurocognitive factors in sensory restoration of early deafness: a connectome model[END_REF], deafness at mid-age is the main modifiable factor for dementia [START_REF] Livingston | Dementia prevention, intervention, and care[END_REF], and profoundly deaf rehabilitation by cochlear implant (CI) in elderly reduces dramatically cognitive impairment and its evolution to dementia [START_REF] Mosnier | Long-term cognitive prognosis of profoundly deaf older adults after hearing rehabilitation using cochlear implants: cognitive prognosis after hearing rehabilitation[END_REF]. Multichannel CI available since the late 1980s has changed the life of severe to profoundly deaf people [START_REF] Wilson | Cochlear implants: a remarkable past and a brilliant future[END_REF], who would further gain benefit from new technologies by upgrading the external sound processor (SP). New SP generation has improved dramatically CI performance in noise intelligibility for most patients [START_REF] Mosnier | Benefits from upgrade to the CP810 Sound Processor for Nucleus 24 Cochlear implant recipients[END_REF][START_REF] Mosnier | Benefit of the UltraZoom beamforming technology in noise in cochlear implant users[END_REF][START_REF] Franco-Vidal | The Oticon Medical Neuro Zti cochlear implant and the Neuro 2 sound processor: multicentric evaluation of outcomes in adults and children[END_REF]. However, some oldest CI implanted more than 20 years ago could not be upgraded with new SPs due to the lack of compatibility for some brands of CI. It raised the question of surgical reimplantation of the internal part in this population which would have concerned at most more than 600 recipients implanted before 2000 in France. The aim of this study was to show the benefit of a recent generation of SP in the oldest CI population.

Methods

A monocentric study was prospectively conducted on CI recipients, implanted between 1989 and 1997 with the first generation of CI (CI22M, Cochlear, Sydney, Australia), to test the benefit of the more recent compatible SP (CP 900). Among 63 implanted patients, 17 were either lost of follow-up or deceased, 11 others were not eligible for a reimbursement of a new SP, and 2 had incomplete data.

Finally 33 daily CI22M users for 21±2.2 years (range: 19-26) were included. Among them, 16 pre-lingual, 2 peri-lingual, and 15 post-lingual subjects with a profound to total hearing loss were implanted at 5 years (range: 3-16), at 37 and 38 years, and at 36 years (range: 16-60), respectively. Demographic data of the population are described in Table 1. The current SP at the time of the upgrade was a Freedom and an ESPrit 3G SP for 31 and 2 subjects respectively, and all of them used the SPEAK coding strategy. Characteristics of the three generations of SP were indicated in Table 2. The CP 900 SP was fitted with the same SPEAK strategy and the same MAP characteristics as was used in the original processor. Three programs were available: 2 new automatic signal processing algorithm programs (SCAN) and a non-SCAN program with the same options used as in the original SP (Table 1). CI users were tested with their current revised SP, and at 2 months post-upgrade with the CP900 with their preferred program, in a soundproof room with two lists of monosyllabic words presented at 60 dB SPL in quiet (Lafon list) and for the best performers, with an automated adaptive test in quiet and noise (Framatrix [START_REF] Jansen | Comparison of three types of French speech-in-noise tests: A multi-center study[END_REF]), with the current revised SP and the CP900 at 2 months post-upgrade. Subjective perceptions were evaluated using the Abbreviated Profile of Hearing Aid Benefit (APHAB) questionnaire [START_REF] Cox | The abbreviated profile of hearing aid benefit[END_REF]. The results for each test session were compared independently. Scores for words in quiet, the Framatrix test and the APHAB questionnaire were not normally distributed, so a non-parametric Wilcoxon paired test were used. A p value of less than 0.05 was considered to be significant.

Results

Two months after upgrading first generations of internal receiver (CI22M) with the CP900 SP, 31 out of the 33 patients used the SCAN programs all day long (Table 1).

No change was observed in quiet for median scores of words (41% vs 43%) and phonemes (72% vs 67%) recognition of the monosyllabic words list. Using adaptive test,speech reception threshold measurement in quiet was feasible before upgrading, in 18 out of 33 CI 22M recipients in quiet (54%) and 13 out of these 18 patients in noise (39%). In this population, the median speech reception threshold was decreased by 6dB (range: -27, +5; p<0.05) with CP900 as compared to previous SP (Fig. 1A). In noise, it decreased with CP900 by 5.3dB (range: -16, -0.4; p<0.0005) as compared to older SP, and even became measurable at +12.4dB in 6 other subjects (Fig. 1B). The benefit in noise was observed in pre-lingual and in post-lingual CI22M recipients. APHAB scores were available before and after upgrading in 30 out of 33 CI 22M recipients and was unchanged with the CP900 SP (Fig. 2).

Discussion

Over the last 15-20 years, the series of SPs from Cochlear Limited has introduced a number of refinements designed to enhance CI performance mainly in noise. The processor CP810, introduced with the 5 system implant was a technological breakthrough for hearing in noisy conditions by an electrical stimulation of the deaf cochlea. Being compatible with previous CI, upgrading Nucleus® 24 with CP810 noise program improved by more than 20% performance in 77% of users, thanks to the two adaptive omni-directional microphones and new front-end processing options [START_REF] Mosnier | Benefits from upgrade to the CP810 Sound Processor for Nucleus 24 Cochlear implant recipients[END_REF]. In the Nucleus 6 (CP 900 series), an automatic signal processing algorithm (SCAN) has been added allowing automatic transition between six scenes based on the analysis of environmental signal. Because CI is a highly reliable implanted device (cumulative survival percentage of 92.1% over 29 years for CI22M reported by the company) [START_REF]Cochlear ™ Nucleus ® reliability report[END_REF], reimplantation would not be necessary to upgrade ancient CI22M devices which have been implanted more than 20 years ago with new technologies if also compatible. Once made possible, upgrading the old CI with the compatible new generation CP900 SP, yielded to significant improvement in noise intelligibility for more than half of CI22M users in a short period of time. The number of CI22M implantees who reached the level of intelligibility in noise using the Framatix test increased from 39% to 58% with previous SPs and CP900 SP, respectively and most patients prefer new automatic algorithms. No change was observed in quiet regardless of the test used. These results confirm a study including two groups of 15 and 24 CI22M recipients recruited from six North American clinics, showing no change in quiet for CNC words but speech intelligibility improvement in noise with the CP900 compared with the Freedom SP, for AzBio sentences at fixed signal-to-noise ratio [START_REF] Biever | Upgrade to Nucleus 6 in previous generation Cochlear sound processor recipients[END_REF]. This suggests that after such a long duration of cochlear implantation, the frequency mapping should be located in the auditory cortex which may allow a rapid adaptation to new coding strategy and audiological technology. However, the patients implanted more than 20 years ago were different from those who have been implanted more recently and have already benefit from new technologies. Early CI candidates with a prelingual deafnesswere implanted later than nowadays, a critical period for hearing maturation [START_REF] Bruijnzeel | Evaluation of pediatric cochlear implant care throughout Europe: Is European pediatric cochlear implant care performed according to guidelines?[END_REF]. Further post-lingual CI users had a more profound deafness and for longer time than during the last 10 years [START_REF] Holder | Current profile of adults presenting for preoperative cochlear implant evaluation[END_REF]. This may explain that, before upgrade, only 54% of CI recipients in quiet and 39% in noise were able to perform adaptive tests which are more sensitive to show a benefit with the new SP but also more difficult. Moreover, the benefit in noise of upgrading CI22M with CP900 appears to be less striking than upgrading Nucleus® 24 with CP810 [START_REF] Mosnier | Benefits from upgrade to the CP810 Sound Processor for Nucleus 24 Cochlear implant recipients[END_REF] and APHAB questionnaire was probably not sensitive enough to capture a weaker subjective benefit.

Conclusion

More than half of the upgraded old CI recipients with a compatible recent SP performed better the difficult Framatix test in quiet and noise in a short period of time. Although it was thought that reimplantation of the internal receiver should be necessary several times during life-spanning of the majority of profound to total deaf patients, this study evidences that cochlear implantation would be performed only once, especially in children operated on nowadays before 1 year old. Autosensitivity: automatic sensitivity control which reduces the sensitivity of the microphone based on the level of detected background noise; Whisper: algorithm which extend the lower limit of the IIDR in order to improve access to quiet speech; ADRO: adaptive dynamic range optimization algorithm which regulates individual channel gains to improve comfort and intelligibility Beam: adaptive algorithm which modifies its polar characteristics according to direction of the dominant sound source; Zoom: strong directional response, similar to Beam, but fixed rather than adaptive SNR-NR (Signal-to-Noise Ratio based Noise Reduction): attenuate steady-state background noises irrespective of the direction to reduce instantaneously background noise levels WNR (Wind Noise Reduction): algorithm to reduce the low frequency noise from wind SCAN: automatic analyze of the environmental signal classified in 6 scenes (quiet, noise, speech, speech in noise, wind or music) allowing an automatic transition between these 6 scenes.

Speech processors

Esprit 3G

Figure 1 :

 1 Figure 1: Performance using Framatrix in quiet (A, n=18) and in noise (B, n=19) of

Figure 2 :

 2 Figure 2: APHAB scores of CI 22 M recipients upgraded with CP900. The box plots

Table 2 : Sound processing features of the processor types evaluated in the study.

 2

No funding. Written informed consent was obtained from each subject before enrolment into the study (CNIL n° 2040853v0 ).


Benjamin Vincent
email: benjamin.vincent@cnrs-orleans.fr

Sedina Tsikata
email: sedina.tsikata@cnrs-orleans.fr

George-Cristian Potrivitu

Laurent Garrigues

Gaétan Sary

Stéphane Mazouffre

George-Cristian Potrivi

Ga Étan Sary

St Éphane Mazouffre

Electron properties of an emissive cathode: analysis with incoherent Thomson scattering, fluid simulations and Langmuir probe measurements

Keywords: hollow cathode, numerical modeling, incoherent thomson scattering, langmuir probe, thermionic emission, fluid simulation, plasma diagnostics

teaching and research institutions in France or abroad, or from public or private research centers.

Introduction

Hall plasma thrusters, widely used for in-space propulsion of telecommunications satellites, require the use of external emissive cathodes. These act as electron sources for the ionization of the neutral gas propellant and for downstream neutralization of the resulting ion beam. Research into the operation modes and design optimization of emissive cathodes for thrusters has been ongoing for the past several decades [START_REF] Csiky | AIAA 7th Electric Conf[END_REF][START_REF] Csiky | [END_REF][3]. Recent work concerns not only cathodes in the high-current (100 A) range, destined for thruster powers of several tens of kW [4,5], but also low-power designs [START_REF] Nürmberger | AIAA/SAE/ASEE Joint Conf[END_REF][START_REF] Lev | [END_REF][8] for miniature thrusters, suitable for micro-and nanosatellite applications.

In recent years, a major research effort in simulations, experiments and theory has been devoted to the study of the crossed-field discharge of the main Hall thruster plume. This has allowed the analysis of features such as plasma waves [9] which potentially affect thruster operation. Technological advances, such as the development of magnetically-shielded thruster architectures [10], offer the opportunity of extending thruster lifetimes significantly.

Increased thruster development efforts have contributed to further research into cathodes, known to be the lifetimelimiting factor for Hall thrusters. Emissive cathodes are subject to bombardment and erosion by accelerated ions [11], a phenomenon potentially linked to the development of ion-acoustic turbulence in the cathode plume [11][12][13]. Coupling between the emissive cathode and the main thruster plasma plumes [14] may affect global features such as discharge oscillations. Understanding such complex cathode physics and its influence on thruster operation requires reliable diagnostics for particle properties and the development of advanced cathode codes.

Langmuir probes [15,16] are conventional diagnostics which have long been used for characterization of different cathode plasmas, such as in references [17][START_REF] Siegfried | 13th Int. Electric Conf[END_REF][START_REF] Krishnan | [END_REF]. The shape of the current-voltage characteristic obtained from an electrode inserted in the plasma is used to extract electron density and temperature, floating potential, and plasma potential. The simplicity of this diagnostic is associated with drawbacks, such as discharge perturbation [20] and restricted spatial resolution (dependent on the probe dimensions) and temporal resolution.

Incoherent Thomson scattering (ITS) is less-commonly used for electron property characterization in low-density plasmas. This diagnostic technique uses the intensity and spectral shape of the light elastically scattered by the free electrons of the plasma. The main challenges of the diagnostic are its complexity (relative to probes) and the low scattered intensity available at low plasma densities, but it has the advantage of avoiding plasma perturbation, if certain criteria are respected [21,22] and high spatial and temporal resolution are attainable.

In addition to experimental tools, numerical simulations are the source of insights into cathode physics. Early analytical studies [23][24][START_REF] Salhi | th Joint Conf. and Exhibit[END_REF] on aspects of hollow cathode physics such as electron emission gave rise to increasingly sophisticated numerical models, and in particular, the first comprehensive 2D axisymmetric code OrCa2D [START_REF] Mikellides | [END_REF], which has been applied in recent years to the study of several features, including cathode erosion and the development of oscillations and spotplume regime transitions [27][28][START_REF] Mikellides | AIAA/SAE/ASEE Joint Conf[END_REF]. Over the last few years, another 2D axisymmetric fluid code [START_REF] Sary | [END_REF]31] was developed and evaluated for the simulation of thermal and plasma physics of the NSTAR NASA hollow cathode. This code proved capable of accurately reproducing profiles for the interior and exterior electron properties measured using probes.

The availability of a sensitive, non-invasive and reliable ITS diagnostic offers the opportunity to evaluate the performance of numerical codes and determinine their ability to reproduce typical cathode operation and electron property profiles.

With this in mind, in this work, we adapt the Sary and Garrigues fluid code [START_REF] Sary | [END_REF]31] to studies of a 5 A MIREA cathode equipped with a LaB 6 pellet emitter, and investigate a range of operating regimes. The electron properties from these regimes are compared directly with ITS electron property measurements on the same cathode. This, to our knowledge, is the first work combining ITS and numerical simulation studies for emissive cathodes. A limited comparison of results to Langmuir probe measurements is also discussed.

This paper is structured as follows. The emissive cathode source and diagnostics (ITS, Langmuir probe) used for the investigations are desribed in section 2. Section 3 is devoted to the presentation of the main features of the 2D fluid model. In section 4, diagnostic and simulation results for electron properties, obtained under various conditions, are presented and discussed. A summary of the main insights gained is presented in section 5.

Experimental setup

Cathode details and operation

The cathode was tested under vacuum in a 1.8 m long, 0.8 m diameter stainless steel vacuum vessel. A 110 l s -1 dry pump combined with a 350 l s -1 turbomolecular pump and a cryogenic pump cooled down to approximately 30 K gave a base pressure of 10 -5 Pa and an operating pressure below at least 4 × 10 -3 Pa for all the cathode discharge operating conditions.

The cathode considered in this work is based on the model developed by the MIREA Institute and features a thermoemissive LaB 6 pellet heated by a tungsten filament to supply electrons to the discharge. The cathode internal body, made of molybdenum, minimizes radiative losses dissipated by the tungsten filament. In addition, the high electrical conductivity of the internal body and LaB 6 pellet guarantees electrical closure of the heating circuit connections shown in red in figure 1. The pellet is a disc 8 mm in diameter and 2 mm thick, however, the emissive surface is reduced to 6 mm in diameter due to the internal body architecture.

Both the experimental campaign and simulation focus on the behavior of the cathode operation in diode mode, operated independently of a thruster and using an externally-mounted anode. For the experimental campaign, the cathode was operated in self-heating mode. The experimental setup uses a stainless steel disc anode to establish the potential drop through which the electrons are accelerated, of diameter 110 mm and with a central hole of diameter 21 mm. This anode directly faces the 3 mm-diameter cathode orifice and is located at an axial distance of 40 mm. Figure 1 shows the basic architecture of this cathode and anode mounting. The cathode is operated over a range of discharge currents (2-16 A) and xenon flow rates (0.2-0.8 mg s -1 ).

Two operating modes named the 'spot' and 'plume' modes are commonly identified for cathode discharge operation in the diode configuration [START_REF] Csiky | AIAA 7th Electric Conf[END_REF]. The spot mode is usually the operation mode desired and is associated with high electron density, low electron temperature (Csiky identified this mode with temperatures in the range of 0.5 eV [START_REF] Csiky | [END_REF]) and low discharge oscillation. The plume mode is associated with lower electron density, higher electron temperature (in the range of 1 -2 eV in reference [START_REF] Csiky | [END_REF]) and higher oscillations of the discharge voltage and current (a ratio I d,std /I d,mean above few percent is usually attributed to the plume mode [32]), leading to a lower cathode lifetime [33]. The experimental work presented in this article focuses on discharge conditions for which the cathode discharge is in the so-called spot mode.

Incoherent Thomson Scattering diagnostic

A complete description of the design of the incoherent Thomson scattering diagnostic, as well as the procedures used for signal detection and calibration, can be found in detail in reference [34]. This section makes a summary of the main features of the diagnostic and the specificities related to its implementation on the cathode.

Summary of diagnostic design.

Electron property measurements were performed using the recently-developed incoherent Thomson scattering diagnostic named THETIS [34]. A key feature of this diagnostic is the use of a recently-developed volume Bragg grating (VBG) notch filter to efficiently limit stray light transmission near the incident laser wavelength, while preserving a high transmission efficiency outside the band-stop region. This diagnostic design allows a simpler detection branch in comparison to diagnostics using triple grating spectrometers and significantly reduces transmission losses, resulting in a high degree of sensitivity.

The diagnostic setup is shown in figure 2. It is comprised of (i) a transmission branch, used to produce the Thomson scattering signal (this branch is shown in dark green in figure 2), and (ii) a detection branch, used to collect the scattered light over a reduced solid angle, filter the sources of stray light, spectrally disperse the photons collected, and measure the scattered signal with an intensified CCD camera (this branch is shown in light green in figure 2).

The light source of the transmission branch is a 10 Hz pulsed, frequency-doubled Q-switch Nd:YAG laser (wavelength of 532 nm, pulse energy of 430 mJ and pulse width of 5 ns). The initial vertically polarized 9 mm diameter beam is expanded with a magnification of 3 and focused with a 2 m-focal length lens in the plasma volume. To reduce stray-light propagation the beam passes through two Brewster windows mounted at the end of tubes inside which apertures are placed. The beam ends its path in a large-aperture beam dump.

Scattered light is collected by the detection branch and imaged onto the entrance of a 5 × 3 matrix of fibers (0.3 mm in individual diameter and NA = 0.22). Light is focused onto the fibers with a 200 mm focal length, 100 mm diameter biconvex lens positioned at 520 mm. The fiber bundle conducts the light to the optical bench where it is collimated and filtered by the VBG notch filter. The filtered light is then transmitted to the spectrometer through a 0.3 mm wide slit and dispersed with a 2400 lines mm -1 holographic grating. An intensified camera, triggered on the laser Q-switch signal, measures the spectrum of the remaining photons over a 15 ns temporal window.

A translation stage allowed measurement from 1 mm to 10 mm from the orifice. Due to the collection optics magnification, the size of the fiber bundle and the size of the laser spot, the observation volume over which electron properties are averaged can be considered as a cylinder 0.3 mm in diameter and 2.5 mm in length. Given the orientation of the detection branch relative to the transmission branch, electron properties are probed in the (yz) plane at an angle 45

• relative to the y direction.

Procedure for data analysis.

To obtain quantitative data, the detection branch of the THETIS diagnostic was calibrated with the same acquisition parameters as for Thomson scattering measurements.

The spectral calibration was performed with a low-pressure neon gas lamp. The instrument function (which characterizes the light redistribution induced by the detection branch) was estimated using the spectrum measured from the background signal generated by the laser reflections. These sources of light from the laser were considered perfectly monochromatic in comparison to the spectrometer resolution. These spectral calibrations were performed with the VBG filter removed, and used to quantitatively analyze the spectral distribution of the signals acquired.

The intensity calibration was performed using the fitting of the Raman scattering spectrum, measured from a known pressure of nitrogen, and with the energy of the laser, monitored after each acquisition. The VBG notch filter was tuned for stray light rejection at 532 nm. With this calibration, the total intensity of the signal from a given scattering process can be used to determine the density of the scattering particles when the corresponding scattering cross-section and intensity of the incident light are well-characterized.

Examples of Thomson scattering spectra obtained from cathode experiments for two different discharge currents are shown in figure 3, with corresponding electron properties indicated in the caption. From these spectra, the electron temperature is deduced from the spectral width and the density is estimated from the total intensity scattered. These spectra show the spectral region of rejected data (corresponding to the residual stray light within the notch filter bandstop, in orange), the full set of experimental points used for the analysis procedure (in blue) and the green fit line from which Gaussian fit parameters are determined to deduce the electron properties. As is clear from figure 3(a), the high sensitivity of this diagnostic allows detection of Thomson signal amplitudes of only a few counts, which corresponds to electron densities in the 10 16 m -3 range.

Langmuir probe

Probe measurements were performed with a single cylindrical Langmuir probe made with a 0.10 mm-diameter tungsten wire, protruding 1 mm from a single-bore alumina tube. Measurements were performed on the cathode axis, with the probe oriented along the x-direction (cf figure 2). In addition to the three-axis displacement stages of the cathode, the probe was placed on a translation plate for quick displacement along the x-axis. The resolution of electron properties along the x-axis with the Langmuir probe is assumed to be approximately the size of the probe itself, with the influence of the sheath and imprecision of translation stages neglected. The acquisition of the I -V characteristics was achieved with an Impedans ALP system. A 1:11 attenuator was used when the collected current level exceeded that of the analog-to-digital converter dynamics.

Because the plasma outside the cathode orifice is collisionless (electron mean free path much greater than the probe dimensions) and unmagnetized, the Druyvesteyn method [35] was used to extract information on electron properties from the I -V characteristic. To do so, the I(V) characteristic was first filtered with a fourth order Butterworth filter and then smoothed using a Savitzky-Golay filter [36]. The plasma potential was obtained from the root of the second derivative of the I(V) curve (calculated by the Savitzky-Golay filter). The Druyvesteyn formula was applied to this second derivative to obtain an EEDF (not normalized):

F(ϵ) = 4 e 2 A p ( mϵ 2e 2 ) 1/2 d 2 I dV 2 (J -1 m -3 ) (1)
where A p is the probe collection area, ε = e(V p -V) the electron energy and V p the plasma potential estimated with the position of the maximum of the first derivative of the I(V) function.

The normalization coefficient of equation ( 1) within the interval boundary starting from zero electron energy (when the probe potential is at plasma potential) up to the maximum electron energy (in principle, obtained at the minimum probe potential scanned) give the electron density:

n e = ˆ+∞ 0 F(ϵ)dϵ(m -3 ) (2)
Using the hypothesis of a Maxwell-Boltzmann Electron Velocity Distribution Function (EVDF) and isotropic electron behavior, the normalized EEDF f (ε) = F(ε)/n e is a chi-squared distribution with three degrees of freedom. The calculation of the first moment of this distribution gives:

E[ϵ] = ˆ+∞ 0 ϵf(ϵ)dϵ = 3 2 k B T e (J) (3)
Under the hypothesis of a chi-squared experimental EEDF we obtain :

T e = 2 3k B E[ϵ](K) (4)

Uncertainty analysis

Systematic and random uncertainties may be considered relevant to this work. For the uncertainty analysis of the experiments presented, only systematic uncertainties were taken into account. To determine the random uncertainty, measurements would have to be repeated many times. For example, due to fast probe aging in the plasma studied, a new probe for each Langmuir measurement would have to be used; this would not have been feasible within the restricted time available. For the Thomson measurements, each spectrum requires 20 min of acquisition (10 min for the Thomson scattering signal and 10 min for the background signal accumulations for subtraction). The time constraint is again the limiting factor. The uncertainties were calculated using a propagation of error calculation involving the quadratic sum of systematic uncertainties in the parameters used to calculate n e and T e .

Thomson scattering uncertainties.

Uncertainties in T e and n e from Thomson measurements have few contributions from the experiment equipment. The main experimental contribution is from laser power root mean square (r.m.s) fluctuations. Statistics on the laser energy measured over 10 min during background acquisitions confirmed an r.m.s variation below 1%. The use of a high precision Inficon CDG025D capacitive diaphragm pressure gauge enabled a relative precision of 0.2% in the pressure measurement. The two main sources of uncertainties are from the data analysis phase. For data analysis of the calibration spectrum, a Raman cross-section with an uncertainty of 8% is used [37]. In addition, each spectrum fitting is obtained with an uncertainty whose value is linked to the signal-to-noise ratio level; to take this contribution into account, the covariance matrix was extracted for each spectrum fit.

Langmuir probe uncertainties.

For probes, accounting for errors associated with the EEDF method used to determine electron properties is a rather complex task. Herman and Gallimore [START_REF] Herman | AIAA/ASME/SAE/ASEE Joint Conf. & Exhibit[END_REF] predict an error of ± 8% in the determination of the EEDF using the Druyvesteyn method, but only when all the constitutive parameters that enter into the EEDF computation are known to a high degree of accuracy. The plasma potential magnitude remains the main source of uncertainty in computing the EEDF. The characteristic times and filtering circuits of the Impedans ALP System, together with data processing (filtering, numerical differentiation) induce systematic uncertainties in the plasma property estimations. If the first derivative of the probecollected current is used for determining the plasma potential, the accuracy of the potential measurement is approximated to ± 10% [START_REF] Popov | [END_REF]. The following results use an uncertainty in plasma density of ± 50% and ± 30% for the electron temperature [START_REF] Reid | AIAA/ASME/SAE/ ASEE Joint Conf. & Exhibit[END_REF].

Numerical model

As mentioned in the introduction, this work seeks to compare numerical results on electron properties for the 5 A MIREA cathode with properties measured directly in experiments. A particular focus is placed on ITS as a new diagnostic. This work follows (i) application in 2017 of our code [START_REF] Sary | [END_REF]31] to a different standard cathode, and (ii) recent validation of the performance of the ITS diagnostic [34]. This research seeks to build on the existing body of experiment-versus-simulation comparisons, many of which have been performed in studies using the well-established OrCa2D [START_REF] Mikellides | [END_REF] code. Full details on the numerical method are presented in reference [START_REF] Sary | [END_REF], and are only briefly summarized here.

Governing equations

The numerical scheme considers three species: Xe atoms, Xe + ions, and electrons, with the multiply-charged ion population are neglected. Maxwellian velocity distributions are considered to apply for all species in this plasma. Inside the cathode, the transport of species is dominated by collisions between electrons and neutrals (e -n), between ions and neutrals (i -n), and Coulomb collisions between electrons and ions (e -i). The plasma is assumed to be quasineutral with the ion density n i equal to the electron density n e .

The mass conservation equation for the ions can be written as:

∂n ∂t + ∇ • (n⃗ u i ) = S (5)
where ⃗ u i is the ion velocity. The mass conservation equation for the neutrals can be written as:

∂n n ∂t + ∇ • (n n ⃗ u i ) = -S (6)
where n n is the neutral density. In equations ( 5) and 6, S = nn n k io (T e ) is the source term and k io the ionization rate coefficient.

The momentum conservation equations for Xe + ions, neglecting viscosity and magnetic field, are written:

M ( ∂ ∂t (n⃗ u i ) + ∇ • (n⃗ u i ⊗⃗ u i ) ) = -ne∇ϕ -∇p i + MS⃗ u n (7)
where ϕ is the plasma potential.

The following electron drift diffusion equation is used:

⃗ j e = -enµ e ∇ϕ + µ e ∇p e (8) with the electron pressure p e = nk B T e and the electron mobility µ e = e/mν e , and charge conservation is written as:

∇ • ( ⃗ j i + ⃗ j e ) = 0 (9)
giving the expression for the plasma potential

∇ • (enµ e ∇ϕ) = ∇ • (µ e ∇p e + en⃗ u i ) (10)
Lastly, energy conservation for the different species is written as follows: for neutrals:

∂ ∂t ( 3 2 p n ) + ∇ • ( 5 2 p n ⃗ u n +⃗ q n ) = ⃗ u n • ∇p n - 3 2 k B T n S -⃗ u s • ∇ • ⃗ T + ∇ • ( ⃗ T⃗ u s ) +Q n -( ⃗ R iso,i→n + ⃗ R e→n ) •⃗ u n ( 11
)
where the neutral thermal conduction flux is

⃗ q n = -λ n ∇T n , Q n = nν iso k B (T i -T n ) + 3 m M nn n k m e,n k B (T e -T n )
and ⃗ R s→t is the momentum exchange due to elastic collisions between species s and t.

for ions:

∂ ∂t ( 3 2 p i ) + ∇ • ( 5 2 p i ⃗ u i +⃗ q i ) = ⃗ u i • ∇p i + 3 2 k B T n S + 1 2 MS(⃗ u i -⃗ u n ) 2 (12)
+Q i -( ⃗ R iso,n→i + ⃗ R e→i ) •⃗ u n ( 13
)
where

Q i = nν in k B (T n -T i ) + 3 m M nν ie k B (T e -T i ) and ⃗ q i = -λ i ∇T i ,
for electrons: (14) where

∂ ∂t ( 3 2 p e ) + ∇ • ( 5 2 p e ⃗ u e +⃗ q e ) = ⃗ u e • ∇p e + Q e + mν e n⃗ u 2 e -Sϵ io
Q e = 3 m M nν ie k B (T i -T e ) + 3 m M nn n k m e,n k B (T n -T e ),
⃗ q e = -λ e ∇T e and ε io is the energy lost per ionization, i.e. 12.1 eV, if only the first ionization energy of xenon is considered.

Simulation domain and boundary conditions

The existing simulation domain is adapted to match that of the MIREA cathode under study. Figure 4 shows the mesh with respect to the cathode dimensions and architecture. The figure shows the cylindrical simulation domain applied to half of the axisymmetric cathode, along with the dimension and placement of the anode (blue) and LaB 6 pellet emitter (red).

The plasma simulation domain is discretized using a Cartesian mesh of 1880 nodes (710 in the internal region of the cathode). Further details on the discretization method are provided in reference [START_REF] Sary | [END_REF]. The mesh size is at its smallest between the emitter and the orifice with a square cell size of 0.215 mm.

For the plasma conditions inside the cathode orifice, the Debye length, at the micrometer scale, cannot be resolved by the mesh. Instead, as discussed in detail in reference [START_REF] Sary | [END_REF], the wall sheaths are modeled using boundary conditions accounting separately for (i) thermionic emission of the LaB 6 pellet, and (ii) the conducting nature of the cathode body.

The calculations of ion and electron fluxes and energy losses deposited on the metallic parts of the cathode (in brown in figure 4) can be treated analytically, assuming a non-collisional sheath with known plasma properties at the sheath edge (plasma density, electron temperature and sheath potential drop) [START_REF] Lieberman | Principles of Plasma Discharges and Materials Processing[END_REF]. We assume that ions lost on the walls are recombined into neutral atoms at the temperature of the wall. A slip-length boundary condition is used for neutrals, which nullifies the neutral velocity parallel to the wall at a distance below the boundary (this distance can be estimated from the fluid viscosity). The normal component of the flux is set to the recombined ion flux.

At the gas inlet (not shown in figure 4), we assume a radial profile of the neutral velocity normal to the inlet given by a Poiseuille profile. A reference density is chosen based on the density in the closest domain cell and the velocity is then computed from the desired value of the xenon mass flow injected. The temperature of the gas corresponds to the cathode body temperature. Regarding the anode plane (in blue in figure 4), a uniform electron current density is assumed (whose integration yields the given electron current at the anode). Ions are accelerated at the Bohm velocity and recombine as atoms. The external cathode body is floating and its potential is determined through a balance of positively-and negatively-charged current that impacts its surface. A complete description of the specific treatment of the open boundary of the plume region, which is not a physical boundary, can be found in reference [START_REF] Sary | [END_REF].

Figure 5 shows 2D maps of electron temperature and density obtained from the simulation of the cathode discharge with two different thermionic emitter temperatures (1900 and 2050 K ). As shown in reference [31], similar maps can be obtained for the neutral density, neutral temperature, ion temperature and plasma potential, for example. The maps of figure 5 correspond to snapshots of the steady state solutions obtained after a simulation run time of about 5 ms. The color scale for density and temperature has been adapted in order to highlight variations in the electron properties outside the cathode orifice, where the laser scattering measurements were performed.

Thermionic emission description

Particular attention must be paid to the treatment of thermionic emission of electrons from the emitter (in red in figure 4). The electron current density emission is given by the Richardson-Dushman law:

J em e = A 0 T 2 w exp ( -eϕ b k B T w ) ( 15
)
where A 0 = 119.58 A cm -2 K -2 is Richardson's constant, T w is the emitter temperature, and ϕ b = ϕ WF -ϕ SH a potential barrier that includes ϕ WF , the temperature-dependent work function, and ϕ SH , the work function lowering due to the Schottky effect.

To calculate ϕ b , a self-consistent model has been established to describe the continuous transition from the space charge-limited to thermionic emission regimes that can occur in the sheath in front of the emitter. This model, detailed in reference [START_REF] Sary | [END_REF], is an extension of the work of Lin and Eng [START_REF] Lin | [END_REF].

It requires, as input plasma properties at the sheath edge, electron density n e , electron temperatures T e and T w (of the bulk plasma and at the wall, respectively), electric potential ϕ 0 , and electron emission at the walls.

For the electron emission, the temperature dependence of the work function ϕ WF (T w ) uses the experimental data of reference [12]:

ϕ WF (T w ) = 2.66 + 1.23 × 10 -4 T w [eV] (16)
In practice, we calculate a priori the emitted electron current density J em e as a function of the emission model parameters (n 0 , ϕ 0 , T e , T w ) for a large range of parameter variations. During the 2D simulations, the emitted electron current density is simply extrapolated from pre-tabulated values.

In this work, we studied the impact of the emitter temperature on plasma properties.

In recent numerical studies of Sary et al, two cases were considered: (i) a validation case, first using an imposed emitter temperature profile [START_REF] Sary | [END_REF] (an uncoupled plasma-thermal model), based on experiments on the NSTAR cathode performed by Katz et al [43], and (ii) a coupled plasma-thermal model [31].

reference [START_REF] Sary | [END_REF] introduced a new model to describe electron emission in both the space charge limited regime and thermionic emission limited regime. The model of reference [START_REF] Sary | [END_REF] proved capable of reproducing electron properties for the imposed temperature profile.

In this work, a slightly different approach to the work of references [START_REF] Sary | [END_REF]31] is taken. Firstly, as the exact thermal design of the MIREA cathode was not available, we opted not to investigate a full coupled plasma-thermal model as yet. Secondly, for the uncoupled model, rather than imposing a temperature profile for the emitter, we imposed a single uniform temperature for this element. This simplified approach is, to a certain extent, justified when one considers the very different architectures of the NSTAR cathode (which uses an emissive cylinder which develops a natural temperature gradient along its length) and the present MIREA cathode, which uses a simple emissive pellet.

The emitter temperature is treated as a parameter indirectly characterizing the emission properties of this cathode. In the simulation, the variation of this temperature is used to scan different modes of cathode operation and electron properties. This variation is performed using a lower bound for the emitter temperature which was previously measured with a thermocouple at similar discharge conditions.

In recent cathode numerical studies, Guerrero et al [START_REF] Guerrero | 36th Int. Electric Conf[END_REF] have used the OrCa2D plasma code of reference [START_REF] Mikellides | [END_REF] coupled self-consistently with a thermal model derived using COM-SOL. In order to evaluate the performance of our plasmathermal code, a thermal model for the MIREA should be developed in the future.

Ion acoustic instability

In the NASA NSTAR cathode, the experimental evidence of fluctuating potential responsible for ion backflow on the keeper cathode has been linked to an enhanced resistivity in the cathode plume, first proposed in reference [START_REF] Mikellides | [END_REF]. The drop of electron density induces an increase of electron velocity to maintain the desired electron current. As a result, the relative velocity between electrons and ions can exceed the ion acoustic velocity, exciting the ion acoustic instability (IAI). Sary et al have implemented a numerical model of the acoustic instability for the cathode plume (largely inspired by the work of Davidson and Krall [46]). The IAI is modeled through an energy equation that is solved for the fluctuating electric field ε f :

∂ε f ∂t + ⃗ ∇ • ( ⃗ V d | ⃗ V d | c s,i ε f ) = 2γ M ε f ( 17
)
where ⃗ V d = ⃗ u e -⃗ u i is the drift velocity and c s,i = √ kBT e /M is the ion acoustic velocity. γ M the maximum growth rate (only the fastest-growing mode is retained and detailed in reference [START_REF] Sary | [END_REF]). Assuming that the Mach number M e = | ⃗ V d |/c s,i ≫ 1, and γ M ≪ ω pe , ω pe being the electron plasma frequency, we can then link the anomalous electron collision frequency to the electron temperature, plasma density, and fluctuating electric field ε f :

ν eAN = 2 √ πω pe ε f nk B T e ( 18
)
Thermal terms in the energy conservation equation of ions and electrons have also been included. The system of equations here described is valid only for the description of the growth rate of the IAI in the quasi-linear regime, when ion and electron energy distribution functions are almost Maxwellian. It cannot describe the saturation mechanism of the instability due to kinetic effects. Using a derivation from particle-incell simulations, Davidson and Krall have proposed to limit the maximum of ε f [46]:

[ε f ] MAX = α nTe nk B T e , where α nTe ∈ [10 -3 , 10 -2 ] (19)
Results presented in reference [31] for a BaO NASA NSTAR cathode have shown that α nTe = 3 × 10 -3 gives a reasonable agreement with experiments. We use the same value for this study.

The time-averaged values of electron temperature simulated in reference [31] were weakly affected by the saturation threshold of the acoustic ion instability. With a higher saturation threshold, the electron temperature values oscillated more strongly, but the time-averaged values remained unchanged. The focus of this work is temporally-averaged values of electron parameters, for which we do not expect the ion-acoustic instability effect on temperature to be necessarily visible.

A brief discussion of the influence o f t he e mitter temperature on the appearance of discharge current oscillations compatible with transit time oscillations and the ion acoustic instability can be found in reference [START_REF] Joussot | 35th Int. Electric Conf[END_REF]. Qualitative agreement between experiments and calculations is also reported in reference [START_REF] Joussot | 35th Int. Electric Conf[END_REF] when the cathode is operated in non-self-heated mode. It is important to note that the fluid treatment of the species in this code would limit, to some degree, the ability of the model to quantitatively describe features in the plume.

Results and discussion

Electron properties from ITS measurements

Electron density and temperature are determined for various flow rates and discharge currents. Generally, as the dense plasma created inside the cathode expands outside over a wide solid angle, a decrease of the neutral, ion, and electron density is expected, along with the electron temperature.

Xenon flow rate exploration.

Axial measurement of electron properties with a discharge current of 2 A and xenon flow rate between 0.2 mg s -1 and 0.8 mg s -1 were made using the ITS diagnostic. Results are shown in figure 6. It is observed that an increase of xenon flow i nduces an increase of electron density. This increase is more pronounced close to the orifice and can be explained under the hypothesis of a steady ionization rate. At higher xenon flow rate, neutral density, and consequently electron density, rise if the ionization rate remains constant.

For the electron temperature, a slight decrease is observed at a higher xenon flow r ate. T his i s d ue t o a n i ncreased rate of electron collisions with neutrals or ions with higher xenon flow.

Discharge current exploration.

Axial measurements were performed with a fixed xenon flow rate of 0.8 mg s -1 and discharge currents between 2 A and 16 A. Results are shown in figure 7.

Here, an increase in discharge current induces an increase in electron temperature. At low currents (2 A to 6 A) electron temperature decreases monotonically with distance. However, for higher currents (6 A to 16 A), this trend changes: an increase of electron temperature at distances greater than 5 mm is observed. By extrapolation, for distances greater than 10 mm (where electron density is too low to probe electron properties using the diagnostic) it might be expected that at low discharge current, electron temperatures would remain low ( < 0.5 eV), while at high discharge current electron temperatures might reach high values ( > 2 eV).

The separation between the spot and plume regimes has been shown to be dependent on a number of factors, including discharge parameters and anode distance [START_REF] Csiky | [END_REF]32,[START_REF] Hall | 36th Int. Electric Conf[END_REF]. The spot-plume transition may not necessarily be a sharply-defined boundary, as shown by Csiky [START_REF] Csiky | [END_REF], and the two regimes may be separated by a transition regime, as shown by Kaufman [START_REF] Kaufman | Technology of electron-bombardment ion thrusters[END_REF]. In characterization of these regimes, we will be guided by the trends in electron properties observed from the ITS measurements and also, discharge metrics. The plume mode is associated with regimes of (i) higher discharge current/flow rate ratios, and (ii) higher voltage oscillations. Recently, Mikellides et al investigated transitions between these regimes in the OrCa2D code on the basis of similar metrics [START_REF] Mikellides | AIAA/SAE/ASEE Joint Conf[END_REF].

The observation of the highest electron temperatures at some distance from the orifice is compatible with the plume mode for 8 A and 16 A discharge currents. The observation of a monotonically decreasing electron temperature profile with a maximum near the orifice is compatible with the spot mode.

The plume mode has been associated with a sheath extending well beyond the cathode orifice, while the spot mode is characterized by a contracted sheath primarily within the orifice. Although investigations within the orifice are not accessible in the experimental study, the rate of decrease of the electron density beyond the orifice as a function of axial distance (clearly slower for 16 A than at, for example, 2 A) further suggests that the higher current regime at 16 A has features of a plume regime.

The nature of the regimes may be further investigated by a consideration of the discharge voltage standard deviation. In this work, no measured voltage oscillations exceed 2% of the mean voltage; this was the threshold described in reference [32] for a spot-plume transition. For all the discharge conditions explored experimentally, the discharge appears to remain in the spot mode, according to this criterion. However, other features (such as the trend in the experimentally-determined electron temperature and density of figure 7 at 8 A and 16 A) may develop before the discharge fully transitions from the spot to the plume mode at higher currents.

An increase in electron density is observed at higher discharge currents. As electron current is the main discharge current contribution inside a cathode discharge, the electron current density is expected to follow the total discharge current trend. Because the electron current density is proportional to both the electron density and the drift velocity, its increase is expected to be correlated with the increase of at least one of these quantities. The results presented show that the increase in discharge current is at least due to an increase in electron density. The electron velocity component along the probed radial direction was too low to be resolved and the component along the cathode axis could not be probed. With these experimental constraints, no comments can be made regarding the contribution of electron drift velocity to the electron current density.

Electron properties from 2D simulations

As mentioned earlier, in these simulations, a variation of the LaB 6 pellet temperature is studied. This temperature variation corresponds to a variation of the emission current density [4,12] (see equation ( 15)).

4.2.1.

Emitter temperature exploration. The 2D fluid model was used to simulate the discharge at identical discharge conditions to those explored by ITS. Because the temperature of the LaB 6 emissive pellet is a key parameter for the simulation, and one whose measurement with high precision is challenging [START_REF] Guerrero | 36th Int. Electric Conf[END_REF], particular attention was paid to its influence on the electron density and temperature values outside the cathode. The lowest temperature of the pellet was chosen close to values previously measured in our team with a thermocouple. In the simulation, the temperature was varied by 5 K steps over a temperature range of 100 K (where possible). Discharge currents from 2 to 16 A with a xenon flow rate of 0.8 mg s -1 were simulated.

For each set of discharge parameters simulated, the simulation was run for 5 ms, and convergence was verified. Electron temperature and density along the cathode centerline were extracted from the 2D maps of electron properties calculated by the simulation. Snapshots of these are provided as examples in Fig 5, for two values of the electron emitter temperature. The axial profiles extracted to perform a parametric study were obtained from a temporal average over the last 1 ms of the simulation and from a spatial average over the radial length which corresponded to that probed in scattering experiments (up to 1.2 mm from the cathode axis).

The axial profiles for a selected set of emissive pellet temperatures are shown in figures 8, 9 and 10 for a discharge current of, respectively, 2 A, 4 A and 8 A. The set of emissive pellet temperatures displayed was chosen to maximize readability on the figures. With a 16 A discharge current, the range of tested temperatures was reduced to 50 K because the discharge systematically turned off above 2060 K. The corresponding results are not presented in this paper.

Some general trends apply to the discharge currents shown in figures 8 t o 1 0. T he e lectron d ensity a lways decreases monotonically with increased distance from the cathode orifice. When the temperature of the emissive pellet is increased, electron density outside the cathode orifice progressively decreases before reaching a steady low-density profile at high emissive pellet temperature.

The trend in electron temperature changes as the emissive pellet temperature rises. For the presented results, electron temperature decreases with distance from the orifice a t low emissive pellet temperature. However, at higher emissive pellet temperatures, the electron temperature first increases at distances close to the cathode orifice before approaching a plateau (clearly visible for the cases at 4 A and 8 A). The position at which this plateau starts to develop approaches the orifice a s e mitter t emperature i ncreases. A s s uggested w ith the simulation at 2 A discharge current, it seems that the electron temperature profile returns to a decreasing trend once the plateau has reached the orifice. F or t he 8 A d ischarge condition, the transition from a low T e , high n e regime at low emissive pellet temperature to a high T e , low n e regime at pellet temperature above 1925 K may correspond to a transition of the discharge behavior from a spot mode to a plume mode, if we consider electron density and temperature values as the criteria for judging such a mode transition. As shown in figure 5 , t he r egion w here a l arge e lectron t emperature is predicted fills most of the simulated domain outside the cathode orifice. These large values of electron temperature would probably translate into a luminous large plume outside the cathode orifice, commonly observed with a cathode in plume mode.

In the simulation, the discharge current is fixed a nd the varying discharge voltage U d may be examined for information on the different regimes explored. The discharge voltages obtained in the simulation at 8 A for different emissive pellet temperatures (1900 K, 1925 K, 2000 K and 2050 K) are considered.

The values of mean discharge voltage and standard deviation of these temperatures are, respectively, 44.4(3) V, 43.9(4) V, 39.4(1) V and 39.1(18) V. In experiments, a mean and standard deviation of the discharge voltage of 37.4(5) V were measured at 8 A in self-heating mode; the simulation therefore produces voltage values close to those measured experimentally. We can attempt to classify the discharge mode as either plume or spot using the discharge voltage oscillations, as done earlier using the discharge current (where a plume mode was associated with a standard deviation above 2%). The increase in pellet temperature from 2000 K to 2050 K coincides with a jump in the standard deviation from 0.3% to 4.7%. This suggests a possible spot-plume mode transition occurring above 2050 K. This transition seems to appear at a different emissive pellet temperature than that previously predicted from the n e and T e profiles (i.e. a pellet temperature of 1925 K in figure 10).

The sensitivity of the discharge mode to the emissive pellet temperature is evident in the simulations. Factors contributing to the thermal balance (the power applied to the heating circuit and the heat shield design) are clearly key elements in the design of a cathode for operation in the spot mode, for a specific d ischarge current range.

Comparison with ITS measurements.

The axial profile of electron properties from simulations can first be compared to ITS experimental results for the electron properties obtained for the same discharge current and xenon flow rate. Electron properties derived from the simulation at several emitter temperatures were considered for comparisons.

For a xenon flow rate of 0.8 mg s -1 and discharge currents of 2 A, 4 A and 8 A, the emissive pellet temperatures leading to the best agreement between the electron properties obtained from ITS measurements (shown in figure 7 ) a nd simulation result (shown in figures 8 to 10) are, respectively, 1800 K, 1860 K and 1925 K.

The emitter temperatures leading to the best agreement in electron properties (for the ITS-simulation comparisons) do not correspond to the temperatures that give best agreement in discharge voltage (for ITS-simulation comparisons). In other words, not all the discharge regime features apparent in experiments are captured by the simulations, even though certain behaviors (such as spot-plume transitions) appear to be seen in the simulation. Refinement of the simulation (by selfconsistent coupling of the MIREA plasma and thermal models) will be needed in the future.

Consistently, for the discharge parameters used, the best agreement between simulations and experiments was observed near the transition from a high n e , low T e axial profile t o a low n e , high T e axial profile. The experimental measurements align with a cathode on the verge of a spot-plume transition. Nevertheless, as seen in the previous section, the demarcation between these regimes is not easy to identify using only the amplitude of the discharge voltage oscillations. A high value of T e seems to be observable with a low discharge voltage oscillation level, particularly near the mode transition.

Electron properties from Langmuir probe measurements

The comparisons may be extended with the inclusion of electron property measurements from Langmuir probes. For this, electron density and temperature measurements obtained by Langmuir probes at 8 A discharge current and a xenon flow rate of 8 mg s -1 were compared to ITS and simulation data for the same operating point. Results are shown in figure 11. For these discharge parameters, the discharge voltage oscillations were low (below 2% of the mean discharge voltage) and the discharge was considered to be in the spot mode. The constant nature of the discharge properties was verified with ITS before and after the Langmuir probe measurements.

Comparison with ITS measurements.

With measurement uncertainties taken into account, figure 11 shows that the electron density obtained from the Langmuir probe and ITS measurements are in agreement and exhibit the same trend. Nevertheless, this is not the case for the electron temperature, where temperatures measured with the Langmuir probe are four times the values obtained from ITS measurements.

The discrepancy between the ITS and Langmuir probe temperatures may be due to a number of factors.

One important reason might be the change in the discharge mode when the probe was inserted and polarized inside the plasma. When the probe was inserted and kept inside the discharge, large discharge voltage oscillations were observed. These oscillations, with amplitudes as large as 5 V, were observed when the probe was placed at distances 10 mm and closer to the orifice. At distances 5 mm and below, these perturbations remained present for up to a minute after extraction of the probe from the plasma. This significant perturbation strongly suggests that probe insertion and measurement were associated with a change of the cathode regime from the spot to the plume mode. This would explain the higher electron temperatures measured with the probe, on the order of those expected in a plume mode [START_REF] Csiky | [END_REF]. Nevertheless, plume modes are usually associated with lower electron density. The relatively high electron density measured with the probe, despite a possible transition to the plume regime (where lower electron density is expected), might be attributable to the sputtering of the Langmuir probe, leading to a local increase of plasma density.

Another possible explanation for the discrepancy between the ITS and Langmuir probe temperatures could be the directional sensitivity of the ITS measurements (along with the probed wave vector), which is lacking with the Langmuir probe measurements. Due to the axial electric field a nd low collisionality of the plasma outside the cathode orifice, a nonisotropic electron temperature can be expected with a higher electron temperature along the electric field d irection. The electron temperature obtained from ITS would correspond to the low electron temperature along the radial direction, while the values obtained from the Langmuir probe would correspond to an average over all the directions. An ITS measurement with an axial probed wave vector is needed to settle this question. This measurement was not performed because the available viewing ports of the vacuum chamber did not permit such a scattering configuration.

Finally, as mentioned in section 2.3, the estimation of electron temperature with Langmuir probes implies the assumption of a Maxwell-Boltzmann EVDF. Deviation of the true EVDF from this theoretical distribution may contribute to the discrepancy observed. Although no deviation from the Gaussian shape was observed in Thomson scattering spectra, the contribution from a highly energetic electron population (above 3 eV) might exist, but would not be clearly visible in these Thomson spectra, due to (i) the short spectral range examined for these investigations, in order to facilitate accurate resolution of low electron temperatures, and (ii) the necessity for the highly-energetic electrons to be present at sufficiently high densities for ITS detection. corresponds to the temperature leading to the best agreement between simulations and ITS measurements, as briefly mentioned in the last section. The second profile, obtained with an emissive pellet temperature of 2050 K, corresponds to a temperature leading to a plume mode of the discharge, with large voltage oscillations. As before, the simulation profiles of electron properties are obtained at steady state (≈ 10 ms into the simulation).

Figure 11 shows that there is a good agreement between the axial profiles o f n e a nd T e o btained b y L angmuir probe measurements and from simulations when the emissive pellet temperature is fixed at 2050 K. The discharge mode obtained with a 2050 K emissive pellet temperature is clearly in the plume mode and seems to better simulate the electron properties obtained from Langmuir probe measurements. This observation seems to validate the explanation advanced in the last section for the discrepancy between Langmuir probe and ITS measurements: insertion of the probe inside the plasma switches the discharge from the spot mode to the plume mode. In contrast, the electron temperature obtained from ITS measurements is better simulated with a discharge in a spot mode, with an emissive pellet temperature of 1925 K.

This comparison illustrates once more the sensitivity of agreement to the emitter temperature, which is fixed in the simulation. While the agreement (the order of magnitude of the values obtained in simulation, compared to those experimentally measured) is quite good, further refinements t o t he c ode f or t hermal m odeling c ould likely improve the match between the simulated and measured properties.

Conclusions

In this work, numerical and experimental tools have been applied to the study of electron properties and regime transitions of a 5 A MIREA LaB 6 cathode.

A recently-developed 2D fluid code has been used to simulate the cathode discharge for various discharge currents and emissive pellet temperatures. These results were used to identify the modification of the discharge behavior during the transition from the spot mode to a plume mode. The emissive pellet temperature has been taken as a parameter allowing a scan of varying emission current densities and a means of exploring the physics associated with different regimes.

The recently-developed THETIS incoherent Thomson scattering diagnostic has been used to measure electron properties of this cathode along the cathode centerline for various discharge currents and xenon flow rates. This marks the first time such a diagnostic has been applied for characterizations of this discharge. At higher discharge currents, electron properties suggesting a transition to a plume mode (a non-monotonic spatial variation in temperature, and a slower decrease in electron density) are observed.

Comparisons between the code and reliable ITS measurements show agreements in general trends, however, improvements to the code will be necessary in the future. The selfconsistent coupling of the thermal model for this cathode to the plasma code (as was achieved for a different cathode architecture [31]) will be important for quantitative comparisons.

Complementary Langmuir probe measurements were performed at one discharge condition used for the ITS and numerical investigations. Here, while the ITS and probe measurements were found to yield similar densities, Langmuir probes gave higher electron temperatures. This discrepancy was attributed to a suspected mode transition induced by the invasive nature of the probe, an effect clearly seen from an increase in discharge voltage oscillations during, and even after, the probe residence time in the plasma. This hypothesis was strengthened by comparison with results from simulations, which suggest that the electron property profiles obtained with the Langmuir probe are more consistent with plume mode features. For an apparently simple, unmagnetized plasma, such as the plume of a cathode, invasive diagnostics may have some shortcomings due to plasma perturbation, rendering their use for code validation risky. The successful application of ITS in this work and its comparison to simulations offers a reliable path forward for advancing cathode modeling.

Figure 1 .

 1 Figure 1. Diagram of the MIREA cathode architecture and electrical connections. The cathode is operated in diode mode with an external disc anode. The plasma region is shown in blue.

Figure 2 .

 2 Figure 2. View of the THETIS incoherent Thomson scattering bench setup for cathode investigations (not to scale). Two views are shown: a top view of the bench and laser beam trajectory (in dark green) and a side view, showing the collection of scattered light (in light green) at the top of the vacuum vessel. L0, L1, L2, and L3 indicate lenses; A0, A1, and A2 indicate apertures; VBG indicates the volume Bragg grating notch filter used for stray light suppression.

Figure 3 .

 3 Figure 3. Examples of Thomson spectra obtained on the cathode, illustrating the diagnostic sensitivity. These examples correspond to (a) 2 A discharge current, 6 mm from the cathode orifice, Te = 0.50(5) eV; ne = 2.0(2)× 10 16 m -3 and (b) 16 A discharge current, 1.3 mm from the cathode orifice, Te = 1.2(1) eV; ne = 9.5(8)× 10 17 m -3 . Experimental points, Gaussian fit line and rejected points are shown in blue, green and orange, respectively.

Figure 4 .

 4 Figure 4. Simulation region analyzed with the fluid code to match MIREA architecture. The disc anode is shown in blue, and the disc LaB 6 pellet emitter in red. Note that the domain is cylindrical.

Figure 5 .

 5 Figure 5. Time-averaged 2D maps of the electron density (left) and temperature (right) obtained from the simulation of the cathode discharge steady state with 8 A of discharge current, 0.8 mg s -1 of xenon flow rate and for a electron emitter temperature of 1900 K (top) and 2050 K (bottom). The x-axis and y-axis coordinates are indicated in mm. The color scale for density and temperature has been adapted in order to highlight variations in the electron properties outside the cathode orifice, where the laser scattering measurements were performed.

Figure 6 .

 6 Figure 6.Electron density and temperature obtained by ITS at various distances from the cathode orifice. Measurements were performed for xenon flow rates from 0.2 mg s -1 to 0.8 mg s -1 with a fixed discharge current of 2 A.

Figure 7 .

 7 Figure 7. Electron density and temperature obtained by ITS at various distances from the cathode orifice. Measurements were performed for discharge currents from 2 A to 16 A with a xenon flow rate of 0.8 mg s -1 .

Figure 8 .

 8 Figure 8. Electron density and temperature obtained with a 2 A discharge current and a 0.8 mg s -1 xenon flow rate for an emissive pellet temperature varying from 1750 K to 1850 K.

Figure 9 .

 9 Figure 9. Electron density and temperature obtained with a 4 A discharge current and a 0.8 mg s -1 xenon flow rate for an emissive pellet temperature varying from 1825 K to 1925 K.

Figure 10 .

 10 Figure 10. Electron density and temperature obtained with an 8 A discharge current and a 0.8 mg s -xenon flow rate for an emissive pellet temperature varying from 1900 K to 2000 K.

4. 3 . 2 .

 32 Comparison with fluid s imulations. Results from the Langmuir probe and ITS measurements are also compared to results from simulations in figure11. Two axial profiles from the simulation have been chosen for the comparison. The first profile, obtained with an emissive pellet temperature of 1925 K,

Figure 11 .

 11 Figure 11. Comparison between electron properties obtained from (i) the fluid simulation at steady state (at two different emitter temperatures, in orange and blue solid lines), incoherent Thomson scattering experiments, in red filled circles, and (iii) Langmuir probe measurements, in green filled triangles, for a discharge current of 8 A and flow rate of 0.8 mg s -1 .

Acknowledgment

The authors acknowledge support from CNES and C Boniface and thank K Hara (Stanford University) for discussions.