# Patent

## Title
Charged particle extraction device and method of design there for

## Abstract

The present invention provides a method for extracting a charged particle beam from a charged particle source. A set of electrodes is provided at the output of the source. The potentials applied to the electrodes produce a low-emittance growth beam with substantially zero electric field at the output of the electrodes.

## Background
This application is a Continuation of U.S. application Ser. No. 12/987,396, filed Jan. 10, 2011, which is a Continuation of U.S. application Ser. No. 10/588,981, having a 371 filing date of Mar. 27, 2008, now U.S. Pat. No. 7,872,242, which is the 371 filing of and claims priority from PCT Application No. PCT/US2004/034984, filed on Oct. 16, 2004, which claims priority from Australian Application No. 2003905709, filed Oct. 17, 2003, which are hereby incorporated by reference.


TECHNICAL FIELD OF THE INVENTION
The present invention relates to the extraction of a charged particle beam from a plasma. The invention is particularly directed toward the extraction of high brightness beams.
BACKGROUND OF THE INVENTION
There exists a range of applications for ion beams particularly in the semiconductor industry. For example, the fabrication and correction of lithography masks involves sub-micron etching capability. This is currently achieved using medium-energy particle beams (10-50 kilo-electron-Volts (keV)), commonly referred to as Focused Ion Beams (FIB). To enable sub-micron feature creation, the FIB must be capable of focus down to a nanometer scale spot size. This requires the extraction of very high brightness beams in excess of 105 Angstroms per steradian per meter squared (A/sr/m2).
Liquid Metal Ion Source (LMIS) technology has been capable of this level of brightness for many years. The technology exploits the capillary effect of liquid Gallium to cover a sharp Tungsten needle onto which a strong electric field is applied, thereby removing ions. The effect of the field is strongest at the needle point and so a beam of ions is created that appears to diverge from a nanometer spot. The beam is then accelerated and focused onto the target where it sputters the surface by collision processes.
Though LMIS technology may present nanometer scale milling capability, it produces unwanted doping effects by introducing Gallium into a substrate or target. To avoid this, a high brightness beam of inert ions would be preferable. Inert ions could be extracted from an inert ion gas plasma. But this has proved difficult and much research is devoted to improvement of extraction mechanisms to extract the ions from the plasma in the form of a beam. For example, electrode extraction optimization by adjustment of aperture ratios and electrode spacing is described by J. R. Copeland, et. al., “A study of the ion beam intensity and divergence obtained from a single aperture three electrode extraction system”, Rev. Sci. Instrum., 44(9):1258, 1973. Other references describing shaping of electrodes include D. E. Radley, “The theory of the pierce type electric gun” J. Electron. Control, 4:125, 1957. E. R. Harrison, “Approximate electrode shapes for a cylindrical electron beam” Brit. J. Appl. Phys., 5:40, 1953, and P. N. Daykin, “Electrode shapes for a cylindrical electron beam” Brit. J. Appl. Phys., 6:248, 1955. Despite these efforts brightness in excess of 105 A/sr/m2 has not been achieved with a plasma ion source. Thus, in the field of extracting an ion beam from a plasma, design of electrodes for extracting a beam of high brightness is desired.
SUMMARY OF THE INVENTION
The present invention provides a method for design of an extraction device to achieve a charged particle beam of high brightness. According to an aspect of the present invention a set of electrodes, each with an aperture, is provided and a beam is drawn from a charged particle source through the apertures of the electrodes. The shapes of the electrodes and the potentials applied to them produce a low-emittance beam.
According to another aspect of the invention, the shapes and potentials applied to the electrodes produce a substantially zero electric field in the vicinity of the aperture of the last electrode furthest from the particle source. The shapes, potentials and positions of the electrodes to achieve zero or low emittance growth are determined from a set of boundary conditions applied at concentric surfaces. The boundary conditions include a substantial, non-zero, electric field in the vicinity of the aperture of the first electrode.
The foregoing has outlined rather broadly aspects, features and technical advantages of the present invention in order that the detailed description of the invention that follows may be better understood. Additional aspects, features and advantages of the invention will be described hereinafter. It should be appreciated by those skilled in the art that the disclosure provided herein may be readily utilized as a basis for modifying or designing other structures for carrying out the same purposes of the present invention. Persons of skill in the art will realize that such equivalent constructions do not depart from the spirit and scope of the invention as set forth in the appended claims, and that not all objects attainable by the present invention need be attained in each and every embodiment that falls within the scope of the appended claims.



BRIEF DESCRIPTION OF THE DRAWINGS
For a more complete understanding of the present invention, and the advantages thereof, reference is now made to the following descriptions taken in conjunction with the accompanying drawings, in which:
FIG. 1 is a block diagram of an embodiment of the invention.
FIG. 2 is an illustration of electrodes in an embodiment of the invention.
FIG. 3 shows a potential distribution for an embodiment of the invention
FIG. 4 shows concentric surfaces to which boundary conditions are applied according to the method of the present invention.
FIG. 5 is a graph showing the effects of non-zero initial gradient
FIG. 6 is a contour of integration of a complex integral.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
The following is a detailed description of example embodiments of the invention depicted in the accompanying drawings. The example embodiments are in such detail as to clearly communicate the invention. However, the amount of detail offered is not intended to limit the anticipated variations of embodiments but, on the contrary, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the present invention as defined by the appended claims. The detailed descriptions below make such embodiments obvious to a person of ordinary skill in the art.
The present invention is particularly directed to the extraction of a high brightness ion beam with low aberrations, and is applicable to the formation of low-emittance charged particle beams. Generally, minimizing emittance ensures the highest possible beam brightness. Emittance is a measure of the parallelism of the individual particle trajectories in a beam. For planar symmetry, particles follow rectilinear and parallel trajectories. For cylindrical and spherical symmetry they move along rectilinear paths that follow radial lines as though diverging from or converging to a single line or point.
Thus, the design of low aberration electrodes starts with the accurate description of the beam to be extracted. Assumed throughout this description is that the desired beam density profile is radially uniform (with respect to the beam axis of symmetry). However, the end user of the device must specify the desired beam shape, beam current density and extraction potential.
The beam shapes described herein are categorized by the aspect of the aperture through which they are extracted, being either a rectangular slit or a circular orifice, and by the angle of divergence of the beam envelope. With these categorizations the following beam types are possible:
1) Parallel beam envelope: cylindrical (circular aperture) and strip (rectangular slit) beams
2) Diverging beam envelope: diverging conical beam (circular aperture) and diverging wedge beams (rectangular slits)
3) Converging beam envelope: converging conical beam (circular aperture) and converging wedge beams (rectangular slits).
These beam types can be thought of as angular sections of current flowing between the concentric surfaces of a standard classical diode as described by Child, Langmuir and Blodgett (CLB). There are three diode configurations. The first is two parallel plates. The second is two concentric cylinders and the third is two concentric spheres. All surfaces are assumed to be perfect conductors. In the case of cylindrical and spherical diodes, current flows from the inner surface to the outer surface to yield a diverging current profile. Or, current flows from the outer surface to the inner surface to produce a converging current profile. If the conducting surfaces are parallel, the current can flow in either direction and results in a parallel current profile.
In the basic CLB analysis the emitting surface is assumed to have an undepleteable source of charged particles with no inherent thermal energy and zero electric field. The emitted particles are taken to enter the extraction gap with no initial velocity. Under these circumstances a current is caused to flow when a potential is applied between the two conducting surfaces (also called electrodes). This potential drop is referred to as the extraction potential.
Because of symmetry considerations, the cylindrical and spherical diodes can be analyzed in two dimensions without any loss of generality. A polar reference frame can be defined with its origin, O, at the center of concentricity and the vector r running along any radius emanating from O. The potential distribution along any radius r is identical from symmetry considerations and is determined by both the applied extraction potential and the presence of charged particles flowing as current between the emitting and collecting surfaces. The function that describes this potential distribution (and consequently the density distribution) of the inter-electrode gap along any radius r, subject to the above mentioned assumptions was first described in the parallel case by Child and Langmuir (independently) and in the cylindrical and spherical cases by Langmuir and Blodgett (together).
It was Pierce in the 1940s who first suggested that electrodes could be designed to form a beam that flowed like the particles in a classical diode. Because of symmetry considerations, azimuthal forces in classical diodes can be ignored because adjacent particles apply equal and opposite forces on each other and hence only radial forces exist. Pierce reasoned that an electric field structure could be set-up in the charge free region adjacent to the beam that compensated azimuthal space charge forces in the same way that charged particles do in diodes. He solved this problem for a parallel strip beam. In 1964, Radley solved Laplace's equation (the relation describing the electric field in the charge free region) in general terms, repeating Pierce's finding in a more rigorous fashion and extending the solutions to encompass a broad range of beam types including those mentioned above.
In principle the beams flowing in the Pierce and Radley electrodes are aberration free, but they only describe current flow between two electrodes. The problem addressed here is that the electric field in the beam at the second electrode as described by the CLB potential profiles is very large. Well known from field theory is that the presence of an electric field on one side of an aperture causes the electric field to balloon out to the field free region, resulting in equipotential lines that are strongly curved in the neighborhood of the aperture. This creates a lens effect that disturbs the beam density profiles (and hence the potential profiles). In this case the CLB equations can no longer be used. The electrode structure no longer adequately compensates the space charge forces in the beam. This results in strong emittance growth and a decline in brightness (a key metric of beam quality).
The problem is solved with multiple stages of electrodes. Thus, the analysis uses more than two concentric surfaces so that there exist two extraction stages. The first stage results in a strong electric field. The second stage serves to bring the strong electric field at the end of the first stage back to zero so that no lens effect occurs at the end of either the first or second stage. The beam potential profile in the second stage is obtained by generalizing the CLB equations. The potential profile can be used to determine the ideal electrode shapes for the second stage by substitution into Radley's equations.
Most ion extraction optics consist of at least three electrodes to block the passage of electrons (created by secondary emission) from the transport region aft of the extractor. The first two serve to extract, focus and accelerate the ions and the third to create a potential barrier to the passage of electrons. Typically this is achieved by applying the extraction potential to the first electrode (in contact with the plasma) and extracting the beam initially to some negative potential applied to the second electrode and then returning the beam potential to ground between the second and third electrodes. The net result is that the beam particles have a final energy equal to the extraction potential. But the electrons in the transport region see a potential barrier equal to the potential difference between the second and third electrodes and hence are blocked from entering the extraction region where their presence has nefarious effects on the beam distribution profiles.
In this simple, common triode system, the electric field in the beam is undesirably still very strong at both the ends of the extraction stage (between electrodes 1 and 2) and the blocking stage (between electrodes 2 and 3). Thus, in an embodiment of the present invention, the electric field at the output of the exit aperture of a three electrode system is brought to zero. This embodiment is sufficient to produce a low, ideally zero emittance growth beam in the absence of electrons.
In another embodiment, a five electrode system is employed. In the five electrode system, the electric field in the beam is brought to zero at (in the vicinity of) the apertures of the third electrode and the fifth electrode. Electrodes one, two and three form a two-stage extraction region in which the electric field magnitude rises from its value at the meniscus to a high value at the aperture of the second electrode and then declines to substantially zero at the aperture of the third electrode. Electrodes three, four and five form a two-stage blocking region in which the electric field magnitude rises from substantially zero at the aperture of the third electrode to an intermediate value at the aperture of the fourth electrode and then declines again to zero at the aperture of the fifth electrode. In this embodiment, the blocking stage inhibits the presence of electrons while maintaining the desired beam profile. The result is a 5 electrode extraction system. This is the number of electrodes used to ensure the electric field is brought to zero at the end of both the extraction and blocking stages while taking into account the presence of electrons. If electrons were not present only three electrodes would be required.
In the analysis presented herein, the beam distributions are generalized to take into account the non-zero electric field at the meniscus and the non-zero velocity at the entry to the blocking region. We use Radley's solutions to determine the electrode shapes. Except for the first and last electrodes, both faces of each electrode are shaped. For the first electrode only the face adjacent to the beam needs to be shaped. For the fifth electrode only the front face (relative to the direction of the current flow) needs to be shaped.
A major assumption of the CLB equations is that the electric field at the emitting surface is zero. This is not true in the case of plasmas where the emitting surface is the plasma meniscus (interface between the plasma and the beam) at which there exists a strong electric field typically on the order of the hundreds of kilo-Volts but also in some cases in the Mega-Volts. This must be taken into account in the description of the beam potential profile in the first stage of the extraction region. Thus, the basic CLB equations cannot be applied to beams extracted from plasmas under the assumptions presented above. To apply them to this analytical technique for the design of low emittance electrodes, they must be generalized.
For the second stage of both the extraction and blocking regions the generalization is simpler requiring only a change in boundary conditions to obtain the CLB solutions in reverse. It is this that allows the potential profile to be brought to zero. Note also, that after the first extraction stage, initial velocity must also be taken into account, but this is done with a simple mathematical transformation.
FIG. 1 shows a focused ion beam system 8 that includes an evacuated envelope 12 in which is located a plasma source and an extraction mechanism 14, to provide a dense plasma for ion beam focusing column 16. An ion beam 18 passes from source 14 through column optics 16 and between electrostatic deflection mechanism 20 toward specimen 22, which comprises, for example, a semiconductor device positioned on movable X-Y stage 24 within lower chamber 26. An ion pump 28, operating in conjunction with pumping system 30 and vacuum control 32, is employed for evacuating the source and maintaining high vacuum in the upper column optics region. The vacuum system provides within chamber 26 a vacuum of typically between approximately 1×10-7 Torr and 5×10-4 Torr, with nominally 10 mTorr in the plasma source and <1×10-6 Torr in the column optics chamber.
High voltage power supply 34 is connected to ion source 14 as well as to appropriate electrodes in focusing column 16 for forming an approximately ion beam 18 and directing the same downward. Deflection controller and amplifier 36, operated in accordance with a prescribed pattern provided by pattern generator 38, is coupled to deflection plates 20 whereby beam 18 may be controlled to trace out a corresponding pattern on the upper surface of specimen 22. In some systems, the deflection plates are placed before the final lens, as is well known in the art.
The ion beam source 14 is brought to a focus at specimen 22 for either modifying the surface 22 by ion milling, material deposition, or for the purpose of imaging the surface 22. A charged particle multiplier 40 used for detecting secondary ion or electron emission for imaging is connected to video circuit and amplifier 42, the latter supplying drive for video monitor 44 also receiving deflection signals from controller 36. The location of charged particle multiplier 40 within chamber 26 can vary in different embodiments. For example, a preferred charged particle multiplier 40 can be coaxial with the ion beam and include a hole for allowing the ion beam to pass. A scanning electron microscope 41, along with its power supply and controls 45, are optionally provided with the FIB system 8.
Signals applied to deflection controller and amplifier 36, cause the focused ion beam to move within a target area to be imaged or milled according to a pattern controlled by pattern generator 38. Emissions from each sample point are collected by charged particle multiplier 40 to create an image that is displayed on video monitor 44 by way of video circuit 42. An operator viewing the image may adjust the voltages applied to various optical elements in column 16 to focus the beam and adjust the beam for various aberrations. Focusing optics in column 16 may comprise mechanisms known in the art for focusing or methods to be developed in the future.
An embodiment of an extraction mechanism for extracting an ion beam from the plasma source is illustrated in FIG. 2, which shows a sequence of electrodes, each with an aperture corresponding to a beam profile. Electric potentials are applied to the electrodes to compensate for azimuthal space charge forces and to ensure that the electric field at the output of the sequence of electrodes is brought to zero.
The potentials shown in FIG. 2 are nominal values for a diverging beam from a Krypton plasma density of about 1013 cm−3. The actual potentials and shapes of the electrodes will depend upon the selected beam type, the desired extraction energy, the total beam current and current density, as well as electrode potentials selected by the designer. As will be seen, given the methodology disclosed herein, zero or nearly zero emittance growth can be achieved with a plurality of combinations of voltages and corresponding shapes. Thus, the potentials shown in FIG. 2 are for exposition.
FIG. 2 shows a first electrode 200 closest to the plasma source, referred to herein as the plasma electrode. The plasma electrode exhibits a high electric potential of 27.4 kilo-Volts (kV) with a substantial non-zero electric field on the order of 100 kV/meter in the vicinity of its aperture. A second electrode 202, referred to herein as the acceleration electrode, exhibits a potential of about 17.3 kV. The electric field in the vicinity of the aperture of the second electrode is on the order of 100 MV/m (Mega-Volts per meter).
A third electrode, referred to herein as the first blocking electrode, exhibits a potential of −200V, with a substantially zero electric field in the vicinity of the aperture of the first blocking electrode. A fourth electrode, referred to herein as the second blocking electrode, exhibits a potential of −100V with a substantial, but intermediate, value of electric field in the vicinity of its aperture. A fifth electrode, referred to herein as the ground electrode, exhibits a potential of zero volts and exhibits a substantially zero electric field in the vicinity of its aperture.
Stated in more general terms, the embodiment of FIG. 2 show that V1, the potential on electrode 200 exceeds V2, the potential on electrode 202. V2 exceeds V3, the potential on electrode 204. Also, V5, the potential on electrode 208 exceeds V4, the potential on electrode 206. V4 exceeds V3. V3 and V4 are negative and V5 is zero.
Electrodes 200 and 202 form a first extraction stage 211. Electrodes 202 and 204 form a second extraction stage 213. Electrodes 204 and 206 form a first blocking stage 215. Electrodes 206 and 208 form a second blocking stage 217. Note that the blocking stage formed by electrodes 204, 206 and 208 impedes the flow of electrons into the extractor region. In the absence of electrons, the electrodes 206 and 208 would not be necessary.
The potentials applied to the electrodes cause a beam flow that produces a substantially zero electric field in the vicinity of the aperture of the ground electrode 208 and in the vicinity of the aperture of the first blocking electrode 204. Bringing the electric field to zero at the end of the extraction stage and at the end of the blocking stage is desired to stop the beam trajectories from being perturbed in the vicinities of these apertures.
The potentials of the electrodes cause a beam flow that produces a high potential gradient in the vicinity of the aperture of the acceleration electrode; a low, preferably zero, potential gradient in the vicinity of the first blocking electrode; an intermediate potential gradient in the vicinity of the aperture of the second blocking electrode; and a low, preferably zero, potential gradient in the vicinity of the aperture of the ground electrode.
The shapes and potentials applied to the electrodes are determined from a set of boundary conditions applied at concentric equipotential surfaces 201, 203, 205, 207, and 209. Each concentric surface is at a radius corresponding to a position of an electrode. That is, the positions of the electrodes are at the positions of the concentric surfaces, and the position of the concentric surfaces are determined from the boundary conditions. Initially, the designer specifies beam type, and the extraction energy, and hence the voltage of the plasma electrode 200. The designer also specifies the beam current density. This defines the required plasma density, and hence, the electric field at the meniscus. The boundary conditions, as described herein, are applied to each concentric surface to produce a low or zero emittance growth beam.
In preferred forms of the invention the beam profile is one of a strip beam, a wedge beam, a cylindrical beam, and a conical beam. These beams can be thought of as sections of current flowing in complete diodes. Consequently, parallel beams can be treated in the same way as a planar diode. Divergent/convergent beams extracted through a rectangular slit can be treated as a cylindrical diode. Divergent/convergent beams extracted through a circular aperture can be treated as a spherical diode. In each case, the plasma/beam interface or plasma meniscus is the anode (emitter) and the 0 Volt equipotential surface is the cathode (collector). The purpose of the extraction device is to ensure that the meniscus and 0 Volt equipotential surfaces are parallel in the case of parallel beams and concentric cylinders or spheres in the case of diverging or converging beams. If this situation is maintained then, neglecting the inherent ion temperature in the plasma, the ions travel in perfect parallelism along radial lines from the meniscus to the 0 Volt equipotential surface and suffer no deflection. In other words the beams have zero emittance growth.
By accelerating the beam, a significant electric field can be produced in the direction of flow. In accordance with the present invention this is compensated to avoid the electric field in the beam from ballooning outwards at the exit of the extractor and deflecting the ion trajectories. In the transport region aft of the extractor, considerable advantage can be gained by the presence of electrons which neutralize the beam space charge. However, their presence in the extractor, especially in the acceleration gap, is highly detrimental to beam quality since they alter the charge density distribution and hence the potential structure through-out the extractor. The existence of an electric field at the exit to the extractor would serve to accelerate electrons into the extraction region in such quantity as to neutralize the electric field at the extractor exit. Accordingly the present invention ensures that the beam potential has zero electric field upon exit of the extractor.
The electrons in the beam plasma aft of the extractor have non-zero temperature and are generally distributed according to Maxwell's law. This means that higher energy electrons from the tail of the distribution are able to enter the extractor if a blocking field is not present. To this end a small blocking potential preferably of some hundred Volts is produced at the end of the extraction region to inhibit the passage of electrons. Again, to avoid ballooning of the electric field this potential must be produced so that the electric field upon exit of the blocking region is zero.
In a preferred form of the invention the potential distribution in the beam will then have the form shown in FIG. 3. There are two regions called the Extraction and Blocking regions both of which are divided into two stages. Stage 1, 301, of the extraction region takes the beam from a low gradient to a high gradient and stage 2, 303 takes the beam from the high gradient to a low gradient. In the extraction region, stage 1 is necessary to match the beam potential to the plasma sheath and stage 2 is necessary to bring the electric field to zero. A similar rationale applies to the two stages of the blocking region.
The extractor region is formed by the first three sequential electrodes. Stage 1 and stage 2 of the extractor region are separated by the acceleration electrode. The blocking region is formed between the first blocking electrode and the last electrode. Stage 1, 305, and stage 2, 307, of the blocking region are separated by the second blocking electrode. In an embodiment of the present invention the last electrode is a “ground” electrode with a zero electric field and zero potential in the beam in the vicinity of the aperture of the last electrode.
In order to describe the plasma extractor device of the present invention it is desirable to first provide some analysis of beam distributions and an overview of solutions to Laplace's equation for different beam profiles. This will be done in separately headed sections which precede an explanation of electrode design and specific examples of the invention.
FIG. 3 shows the general shape of the potential along the beam inside the extractor region and blocking region for low, ideally zero, emittance growth extraction. The overall potential profile is in fact four distributions concatenated. Each stage. 301, 303, 305, and 307, is bounded by two plane or two concentric surfaces depending on the beam type, so that each stage can be treated as separate, ‘complete’ diodes, each with its own set of boundary conditions. To get the complete distribution, the solutions are stitched together by matching the boundary conditions at each surface.
Thus, referring to FIG. 4, the first extraction stage 401 is bounded by surface 401 and surface 402. The second extraction stage 403 is bounded by surface 402 and surface 404. The first blocking stage 405 is bounded by surfaces 404 and 406. The second blocking stage 407 is bounded by surfaces 406 and 408. The positions of each surface correspond to the positions of the respective electrodes.
Referring again to FIG. 3, at point 300, which corresponds to the surface in the proximity of the aperture of the first electrode 200, the electric potential is a high value. Also, the potential gradient (electric field), which is the slope of the potential curve, at point 300 is a substantial non-zero value. The present invention takes into account this non-zero electric field. At point 302, which corresponds to the surface in the proximity of the aperture of acceleration electrode 202, the slope (potential gradient) is a large value. Thus, the electric field magnitude is high in the vicinity of the aperture of the acceleration electrode.
At point 304, which corresponds to the surface in the proximity of the aperture of first blocking electrode 204, the electric potential is a relatively small, but non-zero, negative value, and the electric field magnitude is substantially zero. At point 306, which corresponds to the surface in the proximity of the aperture of the second blocking electrode 206, the electric potential is an intermediate, non-zero, value. The electric field magnitude there is a non-zero intermediate value. At point 308, which corresponds to the surface in the proximity of the aperture of the last electrode 208, the potential and potential gradient are both zero.
The potential and charge distributions in the beam at each stage of extraction are governed by Poisson's law. The basic solution to this problem, using simple boundary conditions, in plane symmetry was solved by Child and Langmuir in 1911 and 1914 respectively, and in cylindrical and spherical symmetry by Langmuir and Blodgett in the 1920s. These initial solutions assumed only one charged particle species (notably electrons), and ignored initial velocity.
Applicants introduce several generalizations. In particular, we consider distributions describing non-zero initial gradient, distributions describing non-negligible initial velocity, distributions tapering from a strong gradient to a zero gradient (the reverse of the standard Langmuir-Blodgett solution) and distributions describing the presence of Maxwellian electrons. A further requirement on these distributions is that they have the same form as the standard Langmuir-Blodgett series solutions since the Radley solutions to Laplace's equation are dependent on this form.
The work published by Langmuir-Blodgett in the 1920s forms the basis for the beam distribution analysis that follows and is given in terms of simple and compact series solutions. A cursory mathematical overview of Langmuir-Blodgett's contribution starting with spherical symmetry and working through cylindrical to planar symmetry is provided. In the case of ion beam extraction from a circular aperture the charge and potential distributions in the beam are assumed to be analogous to those in a complete spherical diode.
Following Langmuir and Blodgett, Poisson's equation between two concentric spheres can be stated as




  
    
      
        
          
            1
            
              r
              2
            
          
          ⁢
          
            
              d
              
                
              
            
            
              d
              r
            
          
          ⁢
          
            (
            
              
                r
                2
              
              ⁢
              
                
                  d
                  V
                
                
                  d
                  r
                
              
            
            )
          
        
        =
        
          ρ
          
            ɛ
            o
          
        
      
    
    
      
        (
        1
        )
      
    
  




where V is the potential at a point a distance r from the common center and ρ is the ion charge density. The current flowing in the diode can be written in terms of the particle velocity, v:

I=4πr2ρv  (2)

where the velocity can be written in terms of the voltage V using the kinetic energy relation:





  
    
      
        
          
            1
            2
          
          ⁢
          
            Mv
            2
          
        
        =
        
          -
          eV
        
      
    
    
      
        (
        3
        )
      
    
  




Combining equations (1), (2) and (3) yields:




  
    
      
        
          
            
              r
              2
            
            ⁢
            
              
                
                  d
                  2
                
                ⁢
                V
              
              
                d
                
                  r
                  2
                
              
            
          
          +
          
            2
            ⁢
            r
            ⁢
            
              
                d
                V
              
              
                d
                r
              
            
          
        
        =
        
          
            A
            ⁡
            
              (
              
                -
                V
              
              )
            
          
          
            1
            /
            2
          
        
      
    
    
      
        (
        4
        )
      
    
  




where:





  
    
      
        A
        =
        
          
            I
            
              4
              ⁢
              
                
              
              ⁢
              
                πɛ
                o
              
            
          
          ⁢
          
            
              M
              
                2
                ⁢
                e
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




Equation (4) can probably not be integrated directly but a series solution can be found. The form of the solution is a function of the ratio R=r/rs




  
    
      
        
          V
          ⁡
          
            (
            R
            )
          
        
        =
        
          
            
              (
              
                
                  9
                  4
                
                ⁢
                A
              
              )
            
            
              2
              3
            
          
          ⁢
          
            
              f
              
                4
                3
              
            
            ⁡
            
              (
              R
              )
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




where f is the analytic function to be found. The term





  
    (
    
      
        9
        4
      
      ⁢
      A
    
    )
  
  
    2
    3
  




serves to normalize for the constant term A, which is related to the current and hence the plasma density and meniscus curvature. The term





  
    f
    
      4
      3
    
  
  ⁡
  
    (
    R
    )
  




serves to remove the square root and hence to simplify subsequent derivations.

A further transformation is performed by setting:

γ=ln(R)  (7)

so that a solution to equation 6 can be expressed in terms of a MacLauren series as follows:





  
    
      
        f
        =
        
          
            ∑
            
              n
              =
              0
            
            ∞
          
          ⁢
          
            
          
          ⁢
          
            
              a
              n
            
            ⁢
            
              y
              n
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




Now substituting equation (6) into equation (4) and using equation (7) results in:

3ff′+f′2+3ff′′=1  (9)

where the prime denotes the derivative with respect to γ and the double-prime indicates the second derivative with respect to γ. From the form of equation (9) it can be seen that where f=0, we have f′=1. Then, the first six terms of the series solution are:

f=γ−0.3γ2+0.075γ3−0.0143182γ4+0.0021609γ5−0.00026791γ6  (10)

corresponding to V=0 and V′=0.

A similar derivation is made for the case of cylindrical symmetry. Poisson's equation becomes:




  
    
      
        
          
            r
            ⁢
            
              
                
                  d
                  2
                
                ⁢
                V
              
              
                d
                
                  r
                  2
                
              
            
          
          +
          
            
              d
              V
            
            
              d
              r
            
          
        
        =
        
          
            B
            ⁡
            
              (
              
                -
                V
              
              )
            
          
          
            1
            /
            2
          
        
      
    
    
      
        (
        11
        )
      
    
  




where:





  
    
      
        B
        =
        
          I
          
            2
            ⁢
            π
            ⁢
            
              
            
            ⁢
            l
            ⁢
            
              
                m
                
                  2
                  ⁢
                  e
                
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




and 1 is the length of the extraction slit. The solution takes the form:





  
    
      
        
          V
          ⁡
          
            (
            R
            )
          
        
        =
        
          
            
              (
              
                
                  9
                  2
                
                ⁢
                Br
              
              )
            
            
              2
              3
            
          
          ⁢
          
            g
            
              4
              3
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




where g is the analytic function to be found over the desired range of r. Substituting equation (13) into equation (11) and using equation (7) results in:

3gg′′+g′2+4gg′+g2=1  (14)

For g=0, g=1 and the series solution is:

g=γ−0.4γ2+0.0916667γ3−0.01424242γ4+0.001679275γ5−0.0001612219γ6  (15)

corresponding to V=0 and V′=0.

Non-Zero Initial Gradient

The assumption that the electric field at the meniscus surface is zero, central to the Child-Langmuir and Langmuir-Blodgett derivations, is incorrect and so will be herein generalized. Furthermore, the final result will be presented in an easy to use series formulation similar to the original Langmuir-Blodgett relation. FIG. 5 demonstrates the effect of non-zero initial gradient. The dashed line corresponds to the case of a divergent conical beam extracted from a plasma of density n=1014 cm−3 and is compared to the standard Langmuir-Blodgett distribution which assumes zero initial gradient.
To understand why the gradient at the plasma/sheath interface and meniscus is non-zero it is necessary to consider the Bohm sheath criterion which stipulates the minimum ion velocity for entry into the sheath to maintain a stable sheath at a plasma boundary. In conjunction with some distribution relation for electrons, this defines a potential structure within the sheath. In particular, the electric field at the plasma boundary is non-zero and is typically several hundred kilo-volts per meter. For continuity of the electric field across the plasma/beam interface the potential gradient must be equal on both sides of this surface, This is not assured by the original assumptions of Langmuir and Blodgett who were modeling particle flow from thermionic cathodes. In that case the source of particles was assumed to be undeletable and to have no intrinsic electric field so that for equilibrium the boundary condition in the extractor was for the electric field to be zero at the entry to the acceleration gap. In the case of plasmas, the flux of ions is fixed and an intrinsic electric field does exist in the sheath which separates the beam from the bulk plasma. This strongly implies that particle beam extraction from plasmas is not space charge limited but rather is source limited and further implies that the voltage distribution in the beam is not given by the Langmuir-Blodgett relation as is usually stated.
General solutions to Poisson's equation for the case of non-zero initial gradient are first present followed by a discussion of the plasma sheath and how the gradient at the meniscus is obtained.
Spherical Symmetry
Taking the first derivative of equation (6), yields:

V′(γ)=λ√{square root over (p)}  (16)

where





  λ
  =
  
    
      4
      3
    
    ⁢
    
      
        (
        
          
            9
            4
          
          ⁢
          A
        
        )
      
      
        2
        3
      
    
  




and ρ=f′(γ)2f(γ)2/3, for γ=0. In the classic Langmuir-Blodgett derivation, f(γ)=0 in (9) leads to f′(γ)=0, for γ=1, (assuming potential increases as a function of position in the extractor,) so that p=O and hence V′(γ=0)=0. In other words, the Langmuir-Blodgett derivation requires a zero initial gradient in potential. Numerically, however, it is possible to have f approach zero without requiring V′(0) to be zero by setting f′(0) such that (16) holds for the desired value of V′(0). For a given value of V′(0) and A, there is a limit to how small f(0) can be set, but in most practical cases it is several orders of magnitude less than unity.

The solution depends only on p, rather than the individual values of f(0) and f′(0), for the range of p which is of interest. Let the series coefficients, an, be expressed as a quadratic:

an=αn+βn+γn  (17)

where the α, β, and γ are the expansion coefficients found by a least squares method. These terms are presented in the following table.








TABLE 1







Expansion terms for the coefficients of the MacLauren series (see


equation 17) for the spherical case. Note that the αn are very


close to the original Langmuir-Blodgett series coefficients.












n
αn
βn
γn
















1
1.0035
4.049
−10.92



2
−0.3084
−8.008
25.11



3
0.08338
7.791
−25.85



4
−0.01825
−3.96
13.47



5
0.002870
1.0004
−3.448



6
−0.0002227
−0.09904
0.3441










Cylindrical Symmetry

Taking the first derivative of the potential and writing it in terms of the parameter p, we have:




  
    
      
        
          
            V
            ′
          
          ⁡
          
            (
            R
            )
          
        
        =
        
          μ
          ⁡
          
            (
            
              
                g
                
                  4
                  3
                
              
              +
              
                2
                ⁢
                
                  p
                
              
            
            )
          
        
      
    
    
      
        (
        18
        )
      
    
  




where





  μ
  =
  
    
      2
      3
    
    ⁢
    
      
        (
        
          
            9
            2
          
          ⁢
          B
          ⁢
          
            
          
          ⁢
          
            e
            γ
          
        
        )
      
      
        2
        3
      
    
  




and with R=1. Note that in the limit as g approaches zero the term in g disappears, so that V′(R)≈2μ√{square root over (ρ)}. This has the same form as the spherical case discussed above. We write:





  
    
      
        g
        =
        
          
            ∑
            
              n
              =
              0
            
            ∞
          
          ⁢
          
            
          
          ⁢
          
            
              b
              n
            
            ⁢
            
              γ
              n
            
          
        
      
    
    
      
        (
        19
        )
      
    
  




and plotting the series coefficients in terms of the parameter p yields:

bn=αn+βn+γn  (20)

where the α, β, and γ are the expansion coefficients found in the same way as for the spherical case. These terms are presented in the following table.








TABLE 2







Expansion terms for the coefficients of the MacLauren series (see


equation 20) for the cylindrical case. Note again that the αn are


very close to the original Langmuir-Blodgett series coefficients.












n
αn
βn
γn
















1
1.0034
3.989
−10.69



2
−0.4086
−8.223
25.43



3
0.1005
7.974
−26.22



4
−0.01866
−4.038
13.65



5
0.002658
1.0184
−3.491



6
−0.0002011
−0.1007
0.3482










The Plasma Sheath

There are several ways to model the sheath. A method in terms of the Bohm sheath criterion and the Boltzman equation in one dimension, which is easily extended to spherical and cylindrical symmetry is described in full in M. A. Lieberman and A. J. Lichtenberg, “Principles of Plasma Discharges and Materials Processing” John Wiley & Sons, New York, 1st edition, 1994. This model incorporates both the non-zero ion velocity vB at the entry to the sheath, required by the Bohm criterion, and also the presence of electrons. Another popular method is the Child sheath, which is extended to spherical and cylindrical symmetry by the use of the standard Langmuir-Blodgett corrections. In this case, the pre-sheath/sheath boundary and the meniscus are considered to be concentric spheres for extraction from a circular aperture and concentric cylinders for extraction from a slit. Though solving the Boltzman sheath is possible in terms of a series, it requires a somewhat more drawn out analysis and so for simplicity, the Child sheath method will be employed with Langmuir-Blodgett corrections. As such, the pre-sheath can be ignored and we assume that the velocity of ions and the potential at the bulk plasma/sheath edge are zero. However, it should be noted that the Child sheath yields smaller gradients than the Boltzmann sheath.
Sheath Potential at the Meniscus
To solve equation (4), the case of extraction from a plasma requires that I must equal the ambipolar flux for ions

Γ=0.6enevBA  (21)

where ne is the plasma density at the sheath edge and A=4πrs2 is the area over which current is extracted, with rs being the radius of curvature of the sheath. It follows from equation (5) and equation (21) that the solution to equation (4) is strongly related to both ne and rs.

In addition, three boundary conditions are imposed. At the entry to the sheath set




  
    V
    ⁡
    
      (
      1
      )
    
  
  =
  
    
      0
      ⁢
      
        
      
      ⁢
      and
      ⁢
      
        
      
      ⁢
      
        
          d
          
            V
            ⁡
            
              (
              1
              )
            
          
        
        
          d
          R
        
      
    
    =
    0.
  




To determine the sheath potential at the meniscus, equate ion flux, assumed constant throughout the sheath,





  
    
      
        
          Γ
          i
        
        =
        
          
            
              n
              e
            
            ⁢
            
              v
              B
            
          
          
            R
            2
          
        
      
    
    
      
        (
        22
        )
      
    
  




to the electron flux at the meniscus





  
    
      
        
          Γ
          e
        
        =
        
          
            
              
                n
                s
              
              ⁡
              
                (
                
                  
                    v
                    _
                  
                  e
                
                )
              
            
            ⁢
            
              e
              
                
                  e
                  ⁢
                  
                    
                  
                  ⁢
                  
                    V
                    m
                  
                
                
                  kT
                  e
                
              
            
          
          
            4
            ⁢
            
              R
              2
            
          
        
      
    
    
      
        (
        23
        )
      
    
  




where





  
    
      v
      _
    
    e
  
  =
  
    
      (
      
        
          8
          ⁢
          
            
          
          ⁢
          
            eT
            e
          
        
        
          π
          ⁢
          
            
          
          ⁢
          m
        
      
      )
    
    
      1
      2
    
  




is the mean electron velocity and Vm is the potential of the meniscus with respect to the plasma/sheath edge. Thus, upon substitution of the Bohm velocity:





  
    
      
        
          
            
              n
              s
            
            ⁡
            
              (
              
                
                  eT
                  e
                
                M
              
              )
            
          
          
            1
            2
          
        
        =
        
          
            1
            4
          
          ⁢
          
            
              
                n
                x
              
              ⁡
              
                (
                
                  
                    8
                    ⁢
                    
                      
                    
                    ⁢
                    
                      eT
                      e
                    
                  
                  
                    π
                    ⁢
                    
                      
                    
                    ⁢
                    m
                  
                
                )
              
            
            
              1
              2
            
          
          ⁢
          
            e
            
              
                e
                ⁢
                
                  
                
                ⁢
                
                  V
                  m
                
              
              
                kT
                e
              
            
          
        
      
    
    
      
        (
        24
        )
      
    
  




which becomes





  
    
      
        
          V
          m
        
        =
        
          
            -
            
              T
              e
            
          
          ⁢
          
            
          
          ⁢
          
            
              ln
              ⁡
              
                (
                
                  M
                  
                    2
                    ⁢
                    π
                    ⁢
                    
                      
                    
                    ⁢
                    m
                  
                
                )
              
            
            
              1
              2
            
          
        
      
    
    
      
        (
        25
        )
      
    
  




This can be expressed in a more convenient form by substituting the mass of the extracted ion species. Krypton is a typical gas used in ion beam extraction from plasmas and has a mass of M=84au. Therefore, in this case, equation (25) can be rewritten as Vm=−5.05Te, which, assuming Te=3 eV, is approximately −15V. This now leads to suitable boundary conditions for equation (4):





  
    
      
        {
        
          
            
              
                
                  V
                  =
                  
                    0
                    ⁢
                    V
                  
                
                ,
                
                  r
                  =
                  
                    r
                    s
                  
                
              
            
          
          
            
              
                
                  V
                  =
                  
                    
                      -
                      15
                    
                    ⁢
                    
                      
                    
                    ⁢
                    V
                  
                
                ,
                
                  r
                  =
                  
                    r
                    m
                  
                
              
            
          
          
            
              
                
                  
                    
                      d
                      V
                    
                    /
                    
                      d
                      R
                    
                  
                  =
                  0
                
                ,
                
                  r
                  =
                  
                    r
                    s
                  
                
              
            
          
        
      
    
    
      
        (
        26
        )
      
    
  




where rm is the radius of curvature of the meniscus.

Potential Gradient at the Meniscus

Since both the voltage and its first derivatives are zero at rs, f is independent of the parameter p. From the familiar Langmuir-Blodgett relation:

α(γ)=γ−0.3γ2+0.075γ3−0.0143182γ4+0.0021609γ5−0.00026791γ6  (27)

This series expansion in conjunction with (6) and the boundary conditions now determine both the sheath width and the potential gradient at the meniscus edge. The sheath width is taken as the value of rm-rs for which (6) is equal to (25). The potential gradient at the meniscus edge is then equal to the first derivative of (6) taken at this value of rm.

Assuming a constant electron temperature and gas type, equations (6), (7) and (27) show that the potential gradient at the meniscus edge is dependent on the bulk plasma density and the radius of curvature of the meniscus.
Solving Poisson's Law Backwards
Consider the standard Langmuir-Blodgett problem in reverse, solving the differential equations from the exit of the extractor where the gradient is zero to the entry where it is large. This is tantamount to reversing the distribution in the case of a parallel beam or solving for the opposite convergence in the case of divergent or convergent beams. The f and g series as defined by equation (10) and equation (19) remain unchanged. However the definition of R and hence γ is altered. In the case of a diverging beam, R is taken to be the ratio of the current position to the first concentric surface and thus greater than unity. But in the case where Poisson's equation is solved backwards it is redefined to be the ratio of the current position to the second concentric surface and is hence less than one. The inverse is true for a convergent beam.
Presence of Electrons
To account for the presence of a population of electrons arriving from the tail of a Maxwellian distribution an exponential term is added to the original differential equations to account for the Boltzman relation.
Generalizing for the presence of electrons involves only the Right-Hand-Side (RHS) of Poisson's law:




  
    
      
        
          ρ
          
            ɛ
            o
          
        
        =
        
          
            e
            ⁡
            
              (
              
                
                  n
                  i
                
                -
                
                  n
                  e
                
              
              )
            
          
          
            ɛ
            o
          
        
      
    
    
      
        (
        28
        )
      
    
  




Since the electrons belong to a Maxwellian distribution their density as a function of potential is given by Boltzman's law:





  
    
      
        
          n
          e
        
        =
        
          
            n
            o
          
          ⁢
          
            exp
            ⁡
            
              (
              
                eV
                
                  kT
                  e
                
              
              )
            
          
        
      
    
    
      
        (
        29
        )
      
    
  




where no is some percentage ζ of the ion density in the beam aft of the extractor to take into account the fact that neutralization is not always 100%. The differential equations for spherical, cylindrical and planar symmetry then become:





  
    
      
        
          
            
              
                
                  r
                  2
                
                ⁢
                
                  V
                  ′′
                
              
              +
              
                2
                ⁢
                
                  rV
                  ′
                
              
            
            =
            
              0.6
              ⁢
              
                
              
              ⁢
              
                ne
                [
                
                  
                    
                      A
                      ~
                    
                    
                      V
                    
                  
                  +
                  
                    ζ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      exp
                      ⁡
                      
                        (
                        C
                        )
                      
                    
                  
                
                ]
              
            
          
          ,
          
            

          
          ⁢
          where
        
        ⁢
        
          

        
        ⁢
        
          
            A
            ~
          
          =
          
            
              v
              B
            
            ⁢
            A
            ⁢
            
              
                M
                
                  2
                  ⁢
                  e
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              
                rV
                ′′
              
              +
              
                V
                ′
              
            
            =
            
              0.6
              ⁢
              
                
              
              ⁢
              
                ne
                ⁡
                
                  [
                  
                    
                      
                        B
                        ~
                      
                      
                        V
                      
                    
                    +
                    
                      ζ
                      ⁢
                      
                        
                      
                      ⁢
                      
                        exp
                        ⁡
                        
                          (
                          C
                          )
                        
                      
                    
                  
                  ]
                
              
            
          
          ,
          
            

          
          ⁢
          where
        
        ⁢
        
          

        
        ⁢
        
          
            B
            ~
          
          =
          
            
              v
              B
            
            ⁢
            
              
                M
                
                  2
                  ⁢
                  e
                
              
            
          
        
      
    
    
      
        (
        30
        )
      
    
  




where n is the ion density in the beam, A is the area of the anode, and





  C
  =
  
    
      eV
      
        kT
        e
      
    
    .
  




Here the solutions take a very different form to those presented previously because of the exponential term. However, it is still possible to give a solution of the potential in terms of a MacLauren series, which when suitably normalized yields the correct form for implementation in the solutions to Laplace's equation.

Non-Zero Initial Velocity

The problem is set out in the same fashion as in section 2.1 for negligible initial velocity, except that now the kinetic energy relation is written:




  
    
      
        
          
            
              1
              2
            
            ⁢
            
              Mv
              2
            
          
          -
          
            
              1
              2
            
            ⁢
            
              Mv
              o
              2
            
          
        
        =
        eV
      
    
    
      
        (
        31
        )
      
    
  




When rearranged this yields:





  
    
      
        v
        =
        
          
            
              
                
                  2
                  ⁢
                  e
                
                M
              
              ⁢
              V
            
            +
            
              v
              o
              2
            
          
        
      
    
    
      
        (
        32
        )
      
    
  




But since:





  
    
      
        
          v
          o
          2
        
        =
        
          
            
              
                2
                ⁢
                e
              
              M
            
            ⁢
            
              V
              exi
            
          
        
      
    
    
      
        (
        33
        )
      
    
  




where Vexi is the energy of the particles at the exit of the previous stage, the differential equations stipulated by Poisson's law become for spherical, cylindrical and planar symmetries:





  
    
      
        
          
            
              
                r
                2
              
              ⁢
              
                V
                ′′
              
            
            +
            
              2
              ⁢
              r
              ⁢
              
                
              
              ⁢
              
                V
                ′
              
            
          
          =
          
            A
            
              
                V
                +
                
                  V
                  exi
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            
              r
              ⁢
              
                
              
              ⁢
              
                V
                ′′
              
            
            +
            
              V
              ′
            
          
          =
          
            B
            
              
                V
                +
                
                  V
                  exi
                
              
            
          
        
      
    
    
      
        (
        34
        )
      
    
  




And the solutions to these equations become:





  
    
      
        
          V
          =
          
            
              
                
                  (
                  
                    
                      9
                      4
                    
                    ⁢
                    A
                  
                  )
                
                
                  2
                  3
                
              
              ⁢
              
                f
                
                  4
                  3
                
              
            
            -
            
              V
              exi
            
          
        
        ⁢
        
          

        
        ⁢
        
          V
          =
          
            
              
                
                  (
                  
                    
                      9
                      2
                    
                    ⁢
                    B
                  
                  )
                
                
                  2
                  3
                
              
              ⁢
              
                g
                
                  4
                  3
                
              
            
            -
            
              V
              exi
            
          
        
      
    
    
      
        (
        35
        )
      
    
  




Importantly, upon substituting equation (35) into equation (34), equations (9) and (14) remain the same. This means that the initial velocity serves only to translate the solution vertically. The f and g series can still be obtained by the various means set forth above.
Laplace's Equation
In a given region of the extractor according to this invention, determining the electrode geometry amounts to solving Laplace's equation subject to the potential along the beam edge. Since in the three cases of interest—plane, cylindrical and spherical geometry—there is strong symmetry, compact solutions can be obtained. These were presented by Radley in 1957 along with a complete and rigorous mathematical derivation. A cursory overview is provided in the following.
However, before doing so it is noted that in treating the instability issues of the solution, Radley remarks that in as much as small variations of the initial surface can produce large differences in the solution so, conversely, do relatively large variations in electrode shapes away from the beam surface produce only small variations in the form of the beam surface. This will have important implications in the section Electrode Design because electrodes will have to be curtailed to avoid overlap or break-down proximity,
Strip Beam
A strip beam can be thought of as an infinite plane diode in which the cathode is the plane x=0 and all charge in the region y>0 has been suppressed. To determine electrodes that would extract a beam of this sort a family of equipotentials in y>0 must be found such that the conditions in y<0 are unchanged. The basic case of the potential distribution in space-charge limited flow when suitably normalised yields:




  
    
      
        V
        =
        
          
            
              x
              
                4
                3
              
            
            ⁢
            
              
            
            ⁢
            in
            ⁢
            
              
            
            ⁢
            the
            ⁢
            
              
            
            ⁢
            region
            ⁢
            
              
            
            ⁢
            y
          
          ≤
          0
        
      
    
    
      
        (
        36
        )
      
    
  




Since assuming zero ion temperature, the trajectories are rectilinear and, perpendicular to the emission surface defined by the extraction slit, the transverse forces on the beam at the beam edge are zero, so that




  
    
      ∂
      V
    
    
      ∂
      y
    
  
  =
  0




there. Restating the conditions in polar coordinates, gives:





  
    
      
        
          
            
              
                
                  
                    V
                    =
                    
                      r
                      
                        4
                        3
                      
                    
                  
                
              
              
                
                  
                    
                      
                        ∂
                        V
                      
                      
                        ∂
                        θ
                      
                    
                    =
                    0
                  
                
              
            
            }
          
          ⁢
          on
          ⁢
          
            
          
          ⁢
          θ
        
        =
        0.
      
    
    
      
        (
        37
        )
      
    
  




The solution in Cartesian coordinates is:





  
    
      
        
          V
          ⁡
          
            (
            
              x
              ,
              y
            
            )
          
        
        =
        
          Re
          ⁢
          
            {
            
              z
              
                4
                3
              
            
            }
          
        
      
    
    
      
        (
        38
        )
      
    
  




where

z=x+iy=reiθ  (39)

so that the solution is given by the conformal mapping:





  
    
      
        =
        
          
            V
            +
            
              i
              ⁢
              
                
              
              ⁢
              U
            
          
          =
          
            z
            
              4
              3
            
          
        
      
    
    
      
        (
        40
        )
      
    
  




Wedge Beam

For a wedge beam it is assumed that the meniscus and 0V equipotential surfaces are concentric cylinders. The coordinate system is therefore chosen so that the origin is at the vertex of the wedge, and the beam surface lies on θ=0. Since the trajectories are rectilinear, lying along the lines θ=constant, δV/δθ=0 on θ=0. The meniscus is taken to be r=rm and the cathode to lie somewhere in the region r>rn in the case of divergent beams and r<rm in the case of convergent beams. In the following it is assumed that the ratio R=r/rm is greater than unity or that the wedge is divergent. The potential distribution along the beam edge, 5 when suitably normalized, is given by:




  
    
      
        V
        =
        
          
            R
            
              2
              3
            
          
          ⁢
          
            β
            
              4
              3
            
          
        
      
    
    
      
        (
        41
        )
      
    
  




where β is an infinite power series in the variable γ=ln(R). Equation (41) can be written in series form:





  
    
      
        V
        =
        
          
            
              R
              
                2
                3
              
            
            ⁡
            
              [
              
                
                  ∑
                  
                    n
                    =
                    1
                  
                  ∞
                
                ⁢
                
                  
                    a
                    n
                  
                  ⁢
                  
                    γ
                    n
                  
                
              
              ]
            
          
          
            4
            3
          
        
      
    
    
      
        (
        42
        )
      
    
  




Since an=1, this last expression may be expanded by the multinomial theorem to give:




  
    
      
        V
        =
        
          
            R
            
              2
              3
            
          
          ⁢
          
            γ
            
              4
              3
            
          
          ⁢
          
            
              ∑
              
                n
                =
                1
              
              ∞
            
            ⁢
            
              
                b
                n
              
              ⁢
              
                γ
                
                  n
                  -
                  1
                
              
            
          
        
      
    
    
      
        (
        43
        )
      
    
  




The coefficients, bn, will be presented and discussed below. The potential outside the beam is obtained by writing Reiθ for R in equation (41). Thus, γ is replaced by ω=γ+iθ, so that:





  
    
      
        V
        =
        
          Re
          ⁢
          
            {
            
              
                R
                
                  2
                  3
                
              
              ⁢
              
                e
                
                  
                    2
                    3
                  
                  ⁢
                  iθ
                
              
              ⁢
              
                ω
                
                  4
                  3
                
              
              ⁢
              
                
                  ∑
                  
                    n
                    =
                    1
                  
                  ∞
                
                ⁢
                
                  
                    b
                    n
                  
                  ⁢
                  
                    ω
                    
                      n
                      -
                      1
                    
                  
                
              
            
            }
          
        
      
    
    
      
        (
        44
        )
      
    
  




Cylindrical Beam

For this case cylindrical polar coordinates, (r,θ,z), are enlarged with the axis of the beam along r=0. By suitably normalizing the coordinates, the beam surface can be taken as r=1. Since the system has axial symmetry, none of the variables involves θ, so this coordinate may be neglected. The ion flow considered is a cylindrical section of an infinite plane diode so that:




  
    
      
        
          
            
              
                
                  
                    V
                    =
                    
                      z
                      
                        4
                        3
                      
                    
                  
                
              
              
                
                  
                    
                      
                        ∂
                        V
                      
                      
                        ∂
                        r
                      
                    
                    =
                    0
                  
                
              
            
            }
          
          ⁢
          on
          ⁢
          
            
          
          ⁢
          r
        
        =
        1
      
    
    
      
        (
        45
        )
      
    
  




Applying the analysis in Radley, the solution is:





  
    
      
        
          V
          ⁡
          
            (
            
              r
              ,
              z
            
            )
          
        
        =
        
          
            1
            2
          
          ⁢
          π
          ⁢
          
            
          
          ⁢
          
            1
            
              
                (
                
                  -
                  
                    7
                    3
                  
                
                )
              
              ⁢
              
                (
                
                  
                    exp
                    ⁡
                    
                      (
                      
                        -
                        
                          
                            14
                            ⁢
                            i
                            ⁢
                            
                              
                            
                            ⁢
                            π
                          
                          3
                        
                      
                      )
                    
                  
                  -
                  1
                
                )
              
            
          
          ⁢
          
            
              ∫
              C
            
            ⁢
            
              
                
                  exp
                  ⁡
                  
                    (
                    
                      
                        -
                        p
                      
                      ⁢
                      
                        
                      
                      ⁢
                      z
                    
                    )
                  
                
                
                  p
                  
                    4
                    3
                  
                
              
              ⁢
              
                
              
              ⁢
              
                (
                
                  
                    
                      
                        J
                        1
                      
                      ⁡
                      
                        (
                        p
                        )
                      
                    
                    ⁢
                    
                      
                        Y
                        o
                      
                      ⁡
                      
                        (
                        pr
                        )
                      
                    
                  
                  -
                  
                    
                      
                        Y
                        1
                      
                      ⁡
                      
                        (
                        p
                        )
                      
                    
                    ⁢
                    
                      
                        J
                        o
                      
                      ⁡
                      
                        (
                        pr
                        )
                      
                    
                  
                
                )
              
              ⁢
              
                d
                p
              
            
          
        
      
    
    
      
        (
        46
        )
      
    
  




where p is a complex parameter, J and Y are Bessel functions of the first and second kind, and C is the contour defined as shown in FIG. 6.

In reality, the contribution of both straight line segments in C cancel and the contour reduces to a circle of radius ρ. In theory any value of ρ will work but it has been found that values of between 1 and 5 produced the most rapid and accurate results.
Conical Beam
Spherical polar coordinates (_r,θ,φ), are employed with a cone semi-angle of θ=θo, the anode on the sphere R=1 and the cathode in the region R>1 in the case of a diverging beam and R<1 in the case of a converging beam. Again, a diverging beam is assumed for this derivation. Writing γ=ln(R) and suitably normalizing the potential gives:




  
    
      
        V
        =
        
          f
          
            4
            3
          
        
      
    
    
      
        (
        47
        )
      
    
  




where





  
    
      
        f
        =
        
          
            ∑
            
              n
              =
              1
            
            ∞
          
          ⁢
          
            
          
          ⁢
          
            
              c
              n
            
            ⁢
            
              γ
              n
            
          
        
      
    
    
      
        (
        48
        )
      
    
  




The coefficients cn depend on the boundary conditions for the beam and are discussed herein. Again, this series representation is expanded using the multinomial theorem and becomes:





  
    
      
        V
        =
        
          
            γ
            
              4
              3
            
          
          ⁢
          
            
              ∑
              
                n
                =
                1
              
              ∞
            
            ⁢
            
              
            
            ⁢
            
              
                d
                n
              
              ⁢
              
                γ
                
                  n
                  -
                  1
                
              
            
          
        
      
    
    
      
        (
        49
        )
      
    
  




The boundary conditions to be applied are:





  
    
      
        
          
            
              
                
                  
                    V
                    =
                    
                      f
                      
                        4
                        3
                      
                    
                  
                
              
              
                
                  
                    
                      
                        ∂
                        V
                      
                      
                        ∂
                        r
                      
                    
                    =
                    0
                  
                
              
            
            }
          
          ⁢
          on
          ⁢
          
            
          
          ⁢
          θ
        
        =
        
          θ
          o
        
      
    
    
      
        (
        50
        )
      
    
  




and, following Radley, gives:





  
    
      
        
          V
          ⁡
          
            (
            
              r
              ,
              θ
            
            )
          
        
        =
        
          
            1
            2
          
          ⁢
          π
          ⁢
          
            
              
                d
                n
              
              ⁢
              
                sin
                ⁡
                
                  (
                  
                    θ
                    o
                  
                  )
                
              
            
            
              
                (
                
                  
                    -
                    n
                  
                  -
                  
                    4
                    3
                  
                
                )
              
              ⁢
              
                (
                
                  
                    exp
                    ⁡
                    
                      (
                      
                        -
                        
                          
                            14
                            ⁢
                            
                              
                            
                            ⁢
                            i
                            ⁢
                            
                              
                            
                            ⁢
                            π
                          
                          3
                        
                      
                      )
                    
                  
                  -
                  1
                
                )
              
            
          
          ⁢
          
            
              ∫
              C
            
            ⁢
            
              
                
                  exp
                  ⁡
                  
                    (
                    
                      -
                      vz
                    
                    )
                  
                
                
                  v
                  
                    n
                    +
                    
                      4
                      3
                    
                  
                
              
              ⁢
              
                (
                
                  
                    
                      
                        P
                        v
                      
                      ⁡
                      
                        (
                        μ
                        )
                      
                    
                    ⁢
                    
                      
                        Q
                        v
                        ′
                      
                      ⁡
                      
                        (
                        
                          μ
                          o
                        
                        )
                      
                    
                  
                  -
                  
                    
                      
                        Q
                        v
                      
                      ⁡
                      
                        (
                        μ
                        )
                      
                    
                    ⁢
                    
                      
                        P
                        v
                        ′
                      
                      ⁡
                      
                        (
                        
                          μ
                          o
                        
                        )
                      
                    
                  
                
                )
              
              ⁢
              
                
              
              ⁢
              
                d
                v
              
            
          
        
      
    
    
      
        (
        51
        )
      
    
  




where υ is a complex variable, P and Q are Legendre functions of the first and second kind and C is the contour defined in FIG. 6. Again, the straight line segments cancel so that the contour reduces to a circle. However, now, the radius ρ is dependent on R.

The first step in the design of electrodes is the choice of beam type, beam current and final extraction energy. For each choice there is a different electrode design. It must be determined from the first instance whether the beam is to be strip/wedge or cylindrical/conical and whether it is to be parallel, convergent or divergent. It must be kept in mind that the whole extractor is to act as if the beam was part of an entire diode from the plasma to the 0V equipotential. The electrode design is such that the extractor apertures follow the beam shape so that they are just in contact with the beam at their respective positions.
The Extraction Gaps
Stage 1
The extraction gap stage 1 comprises the aft face of the plasma electrode and front face of the acceleration electrode. Determining the beam distribution in this region, first takes into account that the electric field at the plasma/beam boundary is non-zero due to the plasma sheath. Then, given the beam form required, the relevant solution can be arrived at by following the derivations above. This in turn is substituted into the relevant solution of Laplace's equation. Note that in this region the potential distribution is convex and that the electric field at the exit of this region is very large,
Stage 2
The extraction gap stage 2 comprises the aft face of the acceleration electrode and the front face of the first Blocking electrode. The purpose of stage 2 is primarily to bring the electric field at the exit of stage 1 to zero. This is necessary as failure to do so will result in a strong ballooning of the electric field aft of the acceleration aperture. This will induce strong aberration and readjustment of the beam charge distribution leading to non-uniformity in the beam and potentially also reshaping of the meniscus away from the ideal plane, cylinder or sphere. The potential distribution in this region is obtained by solving Poisson's law backwards and assuming a non-negligible initial velocity. Again this expression is substituted into the relevant solution of Laplace's equation.
Laplace's equation, it seems, cannot be solved (in the Real domain) for a boundary condition that changes convexity because this would require that the equipotentials overlap, which in terms of electrodes means that they would need to occupy the same space. In fact, this is only a major problem for the acceleration electrode and a solution is to shape the electrode so that it is the median between the two ideal cases. It should be noted that close to the beam, the two equipotentials are almost identical and that away from the beam they are not strongly disparate.
Though the outer electrodes would eventually also overlap, and in practical terms would approach each other so that the inter-electrode gap would lead to break down, Radley has indicated that the effect of the electrodes away from the beam edge is increasingly negligible. Thus as a best approximation to an ideal extractor, the outer electrodes are made to extend to just outside break-down distance and the intermediate electrode is made to be the average of the two ideal equipotentials.
Blocking Electrodes
The Extraction electrodes have accelerated the beam just beyond the desired extraction energy and have done so in such a way that the electric field at the exit of stage 2 is 0 Volts per meter. In the case were positive ions are being extracted, this means that the voltage at the exit of the Extraction region stage 2 is some negative value. The purpose of the blocking electrodes is then to bring the beam potential back up to 0V while ensuring that the electric field upon exit is 0 Volts per meter. Assuming that this Blocking potential is sufficiently high, this impedes the flow of electrons from aft of the extractor to the Extraction region.
Since the potential gradient is zero upon entry to the Blocking electrode stage 1, the standard Langmuir-Blodgett representation can be used for the beam potential. For stage 2, however, the presence of electrons cannot be ignored since a population of higher energy electrons from the tail of the distribution will be able to penetrate some distance up the potential well. To solve for this distribution the relation given above is employed. This potential distribution is then substituted into the relevant solution of Laplace's equation.
Beam Neutralization
Beam neutralization is necessary aft of the extractor to compensate the considerable beam space charge. In many practical solutions, secondary electron emission from sputtering will be sufficient to provide a population of neutralizing electrons. However, if this were not the case, some electron source would be required such as a hollow cathode. In this situation emphasis should be given to ensuring that the electrons are as low energy as possible to avoid needing a large Blocking potential.
Plasma Density Range
There is a limit to how high the extraction energy can be taken as a function of density and beam form. For example, at a plasma density of 1014/cm3, no wedge beam form can be extracted without aberration because the minimum voltage per meter required for extraction is in excess of the breakdown limit of 107V/m. At 1013/cm3 extraction of all beam forms is possible, but for convergent beams the maximum extraction energy is capped at 20 kV in the stage 1 of the extraction region. Above this energy, the electrodes need to be too close for breakdown.
EXAMPLE
Diverging Wedge Beam
For this example a 5° diverging beam is extracted through a 1 mill-meter (mm) wide, 1 meter (m) long rectangular slit. The plasma source is assumed to be a Krypton plasma of density 1013 cm−3 at the sheath edge, The total extraction energy is 20 kV and the total blocking potential is 200V.
Extraction Gap Stage 1
The problem can be considered analogous to that of a complete cylindrical diode of curvature re=5.737 mm The current density flowing across the meniscus yields B=7.6218·108 according to equation (12). The ratio of the sheath to the pre-sheath is given by first solving equation (6) for the boundary conditions set in equation (25):




  
    
      
        
          β
          ⁡
          
            (
            
              γ
              m
            
            )
          
        
        =
        
          
            
              15
              
                3
                4
              
            
            
              
                
                  9
                  2
                
                ⁢
                
                  Br
                  m
                
              
            
          
          =
          
            1.71827
            ×
            
              10
              
                -
                3
              
            
          
        
      
    
    
      
        (
        52
        )
      
    
  




and then solving equation (19) to find γm=1.71797×10−3. From equation (7) R=1.00172 so that the sheath width is given by:

(1.00172)rs−rs=9.9 μm  (53)

The gradient at the sheath edge in terms of γm is




  
    
      
        
          
            
              d
              V
            
            
              d
              γ
            
          
          ⁢
          
            (
            
              
                γ
                m
              
              =
              
                1.71827
                ×
                
                  10
                  
                    -
                    3
                  
                
              
            
            )
          
        
        =
        
          11644
          ⁢
          
            
          
          ⁢
          V
          ⁢
          
            /
          
          ⁢
          unit
        
      
    
    
      
        (
        54
        )
      
    
  




Since





  
    γ
    =
    
      ln
      ⁡
      
        (
        
          r
          
            r
            s
          
        
        )
      
    
  
  ,
  
    
      d
      ⁢
      
        
      
      ⁢
      γ
    
    =
    
      dr
      r
    
  
  ,




which means:





  
    
      
        
          
            
              
                d
                V
              
              
                d
                γ
              
            
            ⁢
            
              (
              
                r
                m
              
              )
            
          
          ≈
          
            
              1
              
                5.737
                ×
                
                  10
                  
                    -
                    3
                  
                
              
            
            ⁢
            
              
                d
                V
              
              
                d
                γ
              
            
          
        
        =
        
          2.03
          ×
          
            10
            6
          
          ⁢
          
            
          
          ⁢
          V
          ⁢
          
            /
          
          ⁢
          m
        
      
    
    
      
        (
        55
        )
      
    
  




According to equation (18), p=1.43557×10−2, which is well below the limit of 0.15. Combining Table 2 with equation (20), the series expansion of g is found to be:

1.05846γ−0.521406γ2+0.209569γ3−0.0738154γ4+0.0165584γ5−0.00157496γ6  (56)

and hence an expression for the voltage in terms of γ through relation (13). Using the multinomial theorem, equation (56) can be written as:

1.0787−0.708501γ+0.342937γ2−0.140694γ3+0.0419959γ4−0.00842186γ5  (57)

Since the wedge is divergent, γ=ln(R), and we solve equation (44) with




  
    
      
        {
        
          
            
              
                R
                =
                
                  
                    
                      
                        x
                        ~
                      
                      2
                    
                    +
                    
                      
                        y
                        ~
                      
                      2
                    
                  
                
              
            
          
          
            
              
                θ
                =
                
                  
                    tan
                    
                      -
                      1
                    
                  
                  ⁡
                  
                    (
                    
                      
                        y
                        ~
                      
                      
                        x
                        ~
                      
                    
                    )
                  
                
              
            
          
        
      
    
    
      
        (
        58
        )
      
    
  




where the non-normalized coordinates are:





  
    
      
        {
        
          
            
              
                x
                =
                
                  
                    r
                    m
                  
                  ⁢
                  
                    x
                    ~
                  
                
              
            
          
          
            
              
                y
                =
                
                  
                    r
                    m
                  
                  ⁢
                  
                    y
                    ~
                  
                
              
            
          
        
      
    
    
      
        (
        59
        )
      
    
  




For the aft face of the plasma electrode we solve equation (44) for V=0. To solve for the front face of the acceleration electrode we solve for:





  
    
      
        V
        =
        
          
            
              10.1
              ⁢
              
                
              
              ⁢
              kV
            
            
              
                (
                
                  
                    9
                    2
                  
                  ⁢
                  
                    Br
                    m
                  
                
                )
              
              
                2
                3
              
            
          
          =
          0.138575
        
      
    
    
      
        (
        60
        )
      
    
  




which is a direct consequence of equation (41)

Extraction Gap Stage 2

At the entry to stage 2, Vext=10.2 kV. The solution to the g series is given by the standard Lanmuir-Blodgett relation but derived backwards. We substitute

g=γ+0.4γ2+0.0916667γ3+0.01424242γ4+0.001679275γ5+0.0001612219γ6  (61)

into equation (44), but with





  γ
  =
  
    -
    
      ln
      ⁡
      
        (
        
          r
          
            r
            d
          
        
        )
      
    
  




where rd is the radius of the concentric surface at the exit of the stage. At this point in the calculation rd is unknown, but by an iterative process a value can be found such that the potential gradient at ra at the entry to the stage is equal to that of the distribution at the exit of the previous stage. We define:





  
    
      
        
          
            γ
            a
          
          =
          
            -
            
              (
              
                
                  r
                  a
                
                
                  r
                  d
                
              
              )
            
          
        
        ⁢
        
          

        
        ⁢
        
          so
          ⁢
          
            
          
          ⁢
          that
        
      
    
    
      
        (
        62
        )
      
    
  
  
    
      
        
          
            e
            
              -
              
                γ
                a
              
            
          
          =
          
            
              r
              a
            
            
              r
              d
            
          
        
        ⁢
        
          

        
        ⁢
        
          and
          ⁢
          
            :
          
        
      
    
    
      
        (
        63
        )
      
    
  
  
    
      
        
          r
          d
        
        =
        
          
            r
            a
          
          ⁢
          
            e
            
              γ
              a
            
          
        
      
    
    
      
        (
        64
        )
      
    
  




Solving for γa:





  
    
      
        
          
            
              
                d
                V
              
              
                d
                γ
              
            
            ⁢
            
              (
              
                γ
                a
                -
              
              )
            
          
          ⁢
          
            |
            
              stage
              ⁢
              
                
              
              ⁢
              1
            
          
        
        =
        
          
            
              
                
                  d
                  V
                
                
                  d
                  γ
                
              
              ⁢
              
                (
                
                  γ
                  a
                  +
                
                )
              
            
            ⁢
            
              |
              
                stage
                ⁢
                
                  
                
                ⁢
                2
              
            
          
          =
          63586
        
      
    
    
      
        (
        65
        )
      
    
  




noting that the g series in stage 1 is different to that in stage 2. Thus γa=0.35464. The voltage at γa is V=0.239885 according to (41).

The deceleration 1 electrode front face is calculated by solving equation (44) for V=0 and the acceleration electrode aft face is calculated by solving equation (44) for V=0.239885. It should he noted that the acceleration voltage in this gap is 17484V which is significantly higher than 10.1 kV. To achieve a total acceleration energy of exactly 20 KV, an iterative approach will be required. But since this does not benefit the illustration of the method, this will not he done here. Because the particles have an initial energy of 10.1 kV, the basic solution is shifted up by this amount.
Blocking Gap Stage 1
A blocking voltage of 200V is required so it is assumed that the voltage in this stage is 100V. The standard Langmuir-Blodgett relation is solved to determine the gap size. This yields γ=0.00712381. Then to solve for this stage replace the standard Langmuir-Blodgett relation into equation (44). For the deceleration 1 aft face V=0 is solved and for the deceleration 2 front face V=0.00137203 is solved. The final solution is shifted up by 27584V.
Blocking Gap Stage 2
The solutions for the Poisson equation with the presence of Maxwellian electrons has not been resolved. For this reason this analysis does not assume the presence of electrons.
The gradient at the entry to the stage is 18729V/unit. Again the Langmuir-Blodgett relation considered backwards is used. The gap spacing that gives a gradient of 18729V/unit at the gap entry is γ=0.007196. The normalized voltage at this point is V=0.00138794. To solve for the Deceleration 2 aft face we solve equation (44) for V=0. To solve for the ground front face equation (44) is solved for V=0.00138794.
The electrodes defined in the previous sections assumed negative ions being extracted from 0 up to 27584V. The electrodes for the extraction of positive ions are identical except that the polarity of the field is reversed. In this scheme, ions are extracted from 27584V down to 0V.
To summarize, the process of determining the electrode shapes, positions and potentials to provide a minimum-aberration, high brightness beam involves the solution to a multiple boundary value problem. First, the designer selects the desired beam type to be implemented. As noted, the symmetrical beams considered include parallel, diverging and converging rectangular, cylindrical and spherical beams. The designer specifies the extraction energy. This gives the potential on the plasma electrode. The designer also specifies the total beam current and the beam current density. This information yields the plasma density, area of the plasma electrode aperture, and radius of curvature for the meniscus.
Then, one computes the electric field at the meniscus. Thus, the voltage and electric field at the first concentric surface is specified. Now the designer selects a first value for the potential at the second concentric surface. With these boundary conditions, one can determine the potential distribution between the first and second surfaces. From the potential distribution one can compute the electric field at the second surface. One also learns from applying the boundary conditions the radial distance between the first and second surface.
Next we impose the boundary condition that the electric field at the third surface is zero. We choose a relatively small voltage for the potential of the first blocking electrode. Using the first value selected for the potential at the second surface in conjunction with the specified boundary conditions at the third electrode, we determine a potential distribution between the second and third surfaces. From the potential distribution we can determine the electric field at the second surface. We compare this to the previously determined electric field at the second surface. To the extent the two values are not equal, the selected potential for the second surface must be adjusted. Thus, we iteratively select a potential for the second surface until we find the value that results in continuity of the electric field at the second surface. Note that the solution of the boundary conditions also gives the positions of the first three surfaces. Further, once the potential distributions are known, the shapes of the electrodes can be determined.
Note the differences between the three electrode system of the present invention and the three electrode system of the prior art. In the three electrode system of the prior art, the first stage is an acceleration stage and the second stage is a blocking stage. The electric field is not brought to a very low value at the aperture of the third electrode. In the present invention, the first stage is an acceleration stage that brings the electric field to a very high value at the second electrode. The second stage is a deceleration stage, which serves to bring the electric field to a very low value, preferably zero, at the aperture of the third electrode. In the absence of electrons, a three-electrode system, designed as described herein, is sufficient to produce a zero emittance growth beam.
Note that the electric field at the second surface is very large. It is several orders of magnitude greater than the electric field at the meniscus. Thus, even if the electric field is not completely brought to zero at the aperture of the third surface, a substantial reduction in the electric field there will result in a low emittance growth beam. For example, the electric field can be brought to less than 10 kV/m (kilo-Volts per meter) which is several orders of magnitude less than the electric field at the aperture of the second electrode and at least about an order of magnitude less than the field at the meniscus. Indeed, the electric field can brought to a value substantially smaller than the electric field at the meniscus. Clearly, the closer the electric field is to zero at the aperture of the third electrode, the lower the emittance growth. Thus, ideally, the electric field is brought substantially to zero. Simulations, using finite element software, predict that zero emittance growth can be achieved using the methods described herein.
In the presence of electrons, more than three electrodes are required. Thus, in one embodiment, two additional electrodes are provided to form a blocking region with two stages. The design of the blocking stage, which includes the aft face of the third electrode, is similar to the design of the first three electrodes just described. We start with the boundary conditions at the third electrode specified above. We select a potential for the fourth electrode and determine the potential distribution between the third and fourth surfaces. Then we compute the potential distribution between the fourth and fifth surfaces resulting from the boundary conditions of zero volts and zero electric field at the fifth surface. When the potential distributions produce equal electric fields at the fourth surface, the problem is solved, and a low-emittance beam is produced that exhibits an electric field magnitude at the fifth electrode aperture that is substantially less than the electric field magnitude at the meniscus. Ideally, the electric field magnitude at the fifth electrode aperture is zero.
Thus, the method of the present invention enables design of a three or five electrode system with low emittance growth. A four electrode system can also be designed. In the four electrode system, the first three electrode shapes, positions and potentials can be determined as described above. The fourth electrode can be designed as a blocking electrode forming a single blocking stage between the third and fourth surfaces. Indeed, using the methodology disclosed herein plasma extraction mechanisms can be designed with other numbers of electrodes and with other potential distributions. The method of the present invention is not restricted to ion plasma sources but is applicable to form a charged particle beam from a charged particle source generally.
Embodiments of the present invention can produce from a plasma ion source beams having brightnesses of greater than 105 A/sr/m2. The embodiments described above are designed to result in little or no emittance growth. Skilled persons will understand, that in some applications, the emittance may be less importance, and some emittance growth can be tolerated. In such cases, one can depart from the ideal teachings of the above embodiments without departing from the concepts of the invention. For example, where the ideal teaching produces zero electric field, a finite electric field that is substantially smaller than the electric field at the meniscus could be tolerated in some embodiments and still provide substantial benefit over the prior art. For example, an electric field of less than 5 percent, less than 10 percent, less than 15 percent, less than 20 percent, or even less than 30 percent or more of the maximum field at the meniscus could be tolerated in embodiments, and such fields could be considered to be substantially smaller than the field at the meniscus. Similarly, a field that is “substantially zero” is one in that is very small in relation to the meniscus field, and that results in very little emittance growth. Similarly, while it is ideal that the equipotential surface in the vicinity of the last electrode be parallel or concentric to the meniscus, some deviations from parallelism or concentricity can be tolerated in embodiments in which more emittance growth can be tolerated.
The electric field is preferably small enough, and the parallelism or concentricity preferably sufficient, to produce a beam having an emittance small enough to be focused to a submicron, or more preferably a sub-tenth-micron spot, and having sufficient current to be used to etch in applications such as circuit edit and mask repair.
Although the present invention and its advantages have been described in detail, it should be understood that various changes, substitutions and alterations can be made herein without departing from the spirit and scope of the invention as defined by the appended claims. The invention achieves multiple objectives and because the invention can be used in different applications for different purposes, not every embodiment falling within the scope of the attached claims will achieve every objective. Moreover, the scope of the present application is not intended to be limited to the particular embodiments of the process, machine, manufacture, composition of matter, means, methods and steps described in the specification. As one of ordinary skill in the art will readily appreciate from the disclosure of the present invention, processes, machines, manufacture, compositions of matter, means, methods, or steps, presently existing or later to be developed that perform substantially the same function or achieve substantially the same result as the corresponding embodiments described herein may be utilized according to the present invention. Accordingly, the appended claims are intended to include within their scope such processes, machines, manufacture, compositions of matter, means, methods, or steps.

## Claims

1. A low emittance growth extractor for extracting particles from a particle source to form a charged particle beam, comprising:
a first electrode closest to the particle source with an aperture through which particles can be extracted to form a beam, the particle source forming a meniscus in the vicinity of the aperture with a substantial electric field at the meniscus;
a second electrode having an electric potential that is less than the electric potential of the first electrode, wherein the beam flows through an aperture of the second electrode; and
a third electrode having an electric potential that is less than the electric potential of the second electrode, wherein the beam flows through an aperture of the third electrode.

2. The extractor of claim 1, in which the electric potential of the first electrode is greater than 20 kilovolts.

3. The extractor of claim 1, in which the electric potential of the second electrode is less than 20 kilovolts but greater than 10 kilovolts.

4. The extractor of claim 1, in which the electric potential of the third electrode is less than zero.

5. The extractor of claim 1, further comprising:
a fourth electrode having an electric potential that is greater than the third electrode, wherein the beam flows through an aperture of the fourth electrode; and
a fifth electrode having an electric potential that is greater than the fourth electrode, wherein the beam flows through an aperture of the fifth electrode.

6. The extractor of claim 5, in which the fourth electrode has an electric potential that is less than zero.

7. The extractor of claim 5, in which the fifth electrode has an electric potential that is substantially zero.

8. The extractor of claim 6, in which the fourth electrode has an electric potential that is less than zero.

9. The extractor of claim 1, in which the particle source comprises a plasma source.

10. A method for extracting particles from a particle source to form a low emittance growth charged particle beam, comprising:
extracting, through an aperture of a first electrode closest to the particle source, particles to form a beam, the particle source forming a meniscus in the vicinity of the aperture with a substantial electric field at the meniscus;
directing the beam to flow through an aperture of a second electrode having an electric potential that is less than the electric potential of the first electrode; and
directing the beam to flow through an aperture of a third electrode having an electric potential that is less than the electric potential of the second electrode.

11. The method of claim 10, in which the electric potential of the first electrode is greater than 20 kilovolts.

12. The method of claim 10, in which the electric potential of the second electrode is less than 20 kilovolts but greater than 10 kilovolts.

13. The method of claim 10, in which the electric potential of the third electrode is less than zero.

14. The method of claim 10, further comprising:
directing the beam to flow through an aperture of a fourth electrode having an electric potential that is greater than the third electrode; and
directing the beam to flow through an aperture of a fifth electrode having an electric potential that is greater than the fourth electrode.

15. The method of claim 14, in which the fourth electrode has an electric potential that is less than zero.

16. The method of claim 14, in which the fifth electrode has an electric potential that is substantially zero.

17. The method of claim 14, in which the fourth electrode has an electric potential that is less than zero.

18. The method of claim 14, in which the particle source comprises a plasma source.