# Patent

## Title
Nuclear batteries

## Abstract

We introduce a new technology for Manufactureable, High Power Density, High Volume Utilization Nuclear Batteries. Betavoltaic batteries are an excellent choice for battery applications which require long life, high power density, or the ability to operate in harsh environments. In order to optimize the performance of betavoltaic batteries for these applications or any other application, it is desirable to maximize the efficiency of beta particle energy conversion into power, while at the same time increasing the power density of an overall device. Various devices and methods to solve the current industry problems and limitations are presented here.

## Background
RELATED APPLICATIONS
This current application is a continuation of a application Ser. No. 13/042,444, filed Mar. 7, 2011 now U.S. Pat. No. 8,134,216, with the same title, inventors, assignee, and specification, which was recently allowed. Thus, this current application incorporates by reference all of the teachings and specification of its parent case, as also included here.
Ser. No. 13/042,444 in turn is a continuation-in-part of (and related to) U.S. application Ser. No. 12/888,521 filed Sep. 23, 2010 now U.S. Pat. No. 8,017,412, and Ser. No. 12/851,555, filed Aug. 6, 2010 now U.S. Pat. No. 8,487,392, which are based on the provisional applications 61/250,504, filed Oct. 10, 2009, 61/231,863, filed Aug. 6, 2009, and 61/306,541, filed Feb. 21, 2010, with common inventor(s), and same assignee (Widetronix Corporation). All of the above teachings are incorporated by reference here.


BACKGROUND OF THE INVENTION
We introduce a new technology for Manufactureable, High Power Density, High Volume Utilization Nuclear Batteries. Betavoltaic batteries are an excellent choice for battery applications which require long life, high power density, or the ability to operate in harsh environments. In order to optimize the performance of betavoltaic batteries for these applications or any other application, it is desirable to maximize the efficiency of beta particle energy conversion into power, while at the same time increasing the power density of an overall device. Increasing power density is a difficult problem because, while both the active area of the semiconductor used for the beta energy conversion and the layer of radioisotope that provides the betas for this conversion are very thin (100's of nanometers), the thickness of the substrate supporting the radioisotope layer and the overall thickness of the semiconductor device wafers are on the order of 100's of microns.
In another embodiment for this technology, there are several technical constraints that must be considered when designing a low cost, manufacturable, high volume, high power density silicon carbide (SiC) betavoltaic device. First, consideration must be given to the energy profile of radioisotopes to be used, and the volume at which such material can be produced. For example, tritium is one of the several viable radioisotope candidates, since it can be produced in sufficient quantities to support high volume device manufacture, and its energy profile fits well with a range of power generation design parameters.
Secondly, in order to produce high power density in betavoltaics, a large device surface area is required. There are issued and pending betavoltaic patents that mention patterning methods for pillars, pores or other structures which yield such high surface area—patent application Ser. No. 11/509,323 is an example, and can be used as a reference for pillared betavoltaic device construction. These methods must be optimized appropriately in order to meet fabrication objectives, while controlling costs.
Thirdly, SiC has been shown to be the ideal material for betavoltaic devices, e.g. see reference patent application Ser. No. 11/509,323. However, SiC has unique processing, fabrication and design requirements which must be met in order to produce a workable device. For example, fabrication of SiC devices requires high temperature epitaxial processes. Because of such high temperature requirements, these epitaxial processes add an element of complexity and cost, not seen with processes relating to other semiconductors, such as Si, and must be taken into account accordingly, or fabrication techniques must be developed to remove such complex and costly processes entirely.
Fourthly, it is desirable to integrate betavoltaic devices directly with Silicon (Si)-based electronics, including, but not limited to, microprocessor and memory devices. Thus, there is a need for designs and fabrication processes which anticipate such integration.
Devices which address or anticipate the aforementioned design considerations are disclosed in this current or co-pending applications, as mentioned above. Methods for fabricating same are also disclosed.
SUMMARY OF THE INVENTION
The small (submicron) thickness of the active volume of both the isotope layer and the semiconductor device is due to the short absorption length of beta electrons. The absorption length determines the self absorption of the beta particles in the radioisotope layer as well as the range, or travel distance, of the betas in the semiconductor converter which is typically a semiconductor device comprising at least one PN junction. We define a volume utilization factor, Volutilization, to quantitatively track how well a betavoltaic device is using the volume of the radioisotope source and the volume of the semiconductor converter (equation 1). To illustrate this, consider the simple betavoltaic structure shown in FIG. 1. There are three important length scales for optimization of such a device:
1) the self absorption length of the beta electrons in the radioisotope
2) the range of the beta electrons in the semiconductor converter material
3) the diffusion length of minority carriers in the semiconductor, Ldiff.
Ldiff determines the maximum thickness of any doped region (p-type or n-type) forming the PN junction. Note that although these design principles apply to any semiconductor material, including, but not limited to Si, GaAs, GaN, and diamond, herein, we focus on SiC because SiC has been shown to be the ideal material for a beta converter.
Also, this invention can be implemented using any beta emitting radioisotopes. Herein, we will consider the three isotopes Nickel-63 (N63), tritium (H3) and the tritides (Scandium Tritide, Titanium Tritide, etc.), and promethium-147 (Pm147). These isotopes have properties as listed in table 1. In this illustration for a simple structure shown in FIG. 1, the radioisotope is supplied by means of a foil. This foil could be carrying either N63, a tritiated metal such as scandium Tritide, or Pm147. We denote the range of the betas in SiC as LSiC and the self absorption length in the radioisotope as Lisotope. The volume utilization in this geometry, neglecting the contacts and isotope volume, is calculated as:




  
    
      
        
          Vol
          utilization
        
        =
        
          
            
              
                (
                
                  t
                  cell
                
                )
              
              ⁢
              Area
            
            
              
                (
                
                  
                    t
                    substrate
                  
                  +
                  
                    t
                    cell
                  
                
                )
              
              ⁢
              Area
            
          
          =
          
            
              (
              
                t
                cell
              
              )
            
            
              (
              
                
                  t
                  substrate
                
                +
                
                  t
                  cell
                
              
              )
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




Where
Area=the total device area, and
tsubstrate=the thickness of the SiC substrate
tcell=the thickness of the active SiC region.
Note that the value of Volutilization is between zero and one.
In order to maximize the power output, this planar style betavoltaic device has to be designed to capture as close to all of the beta electrons leaving the surface of the foil as possible. This means that tcell must be at least greater than the diffusion length of the minority carriers (tcell>Ldiff). However, any material thicker than this limit will not actively participate in energy conversion, so while tcell>Ldiff must be true, tcell must be as close as possible to Ldiff so as to maximize volume utilization. Further, the location of the PN junction depth from the surface of the device must be <Ldiff in order to collect the maximum number of electron hole-pairs.
In addition, one embodiment of this invention is a novel SiC betavoltaic device which comprises one or more “ultra shallow” P+N− SiC junctions and a pillared or planar device surface. Junctions are deemed “ultra shallow”, since the thin junction layer (which is proximal to the device's radioactive source) is only 300 nm to 5 nm thick. In one embodiment of this invention, tritium is used as a fuel source. In other embodiments, radioisotopes (such as Nickel-63, promethium or phosphorus-33) may be used. This is also addressed in our co-pending applications, mentioned above.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows schematic of beta voltaic converter, corresponding to FIG. 5.
FIGS. 2a-c show: Schematic illustration of one embodiment of the invention, corresponding to FIGS. 6a-c. The drawing shows a slap converter geometry being replaced by a number of cube-based converters.
FIG. 3 shows: Schematic of a beta voltaic device embodiment, corresponding to FIG. 7.
FIG. 4 shows a 3D representation, corresponding to FIG. 8. For clarity, space is inserted between the isotope vertical slabs. Ohmic contacts are formed in the rear of the device and on the devices bottom side.
FIG. 5 shows schematic of beta voltaic converter: green region is the SiC power converter, the blue region is the radio isotope, while the black regions are the ohmic contacts.
FIGS. 6a-c show: Schematic illustration of one embodiment of the invention. The drawing shows a slap converter geometry being replaced by a number of cube-based converters.
FIG. 7 shows: Schematic of a beta voltaic device embodiment: Green region is the SiC power converter, the blue region is the radio isotope, while the black regions are the ohmic contacts.
FIG. 8 shows a 3D representation. For clarity, space is inserted between the isotope vertical slabs. Ohmic contacts are formed in the rear of the device and on the devices bottom side and these contacts are shown in black.
FIG. 9 shows the diagram of n+-p−-n+ embodiment of the Endfire structure.
FIG. 10 shows drawing for n-p-n Comb Endfire device.
FIG. 11 shows: MOS capacitor formed on sidewall of the Endfire Betavoltaic device.
FIGS. 12 (a-b) shows: P-type MOS capacitor (a) with Vg=0, biased in the flatband mode (b) with Vg<0, biased in the accumulation mode.



DETAILED EMBODIMENTS OF THE INVENTION
Here are some embodiments of this invention:
In order to maximize the power output, this planar style betavoltaic device has to be designed to capture as close to all of the beta electrons leaving the surface of the foil as possible. This means that tell must be at least greater than the diffusion length of the minority carriers (tcell>Ldiff). However, any material thicker than this limit will not actively participate in energy conversion, so while tcell>Ldiff must be true, tcell must be as close as possible to Ldiff so as to maximize volume utilization. Further, the location of the PN junction depth from the surface of the device must be <Ldiff in order to collect the maximum number of electron hole-pairs.







TABLE I







β-emitting radioisotope and their


ranges in SiC and self absorption lengths










β-Emitting
Mean
Self absorption length
SiC absorption length


Isotopes
energy
(at mean beta energy)
(at mean beta energy)





N63
17.4 keV
0.67 μm
1.84 μm


Scandium
 5.6 keV
0.27 μm
0.25 μm


Trititide


Promethium
  67 keV
8.59 μm
19.56 μm 









Once the output power has been maximized, the only way to increase the power density is to reduce the thickness of the substrate by wafer polishing. A typical SiC wafer is about 350 microns, so if the thickness of the substrate was reduced to 50 microns, this would result in a seven times increase in power density.
The total power out of this planar betavoltaic device is given by:

PTotal=CtisotopeArea(SSSA)  (2)

If we take into account the substrate thickness tsubstrate, the power density produced by this geometry is given as:




  
    
      
        
          
            
              
                P
                Density
              
              =
              
                
                  P
                  Total
                
                
                  Total
                  ⁢
                  
                    
                  
                  ⁢
                  Device
                  ⁢
                  
                    
                  
                  ⁢
                  Volume
                
              
            
          
        
        
          
            
              =
              
                
                  
                    Ct
                    isotope
                  
                  ⁢
                  
                    Area
                    ⁡
                    
                      (
                      
                        S
                        SSA
                      
                      )
                    
                  
                
                
                  
                    (
                    
                      
                        t
                        substrate
                      
                      +
                      
                        t
                        cell
                      
                    
                    )
                  
                  ⁢
                  Area
                
              
            
          
        
        
          
            
              =
              
                
                  
                    Ct
                    isotope
                  
                  ⁢
                  
                    S
                    SSA
                  
                
                
                  (
                  
                    
                      t
                      substrate
                    
                    +
                    
                      t
                      cell
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




The conversion constant C takes into account the energy per beta electron the semiconductor loses (phonon, recombination etc.), the reflection of beta electrons at the semiconductor interface, the emission spectrum of the foil, and is directly related to the device efficiency. ‘Area’ is the area of the device as viewed from the top, and the thickness of the radioisotope is denoted by tisotope. SSSA is the specific surface activity, and is defined as the number of electrons per unit area which leaves the surface of the foil in the direction of the converter. This quantity is a measured value for a particular foil.
For a particular thickness, tisotope, of the radioisotope, only the betas that are not self absorbed leave the surface and are made available for harvesting by the SiC converter. This thickness of the radioisotope within which all the beta particles generated can leave the surface is called the self absorption length. The self absorption length of the beta particles with average energy is denoted by Lisotope. For the semiconductor, the range of penetration into the SiC of the beta particles with average energy is denoted by LSiC. Both LSiC and Lisotope are calculated from the following relationship.




  
    
      
        
          Range
          ⁡
          
            (
            
              in
              ⁢
              
                
              
              ⁢
              microns
            
            )
          
        
        =
        
          
            4
            
              100
              ⁢
              
                
              
              ⁢
              ρ
            
          
          ⁢
          
            
              E
              ⁡
              
                (
                
                  in
                  ⁢
                  
                    
                  
                  ⁢
                  keV
                
                )
              
            
            1.75
          
        
      
    
    
      
        (
        4
        )
      
    
  




where ρ is the density of either SiC or the radioisotope foil, and an expression for the ratio of the density of the two SiC to radioisotope can be written as:




  
    
      
        
          
            L
            SiC
          
          
            L
            isotope
          
        
        =
        
          
            ρ
            isotope
          
          
            ρ
            SiC
          
        
      
    
    
      
        (
        5
        )
      
    
  




An Embodiment
One embodiment of the invention is shown in FIG. 2. While the invention can be implemented with multiple junctions, this first embodiment will be described using a single junction. The top part of FIG. 2 shows the starting geometry which can be viewed as a combination of two slabs—a radioisotope slab and a SiC converter slab. The top slab (shown in red) is the radioisotope slab, and the bottom slab (shown in blue and yellow) is the PN junction slab. The top surface cross sectional dimensions (not shown) of the semiconductor slab are cellx and celly in the x and y directions respectively, and the z dimension (the thickness of the junction, also not shown) is denoted by tcell. In one example, we introduce additional isotope slabs to completely surround up to all four sides of the PN junction slab plus one isotope slab covering the junction slab's bottom or top surface or two additional slabs covering both the top and bottom junction surface. Multiple, and typically thousands, of these isotope enclosed semiconductor slabs will be fabricated across the wafer, resulting in a total top surface area of semiconductor slabs and isotope slabs equal to the final footprint of the new betavoltaic device. For comparison purposes, in this document, the total surface area of the high volume utilization betavoltaic design will approximate the original planar betavoltaic geometry area denoted as “Area” in the description of that planar device in the section above.
Note that there can be embodiments of this high volume utilization betavoltaic invention that use two isotope slabs, or three, or up to six isotope slabs, or e.g. the maximum number that can be physically added. For a given thickness of the junction, tcell, an increase in the number of isotope slabs will lead to an increase in the amount of beta electrons per unit volume available for harvesting by the betavoltaic, and therefore, an increase in the amount of power out for the overall total area of a device.
The relationship between the total area of the betavoltaic device and the cross sectional area, Acell, of the individual semiconductor slabs can be found by taking advantage of the square cross section of the slab design and creating a unit cell that includes both the semiconductor slab cross section and the isotope slabs surrounding it as shown in FIG. 2b. 
Then the area of the unit cell, Auc, is given by:

Auc=(cellx+2tisotope)(cellY+2tisotope)  (6)

For illustrative purposes, the semiconductor slab dimensions cell and cells shall be equal, however, in some embodiments of the invention this may not be the case. If cell and cells are equal, then:

cellx=cellY 

And Auc becomes:

Auc=(cellx+2tisotope)(cellx+2tisotope)

Auc=(cellx+2tisotope)2  (6b)

The total area, denoted as “Area”, covered by all the N unit cells on the device is equal to:

Area=N(cellx+2tisotope)2  (7)

And N, the number of cells in the active area of the device, can be found from:




  
    
      
        N
        =
        
          Area
          
            
              (
              
                
                  cell
                  x
                
                +
                
                  2
                  ⁢
                  
                    
                  
                  ⁢
                  
                    t
                    isotope
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          7
          ⁢
          a
        
        )
      
    
  




The values of each of the parameters defined above are determined by the material characteristics of both the isotope and the semiconductor. The following is a listing of the parameters and their determining material characteristics:
tcell: This parameter is determined by the minority carrier diffusion length, Ldiff, of the semiconductor material. It is important that all the electron hole pairs that are formed in the device active area can make it back to the junction. Keeping tcell close to Ldiff will ensure the maximum collection of electron-hole pairs. In some embodiments of the invention, the range for tcell can be 1 μm to 150 μm.
cellx: This parameter is determined by the range of the betas in the semiconductor, which means that it is also isotope dependent. Because there are isotope slabs on all four sides of the semiconductor slab in one or more embodiments of the invention, then for these embodiments, the cross section of the semiconductor slab can be substantially square to give equal range to the betas in all directions. In some of these embodiments of the invention, the range for cell can be 0.5 μm to 250 μm.
tisotope: This parameter is determined by the self absorption length, Lisotope, of the betas in their respective isotope sources. In one embodiment, tisotope is at least equal to Lisotope to ensure the most efficient volumetric use of the isotope slab. In some embodiments of the invention, the range for tisotope can be 0.1 μm to 20 μm.
One of the major differences between the planar betavoltaic design as well as designs which use textured active device areas with PN junctions that are conformal to a textured surface geometry, and this new high volume utilization betavoltaic invention is that certain surfaces/faces of as many as four isotope slabs are substantially perpendicular to one or more semiconductor slab PN junctions, thus, a significant amount of the betas whose energy are being harvested and used for power conversion enter the device in both the n-type and p-type regions within a diffusion length, Ldiff, of the junction(s). Using this configuration, we can significantly increase the number of betas per unit volume which can be harvested which will directly impact the total power output of the cell, as well as the power density.
To further illustrate the improvements of the invention over a planar device, we can calculate the relative power, PRel, of the new high volume utilization betavoltaic design relative to the standard planar betavoltaic design. The relative power is the ratio of the power of the high volume utilization geometry to the power of the planar single isotope slab geometry, or:




  
    
      
        
          P
          Rel
        
        =
        
          
            P
            
              multi
              -
              slab
            
          
          
            P
            planar
          
        
      
    
    
      
        (
        8
        )
      
    
  




The following are examples of Prel calculations for 6, 5 and 3 isotope slabs. As mentioned herein, other slab configurations in terms of slab quantity and position are possible.
The power for the high volume utilization betavoltaic invention with six isotope slabs, P6 slabs, is given by

P6 slabs={Ctisotope{[4cellxtcell]+[2(cellx)2]}SSSA}Nαedge2  (9)

Where αedge is an edge effect factor that adjusts for the intrinsic attenuation of the beta current from the isotope slabs around each individual SiC cell.
To calculate Prel we need the output power for the planar betavoltaic which was given in equation (2a) as:

PPlanar=CtisotopeArea(SSSA)  (2a)

Therefore,




  
    
      
        
          P
          
            Rel
            -
            
              6
              ⁢
              
                
              
              ⁢
              sides
            
          
        
        =
        
          
            
              P
              
                6
                ⁢
                
                  
                
                ⁢
                slabs
              
            
            
              P
              planar
            
          
          =
          
            
              
                [
                
                  
                    [
                    
                      4
                      ⁢
                      
                        
                      
                      ⁢
                      
                        cell
                        x
                      
                      ⁢
                      
                        t
                        cell
                      
                    
                    ]
                  
                  +
                  
                    2
                    ⁢
                    
                      
                        (
                        
                          cell
                          x
                        
                        )
                      
                      2
                    
                  
                
                ]
              
              ⁢
              N
              ⁢
              
                
              
              ⁢
              
                α
                edge
                2
              
            
            Area
          
        
      
    
    
      
        (
        10
        )
      
    
  




But from equation (7a) we know that:




  
    
      
        N
        =
        
          Area
          
            
              (
              
                
                  cell
                  x
                
                +
                
                  2
                  ⁢
                  
                    
                  
                  ⁢
                  
                    t
                    isotope
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          7
          ⁢
          a
        
        )
      
    
  




So substituting (7a) in (10), we get,




  
    
      
        
          P
          
            Rel
            -
            
              6
              ⁢
              
                
              
              ⁢
              sides
            
          
        
        =
        
          
            
              (
              
                
                  4
                  ⁢
                  
                    
                  
                  ⁢
                  
                    t
                    cell
                  
                  ⁢
                  
                    cell
                    x
                  
                
                +
                
                  2
                  ⁢
                  
                    
                      (
                      
                        cell
                        x
                      
                      )
                    
                    2
                  
                
              
              )
            
            ⁢
            
              (
              Area
              )
            
            ⁢
            
              α
              edge
              2
            
          
          
            
              Area
              ⁡
              
                (
                
                  
                    cell
                    x
                  
                  +
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    
                      t
                      isotope
                    
                  
                
                )
              
            
            2
          
        
      
    
    
      
        (
        
          10
          ⁢
          a
        
        )
      
    
  




Which gives,




  
    P
    
      Rel
      -
      
        6
        ⁢
        
          
        
        ⁢
        sides
      
    
  
  =
  
    
      
        
          (
          
            cell
            x
          
          )
        
        2
      
      ⁢
      
        (
        
          
            4
            ⁢
            
              
                t
                cell
              
              
                cell
                x
              
            
          
          +
          2
        
        )
      
      ⁢
      
        α
        edge
        2
      
    
    
      
        (
        
          
            cell
            x
          
          +
          
            2
            ⁢
            
              
            
            ⁢
            
              t
              isotope
            
          
        
        )
      
      2
    
  




And finally,




  
    
      
        
          P
          
            Rel
            -
            
              6
              ⁢
              
                
              
              ⁢
              sides
            
          
        
        =
        
          
            
              (
              
                
                  4
                  ⁢
                  
                    
                      t
                      cell
                    
                    
                      cell
                      x
                    
                  
                
                +
                2
              
              )
            
            ⁢
            
              α
              edge
              2
            
          
          
            
              (
              
                1
                +
                
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    
                      t
                      isotope
                    
                  
                  
                    cell
                    x
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        
          10
          ⁢
          
            
          
          ⁢
          aa
        
        )
      
    
  




If we only consider 5 radioisotope slabs, around the SiC cell (remove the bottom isotope), then the ratio for 5 is given by




  
    
      
        
          P
          
            Rel
            -
            
              6
              ⁢
              
                
              
              ⁢
              sides
            
          
        
        =
        
          
            
              (
              
                
                  4
                  ⁢
                  
                    
                      t
                      cell
                    
                    
                      cell
                      x
                    
                  
                
                +
                1
              
              )
            
            ⁢
            
              α
              edge
              2
            
          
          
            
              (
              
                1
                +
                
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    
                      t
                      isotope
                    
                  
                  
                    cell
                    x
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        11
        )
      
    
  




Similarly, for 3 isotope slabs (one on top, two on the sides) the ratio becomes




  
    
      
        
          P
          
            Rel
            -
            
              3
              ⁢
              
                
              
              ⁢
              sides
            
          
        
        =
        
          
            
              (
              
                
                  2
                  ⁢
                  
                    
                      t
                      cell
                    
                    
                      cell
                      x
                    
                  
                
                +
                1
              
              )
            
            ⁢
            
              α
              edge
            
          
          
            
              (
              
                1
                +
                
                  
                    2
                    ⁢
                    
                      
                    
                    ⁢
                    
                      t
                      isotope
                    
                  
                  
                    cell
                    x
                  
                
              
              )
            
            2
          
        
      
    
    
      
        (
        12
        )
      
    
  




The power density of the high volume utilization betavoltaic device is also an importance metric. The equation for the power density of a device with six isotope slabs, for example, is given by:




  
    P
    Density
  
  =
  
    
      {
      
        
          C
          ⁢
          
            {
            
              
                [
                
                  4
                  ⁢
                  
                    
                  
                  ⁢
                  
                    t
                    cell
                  
                  ⁢
                  
                    cell
                    x
                  
                
                ]
              
              +
              
                [
                
                  2
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      (
                      
                        cell
                        x
                      
                      )
                    
                    2
                  
                
                ]
              
            
            }
          
          ⁢
          
            S
            SSA
          
        
        
          
            (
            
              
                t
                substrate
              
              +
              
                t
                cell
              
            
            )
          
          ⁢
          Area
        
      
      }
    
    ⁢
    
      
        
          α
          edge
          2
        
        ⁢
        
          
        
        ⁢
        Area
      
      
        
          (
          
            
              cell
              x
            
            +
            
              2
              ⁢
              
                
              
              ⁢
              
                t
                isotope
              
            
          
          )
        
        2
      
    
  




Single Junction Ni63 Embodiment of Invention

The present invention may have embodiments as a single or multi junction device with either Ni63, tritium, or promethium-147, or other beta emitting isotopes. The following describes an embodiment of the invention which comprises a single junction with Ni63 used as the isotope source. This embodiment is shown in FIG. 3. In this case we have a single P/N junction surrounded by 3 slabs of radioisotopes shown in blue. The isotopes are electrically isolated from the P/N junction by a thin oxide layer (not shown). The N+ region is the SiC substrate.
FIG. 4 shows a 3D representation of this embodiment. For clarity, space is inserted between the adjacent radioisotope vertical slabs, where such space would normally be occupied by PN layers. Ohmic contacts are formed in the rear of the device and on the back of the substrate, and these contacts are shown in black.
Edge Effects and Design Equations
Typically, in designing a betavoltaic device, assumptions can be made regarding beta radiation traveling in a straight line with a density proportional to the specific activity. This is a good approximation for the planar case where the length of the foil is large compared to the absorption length in the SiC. However for the present invention, as one example, for each individual cell, one must take into account the edge effects for each mini cell. For a given beta energy and beta emitter position, the beta emitter will emit betas in all directions (all 360 degrees around). There will be an angle α which defines the edge effects. For angles less than 180 degrees there will be a loss of potential carriers given by α/180. We use the expression αedge in the above equations to represent the edge effects as a dimensionless quantity that takes into account carrier loss.
Fabrication of the High Volume Utilization Structure
One exemplary method for the fabrication of the high volume utilization betavoltaic invention is as follows:
1—Deep Silicon Carbide Etch:

    
    
        The channels for the vertical radioisotope slabs have to be etched first. This etch depth exposes the entire thickness of the active SiC cell to the radioisotope.
    
    


2—Oxide Passivation

    
    
        Thermal oxide will be grown on the SiC to serve as insulation from the shorting of the device junction on the sidewalls of the individual cells.
    
    


3—Amorphous Silicon Deposition

    
    
        A layer of amorphous Silicon (a-Si) will be blanket deposited over the deeply etched SiC wafer to allow for the re-planarization of the top surface.
    
    


4—CMP Planarization

    
    
        To ensure that lithography can be performed on the patterned surface of the SiC sample after etching, the a-Si deposited on the sample in the previous step has to be planarized. This planarization step provides a flat template for the subsequent photoresist and lithographic processes.
    
    


5—Wet Oxide Etch

    
    
        A wet oxide etch is done to remove any residual oxide that might be on the surface of the SiC before the metals for the ohmic contact are deposited. The presence of oxide would compromise the quality of the ohmic contact.
    
    


6—Ohmic Contact Metallization

    
    
        The metallization for the formation of ohmic contacts to p-type SiC is selectively deposited on the top surface of the SiC cells.
    
    


7—Reactive Ion Etch Removal of a-Si in Trenches

    
    
        The a-Si is removed from the surface of the device by Reactive Ion Etching (RIE)
    
    


8—Rapid Thermal Anneal

    
    
        The ohmic contact metallization deposited in step 6 is now annealed using a Rapid Thermal Annealer (RTA). This step forms low resistance contacts to the SiC devices.
    
    


9—Frontside Ni Blanket Metallization

    
    
        After the ohmic contacts are formed and annealed, a final blanket Nickel metallization will be done to connect all the individual SiC betavoltaic cells together and to serve as a seed layer for the eventual electroplated Nickel-63 radioisotope layer.
    
    


10—Backside Metallization

    
    
        The SiC betavoltaic device is a vertical device and as such may have an ohmic contact on the front and back of the device. This step forms the ohmic contact on the backside of the device.
    
    


Summary of Some of the Advantages of this Embodiment for Ni63 
We can summarize some of the advantages of this invention, as one embodiment:

    
    
        1. The VUtilization factor for this structure  ̃1 because all of the material is either emitting or collecting betas
        2. Because of the high volume utilization, the power density will increase
        3. This structure can efficiently allow for series combining of junctions to allow for a higher voltage output
        4. This structure allows for the deposition of Ni63 by electro chemistry because the “seed” layer for the deposition is at the bottom of the isotope channel and does not “shield” the beta emission.
        5. Unwanted beta emissions are easily shielded by the ohmic contacts that may be formed at the bottom of the structure along with, in some embodiments, an additional metal layer deposited on top of the structure.
    
    


Passivation of the Endfire Surface
The advantage of the Endfire betavoltaic concept is the increased area for beta particle input. Therefore, a larger source of energy is available for harvesting, relative to a planar betavoltaic device design. The disadvantage of this approach is that the increase in surface area comes with a potential introduction of surface charges and/or surface traps. Surface charges and/or surface traps can reduce the “effective minority lifetimes” of carriers in the device. The result of these charges is that carrier collection is reduced, which results in lower power output by the device.
Surfaces are literal terminations of crystal lattices and the dangling bonds that are formed as a consequence of this termination create localized energy states that can act as generation-recombination centers. These surface states have the potential to reduce the effective minority carrier lifetimes in devices. When the surface-to-volume ratio of a device increases, as is the case with going from a planar to the Endfire betavoltaic design, the total number of surface states increases, which can reduce the power output.
To mitigate this surface effect in the Endfire design, a novel metal-oxide-semiconductor (MOS) capacitor will be integrated with the betavoltaic device. The MOS device will be formed on the surface between the SiC device sidewalls, the insulating oxide, and the metal radioisotope source. This MOS capacitor will be biased in accumulation mode. (see FIGS. 11 and 12)
The MOS capacitor band diagram shown in FIG. 12(a) illustrates the flat band mode where there is no voltage bias on the metal terminal. This condition is characterized by the absence of band bending in the SiC and by the absence of charge build up at the surface. As a negative charge is introduced to the metal-semiconductor contact (FIG. 12(b)), an electric field is set up across the MOS capacitor. This field attracts the positively charged majority carriers in the p-type SiC to the surface where they quickly accumulate. This particular condition is called the accumulation mode. In the accumulation mode, the majority carrier density is increased at the surface and electric fields are produced which act to repel minority carriers from the surface. The action of the electric field on the minority carriers have the effect of isolating them from the traps. This electric field isolation allows for the Endfire design to be less susceptible to the effects of surface traps.
Biasing the MOS Capacitor:
The integrated MOS capacitor can be biased into accumulation by several sources including, but not limited to, the Endfire betavoltaic's generated voltage and the voltage from fixed oxide charges introduced during the fabrication of the devices.
Since the SiC Endfire betavoltaic will produce an open circuit voltage of 2 Volts, a portion of this voltage can be used to bias the MOS capacitor on the sidewalls. Fixed negative charge can also be implanted into the oxide to permanently bias the MOS capacitor into accumulation. The fixed negative charge will allow the device to remain in accumulation, regardless of the external resistive loads that the device may be connected to and will also simplify the fabrication process of the device, by eliminating the need to connect the negative output of the betavoltaic to the MOS terminal.
Alternate Embodiment of the Endfire Design
The Endfire betavoltaic concept can be implemented in different p-n junction configurations. An alternate configuration is shown in FIG. 9. Rather than just being a simple mini p-n junction slab (as the embodiment shown in FIG. 10), there are two back to back p-n junctions in parallel, built into the device, and both harvest beta energy to contribute to the total power output. The structure can be n+-p−-n+ (as shown in FIG. 10), or the mirror structure of p+-n−-p+. The advantages of this embodiment of the device are as follows:

    
    
        For the n+-p−-n+ structure, the minority carrier lifetimes are larger in p-type material
        The maximum depth of the device can be increased
        The total power output is higher
        Surface passivation is easier to achieve
    
    


In summary, we have the following figures: FIG. 1 shows schematic of beta voltaic converter, corresponding to FIG. 5. FIGS. 2a-c show: Schematic illustration of one embodiment of the invention, corresponding to FIGS. 6a-c. The drawing shows a slap converter geometry being replaced by a number of cube-based converters. FIG. 3 shows: Schematic of a beta voltaic device embodiment, corresponding to FIG. 7. FIG. 4 shows a 3D representation, corresponding to FIG. 8. For clarity, space is inserted between the isotope vertical slabs. Ohmic contacts are formed in the rear of the device and on the devices bottom side.
FIG. 5 shows schematic of beta voltaic converter: green region is the SiC power converter, the blue region is the radio isotope, while the black regions are the ohmic contacts. FIGS. 6a-c show: Schematic illustration of one embodiment of the invention. The drawing shows a slap converter geometry being replaced by a number of cube-based converters. FIG. 7 shows: Schematic of a beta voltaic device embodiment: Green region is the SiC power converter, the blue region is the radio isotope, while the black regions are the ohmic contacts.
FIG. 8 shows a 3D representation. For clarity, space is inserted between the isotope vertical slabs. Ohmic contacts are formed in the rear of the device and on the devices bottom side and these contacts are shown in black.
FIG. 9 shows the diagram of n+-p−-n+ embodiment of the Endfire structure. FIG. 10 shows drawing for n-p-n Comb Endfire device. FIG. 11 shows: MOS capacitor formed on sidewall of the Endfire Betavoltaic device. FIG. 12 shows: P-type MOS capacitor (a) with Vg=0, biased in the flatband mode (b) with Vg<0, biased in the accumulation mode.
Maximizing Charge Collection in SiC Betavoltaics—Influence of Junction Depth
This is also addressed in our co-pending applications, mentioned above: To quantify the extent of the surface, it is necessary to know the penetration depth, or range, RB in μm, of the beta electron in the semiconductor, which is given as:

RB (μm)=[4×E01.75 (keV)/100]/ρ(g/cm3)  (1e),

where E0 is the incident beta energy in keV, and ρ is the density of the semiconductor in g/cm3. The penetration depth is simply a function of the energy spectrum of the β-radiation, which is known. The spectrum, to first order, is given by

f(E0)=K√{square root over (E02+2mc2E0)}(E0(max)−E0)2  (2e)

where f(E) is the energy distribution function, m the electronic mass, c the speed of light, and K a normalization constant, such that we have:





  
    
      
        
          
            
              ∫
              0
            
            
              
                E
                0
              
              ⁡
              
                (
                max
                )
              
            
          
          ⁢
          
            
              f
              ⁡
              
                (
                
                  E
                  0
                
                )
              
            
            ⁢
            
              d
              
                E
                0
              
            
          
        
        =
        1
      
    
    
      
        (
        
          3
          ⁢
          e
        
        )
      
    
  




The energy extends to a maximum, E0(max), that typically lies at  ̃3 times the mean energy. For a given beta emitting isotope, a single E0(max) completely specifies the spectrum, as eq. 2e indicates. There is a Coulombic penetration factor that modifies equation (2e) above. This factor accounts for electrons being retarded by the Coulombic attraction from the nucleus, which skews the spectrum towards lower energies. Considering this factor, equation (2e) becomes:

f(E0)=KF(ZD,E0)√{square root over (E02+2mc2E)}(E0(max)−E0)2  (4e)

where F(ZD,E0), called the Fermi function, takes into account the Coulombic penetration effects. This function is tabulated in relevant semiconductor literature, and is related to the daughter nucleus atomic number, ZD, and the energy of the emitted β particle, E0. It can be approximated by:





  
    
      
        
          
            F
            ⁡
            
              (
              
                
                  Z
                  D
                
                ,
                
                  E
                  0
                
              
              )
            
          
          =
          
            
              2
              ⁢
              
                
              
              ⁢
              π
              ⁢
              
                
              
              ⁢
              v
            
            
              1
              -
              
                exp
                ⁡
                
                  (
                  
                    
                      -
                      2
                    
                    ⁢
                    
                      
                    
                    ⁢
                    π
                    ⁢
                    
                      
                    
                    ⁢
                    v
                  
                  )
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        where
        ⁢
        
          

        
        ⁢
        
          v
          =
          
            1.16
            ×
            
              10
              
                -
                3
              
            
            ⁢
            
              
                Z
                D
              
              /
              
                
                  
                    
                      E
                      0
                      2
                    
                    +
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      
                        mc
                        2
                      
                      ⁢
                      
                        E
                        0
                      
                    
                  
                  
                    
                      
                        m
                        2
                      
                      ⁢
                      
                        c
                        4
                      
                    
                    +
                    
                      E
                      0
                      2
                    
                    +
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      
                        mc
                        2
                      
                      ⁢
                      
                        E
                        0
                      
                    
                  
                
              
            
          
        
      
    
    
      
        (
        
          5
          ⁢
          e
        
        )
      
    
  




The penetration depth is then estimated as described in equation (1e). From (4e),  ̃65% of the spectrum energy lies at or below the mean, 5.5 keV for Tritium, while >80% of the energy lies below E(max)/2, which is  ̃9 keV for Tritium.
Assuming that all the beta-generated electron-holes beyond the surface junction p-type layer are collected, while none of those generated in the surface junction layer are collected, we can estimate the charge collection as a function of energy, or as simply the fraction of the total path length (RB) that lies beyond the junction region (Xj). This fraction at each energy in the beta spectrum is (RB−Xj)/RB. Integrating the total charge collection function, we obtain the total charge collection efficiency. More details and results are given in our co-pending applications, mentioned above, which are incorporated by reference here.
Any variations of the teachings above are also meant to be covered and protected by this current application.

## Claims

1. A nuclear battery device, said nuclear battery device comprising:
a P-N semiconductor junction, located between a P-type semiconductor layer and an N-type semiconductor layer;
two or more contacts;
two isotope foils;
wherein said two or more contacts are connected to said P-type semiconductor layer and said N-type semiconductor layer; and
wherein said P-type semiconductor layer and said N-type semiconductor layer are sandwiched between said two isotope foils;
wherein each of said two isotope foils are covered by an oxide layer;
wherein each of said oxide layer is connected to a metal contact;
a metal-oxide-semiconductor capacitor;
wherein said metal-oxide-semiconductor capacitor comprises said sandwiched said P-type semiconductor layer and said N-type semiconductor layer, surrounded on each side by said two isotope foils, surrounded on each side by each of said oxide layer, which is connected on each side by each of said metal contact.

2. The nuclear battery device as recited in claim 1, wherein said metal-oxide-semiconductor capacitor is biased in accumulation mode.

3. The nuclear battery device as recited in claim 1, wherein surface charges and surface traps are passivated.

4. The nuclear battery device as recited in claim 1, wherein power output of said nuclear battery device is increased.

5. The nuclear battery device as recited in claim 1, wherein surface dangling bonds, surface localized energy states, or surface generation-recombination centers are reduced.

6. The nuclear battery device as recited in claim 1, wherein effective minority carrier lifetimes are increased.

7. The nuclear battery device as recited in claim 1, wherein negative charge is introduced at metal-semiconductor contact.

8. The nuclear battery device as recited in claim 1, wherein an electric field is set up across said metal-oxide-semiconductor capacitor.

9. The nuclear battery device as recited in claim 1, wherein majority carrier density is increased at surface.

10. The nuclear battery device as recited in claim 1, wherein electric fields are produced which repel minority carriers from surface.

11. The nuclear battery device as recited in claim 1, wherein said metal-oxide-semiconductor capacitor is biased by betavoltaic's generated voltage.

12. The nuclear battery device as recited in claim 1, wherein said metal-oxide-semiconductor capacitor is biased by voltage from fixed oxide charges introduced during fabrication of said nuclear battery device.

13. The nuclear battery device as recited in claim 1, wherein fixed negative charge is implanted into oxide.

14. The nuclear battery device as recited in claim 1, wherein said metal-oxide-semiconductor capacitor is permanently biased into accumulation mode.

15. The nuclear battery device as recited in claim 1, said nuclear battery device comprising: an NPN structure.

16. The nuclear battery device as recited in claim 1, said nuclear battery device comprising: a PNP structure.

17. The nuclear battery device as recited in claim 1, wherein said P-type semiconductor layer and said N-type semiconductor layer are SiC semiconductor.

18. The nuclear battery device as recited in claim 1, wherein structure of said nuclear battery device comprises multiple junctions.

19. The nuclear battery device as recited in claim 1, wherein structure of said nuclear battery device comprises an amorphous layer.

20. The nuclear battery device as recited in claim 1, wherein structure of said nuclear battery device comprises at least one of the following: isotopes Nickel-63, Tritium, Scandium Tritide, Titanium Tritide, or Promethium-147.