Progress in modeling and numerical simulation of negative hydrogen ion sources.

Recent progress in numerical modeling of H(-) ion sources has been reviewed. The following two topics are mainly the focus: (1) the origin of the nonuniformity of H(-) production and H(-) beam in large H(-) sources and (2) extraction physics, especially the role of the "weak" transverse magnetic field on H(-) extraction. The present understanding of the origin of H(-) beam nonuniformity is reviewed and discussed with the results by a series of Monte Carlo transport modeling of electrons, neutrals, and H(-) ions. Also, the physical mechanism of the large enhancement of H(-) extraction due to the weak transverse magnetic field is discussed with the particle in cell modeling. In addition, some new results on the effect of surface H(-) production on the electric potential structure near the extraction aperture are presented and discussed.

Analysis of electron energy distribution of an arc-discharge H(-) ion source with Monte Carlo simulation.

For optimization and accurate prediction of the amount of H(-) ion production in negative ion sources, analysis of electron energy distribution function (EEDF) is necessary. We developed a numerical code which analyzes EEDF in the tandem-type arc-discharge source. It is a three-dimensional Monte Carlo simulation code with the effects of cusp, filter, and extraction magnets. Coulomb collision between electrons is treated with Takizuka's model and several inelastic collisions are treated with null-collision method. We applied this code to the JAEA 10 ampere negative ion source. The numerical result shows that the order of electron density is in good agreement with experimental results. In addition, the obtained EEDF is qualitatively in good agreement with experimental results.

Spatial structure of electric potential near the extraction region in Cs-seeded H(-) ion sources.

Structure of electric potential near the extraction region in a negative ion source is investigated analytically with the effect of strong surface H(-) production. The potential profile is analyzed one dimensional by solving the plasma-sheath equation, which gives the electric potential in the plasma region and the sheath region near the wall self-consistently. The potential profile depends on the production rate and the temperature of negative ions. As the production rate becomes large and the negative ion energy becomes small, the potential near the extraction region decreases. The negative potential peak is formed near the plasma grid (PG) surface for the case of large amount and low energy surface production. As a result, negative ions are reflected by this negative potential peak near the PG and returned to the PG surface. This reflection mechanism by the negative potential peak possibly affects the negative ion extraction.

Numerical analysis of H(-) ion transport processes in Cs-seeded negative ion sources.

The H(-) ion transport processes are numerically simulated to understand the extraction process of surface-produced H(-) ions. The three-dimensional transport code using Monte Carlo method has been applied to calculate the H(-) ion extraction probabilities in the model geometry of the JAEA 10 ampere negative ion source. The roles of (1) filter magnetic field and (2) collisions with neutrals (H(0) atoms and H(2) molecules) on the H(-) ion extraction are systematically studied. The results show that H(-) ions are extracted mainly by the filter magnetic field under the low gas pressure condition. The simulation results of extracted H(-) ion beam intensity in the JAEA 10 ampere negative ion source without the magnetic filter tend to be smaller than the experimental results, especially under the low pressure condition. Further model improvements, e.g., modeling and implementation of the effects of the electric field near the extraction aperture, will be required to understand the extraction process of the H(-) ions under the low gas pressure condition.

Effect of energy relaxation of H(0) atoms at the wall on the production profile of H(-) ions in large negative ion sources.

Production and transport processes of the H(0) atoms are numerically simulated using a three-dimensional Monte Carlo transport code. The code is applied to the large JAEA 10 ampere negative ion source under a Cs-seeded condition to obtain a spatial distribution of surface-produced H(-) ions. In this analysis, we focus on the effect of the energy relaxation of the H(0) atoms at the wall on the H(-) ion production from the H(0) atoms. The result indicates that, by considering the energy relaxation of the H(0) atoms at the wall, the production profile of the surface-produced H(-) ion is well reflected in the production profile of the H(0) atom production.

Negative ion production in cesium seeded high electron temperature plasmas.

In experiments on uniformity improvement in a large negative ion source, steep gradients have been observed in the profiles of electron temperature and H(-) ion beam intensity. It has been observed that the gradient in the H(-) ion beam intensity is altered by seeding cesium, though the electron temperature distribution is not affected by Cs. Thus in the Cs seeded condition, the H(-) ion beam intensity is enhanced in local area illuminated by high electron temperature plasmas. A brief analysis suggests possible advantages of high electron temperature plasmas for the negative ion surface production, by enhancement of dissociation to yield proton or atoms as parent particles of the negative ions.

Analysis of electron energy distribution function in the Linac4 H⁻ source.

To understand the Electron Energy Distribution Function (EEDF) in the Radio Frequency Inductively Coupled Plasmas (RF-ICPs) in hydrogen negative ion sources, the detailed analysis of the EEDFs using numerical simulation and the theoretical approach based on Boltzmann equation has been performed. It is shown that the EEDF of RF-ICPs consists of two parts, one is the low energy part which obeys Maxwellian distribution and the other is high energy part deviated from Maxwellian distribution. These simulation results have been confirmed to be reasonable by the analytical approach. The results suggest that it is possible to enhance the dissociation of molecules and the resultant H(-) negative ion production by reducing the gas pressure.

Analysis of the beam halo in negative ion sources by using 3D3V PIC code.

The physical mechanism of the formation of the negative ion beam halo and the heat loads of the multi-stage acceleration grids are investigated with the 3D PIC (particle in cell) simulation. The following physical mechanism of the beam halo formation is verified: The beam core and the halo consist of the negative ions extracted from the center and the periphery of the meniscus, respectively. This difference of negative ion extraction location results in a geometrical aberration. Furthermore, it is shown that the heat loads on the first acceleration grid and the second acceleration grid are quantitatively improved compared with those for the 2D PIC simulation result.

Effect of high energy electrons on H⁻ production and destruction in a high current DC negative ion source for cyclotron.

Recently, a filament driven multi-cusp negative ion source has been developed for proton cyclotrons in medical applications. In this study, numerical modeling of the filament arc-discharge source plasma has been done with kinetic modeling of electrons in the ion source plasmas by the multi-cusp arc-discharge code and zero dimensional rate equations for hydrogen molecules and negative ions. In this paper, main focus is placed on the effects of the arc-discharge power on the electron energy distribution function and the resultant H(-) production. The modelling results reasonably explains the dependence of the H(-) extraction current on the arc-discharge power in the experiments.

Effect of Coulomb collision on the negative ion extraction mechanism in negative ion sources.

To improve the H(-) ion beam optics, it is necessary to understand the energy relaxation process of surface produced H(-) ions in the extraction region of Cs seeded H(-) ion sources. Coulomb collisions of charged particles have been introduced to the 2D3V-PIC (two dimension in real space and three dimension in velocity space particle-in-cell) model for the H(-) extraction by using the binary collision model. Due to Coulomb collision, the lower energy part of the ion energy distribution function of H(-) ions has been greatly increased. The mean kinetic energy of the surface produced H(-) ions has been reduced to 0.65 eV from 1.5 eV. It has been suggested that the beam optics of the extracted H(-) ion beam is strongly affected by the energy relaxation process due to Coulomb collision.

Numerical study of plasma generation process and internal antenna heat loadings in J-PARC RF negative ion source.

A numerical model of plasma transport and electromagnetic field in the J-PARC (Japan Proton Accelerator Research Complex) radio frequency ion source has been developed to understand the relation between antenna coil heat loadings and plasma production/transport processes. From the calculation, the local plasma density increase is observed in the region close to the antenna coil. Electrons are magnetized by the magnetic field line with absolute magnetic flux density 30-120 Gauss which leads to high local ionization rate. The results suggest that modification of magnetic configuration can be made to reduce plasma heat flux onto the antenna.

High current DC negative ion source for cyclotron.

A filament driven multi-cusp negative ion source has been developed for proton cyclotrons in medical applications. In Cs-free operation, continuous H(-) beam of 10 mA and D(-) beam of 3.3 mA were obtained stably at an arc-discharge power of 3 kW and 2.4 kW, respectively. In Cs-seeded operation, H(-) beam current reached 22 mA at a lower arc power of 2.6 kW with less co-extracted electron current. The optimum gas flow rate, which gives the highest H(-) current, was 15 sccm in the Cs-free operation, while it decreased to 4 sccm in the Cs-seeded operation. The relationship between H(-) production and the design/operating parameters has been also investigated by a numerical study with KEIO-MARC code, which gives a reasonable explanation to the experimental results of the H(-) current dependence on the arc power.

Linac4 H⁻ ion sources.

CERN's 160 MeV H(-) linear accelerator (Linac4) is a key constituent of the injector chain upgrade of the Large Hadron Collider that is being installed and commissioned. A cesiated surface ion source prototype is being tested and has delivered a beam intensity of 45 mA within an emittance of 0.3 π ⋅ mm ⋅ mrad. The optimum ratio of the co-extracted electron- to ion-current is below 1 and the best production efficiency, defined as the ratio of the beam current to the 2 MHz RF-power transmitted to the plasma, reached 1.1 mA/kW. The H(-) source prototype and the first tests of the new ion source optics, electron-dump, and front end developed to minimize the beam emittance are presented. A temperature regulated magnetron H(-) source developed by the Brookhaven National Laboratory was built at CERN. The first tests of the magnetron operated at 0.8 Hz repetition rate are described.

Correlation between histopathologic features and magnetic resonance images of spinal cord lesions.

The authors report a correlation between histopathologic features and magnetic resonance images of spinal cord lesions in 19 human spinal cords with various types of lesions. Abnormally high T2-weighted image signal intensities appeared nonspecifically in mildly altered lesions or areas with edema. In the gray matter, a low T1-weighted image in addition to a high T2-weighted image signal intensity appeared in severely altered lesions with necrosis, myelomalacia, or spongiform change. In the white matter, abnormally high T1-weighted image intensities appeared in severely altered lesions. Based on these preliminary results, it appears that a correlation between the degree of histopathologic alteration of the spinal cord and magnetic resonance images in clinical cases can be made. This is the first study to clarify the correlation between histopathologic features and magnetic resonance images of the diseased human spinal cord.

Analysis of plasma distribution near the extraction region in surface produced negative ion sources.

In study of a negative ion source, it is important to understand the plasma characteristics near the extraction region. A recent experiment in the NIFS-R&D ion source has suggested that a "double ion plasma layer" which is a region consisting of hydrogen positive and negative ions exists near the plasma grid (PG). Density distribution of plasma near the extraction region is studied analytically. It is shown that the density distribution depends on an amount of the surface produced negative ions and the double ion plasma layer is formed near the PG surface for the case of strong surface production.

Analysis of rapid increase in the plasma density during the ramp-up phase in a radio frequency negative ion source by large-scale particle simulation.

Numerical simulations become useful for the developing RF-ICP (Radio Frequency Inductively Coupled Plasma) negative ion sources. We are developing and parallelizing a two-dimensional three velocity electromagnetic Particle-In-Cell code. The result shows rapid increase in the electron density during the density ramp-up phase. A radial electric field due to the space charge is produced with increase in the electron density and the electron transport in the radial direction is suppressed. As a result, electrons stay for a long period in the region where the inductive electric field is strong, and this leads efficient electron acceleration and a rapid increasing of the electron density.

Numerical study of the inductive plasma coupling to ramp up the plasma density for the Linac4 H(-) ion source.

In the Linac4 H(-) ion source, the plasma is generated by an RF antenna operated at 2 MHz. In order to investigate the conditions necessary for ramping up the plasma density of the Linac4 H(-) ion source in the low plasma density, a numerical study has been performed for a wide range of parameter space of RF coil current and initial pressure from H2 gas injection. We have employed an Electromagnetic Particle in Cell model, in which the collision processes have been calculated by a Monte Carlo method. The results have shown that the range of initial gas pressure from 2 to 3 Pa is suitable for ramping up plasma density via inductive coupling.

Study of plasma meniscus and beam halo in negative ion sources using three dimension in real space and three dimension in velocity space particle in cell model.

Our previous study by two dimension in real space and three dimension in velocity space-particle in cell model shows that the curvature of the plasma meniscus causes the beam halo in the negative ion sources. The negative ions extracted from the periphery of the meniscus are over-focused in the extractor due to the electrostatic lens effect, and consequently become the beam halo. The purpose of this study is to verify this mechanism with the full 3D model. It is shown that the above mechanism is essentially unchanged even in the 3D model, while the fraction of the beam halo is significantly reduced to 6%. This value reasonably agrees with the experimental result.

Plasma ignition and steady state simulations of the Linac4 H(-) ion source.

The RF heating of the plasma in the Linac4 H(-) ion source has been simulated using a particle-in-cell Monte Carlo collision method. This model is applied to investigate the plasma formation starting from an initial low electron density of 10(12) m(-3) and its stabilization at 10(18) m(-3). The plasma discharge at low electron density is driven by the capacitive coupling with the electric field generated by the antenna, and as the electron density increases the capacitive electric field is shielded by the plasma and induction drives the plasma heating process. Plasma properties such as e(-)/ion densities and energies, sheath formation, and shielding effect are presented and provide insight to the plasma properties of the hydrogen plasma.

Equivalent circuit of radio frequency-plasma with the transformer model.

LINAC4 H(-) source is radio frequency (RF) driven type source. In the RF system, it is required to match the load impedance, which includes H(-) source, to that of final amplifier. We model RF plasma inside the H(-) source as circuit elements using transformer model so that characteristics of the load impedance become calculable. It has been shown that the modeling based on the transformer model works well to predict the resistance and inductance of the plasma.